Predicted mutation
evidence seq id position mutation annotation gene description
RA 1 1,640,500 C→A intergenic (‑189/+144) unknown ← / ← unknown SSU rRNA ## 16S rRNA, small subunit ribosomal RNA/Protoporphyrinogen IX oxidase, oxygen‑independent, HemG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*11,640,5000CA100.0% 130.4 / NA 50intergenic (‑189/+144)unknown/unknownSSU rRNA ## 16S rRNA, small subunit ribosomal RNA/Protoporphyrinogen IX oxidase, oxygen‑independent, HemG
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (27/23);  total (27/23)

AAGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTCAAAAAGCTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATCTTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACA‑TTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGATCCGTTTCTCCGCCTGTGATACGCAT  >  1/1640260‑1640745
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aaGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTCAAAAAGCTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCg                                                                                                                                                                                                                                              <  1:235485/251‑1 (MQ=255)
          ccAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTCAAAAAGCTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAat                                                                                                                                                                                                                                    <  1:735344/251‑1 (MQ=255)
                                        ctTCAATTCAAAAAGCTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGccc                                                                                                                                                                                                      <  1:575543/251‑1 (MQ=255)
                                                       cTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATg                                                                                                                                                                                        <  1:760756/250‑1 (MQ=255)
                                                                    tAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGCACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGAt                                                                                                                                                                          >  1:760200/1‑251 (MQ=255)
                                                                                  aatCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGAt                                                                                                                                                            <  1:633800/251‑1 (MQ=255)
                                                                                  aatCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGAt                                                                                                                                                            <  1:642263/251‑1 (MQ=255)
                                                                                       attatAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATattatt                                                                                                                                                                                                                                                   <  1:341244/159‑1 (MQ=255)
                                                                                       attatAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGCAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAg                                                                                                                                                               >  1:450090/1‑243 (MQ=255)
                                                                                       attatAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGc                                                                                                                                                                                                        <  1:708114/202‑1 (MQ=255)
                                                                                       attatAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAg                                                                                                                                                                                                >  1:39879/1‑210 (MQ=255)
                                                                                                                        cACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTAt                                                                                                                                                                    <  1:516909/205‑1 (MQ=255)
                                                                                                                         aCCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGaaa                                                                                                                     <  1:259269/251‑1 (MQ=255)
                                                                                                                          ccTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATAttt                                                                                                                                                                                                                                <  1:835565/143‑1 (MQ=255)
                                                                                                                          ccTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGc                                                                                                                                                       <  1:536855/216‑1 (MQ=255)
                                                                                                                                  gTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGaaaa                                                                                                                    <  1:772768/243‑1 (MQ=255)
                                                                                                                                      aaTCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAacttac                                                                                                        <  1:819090/251‑1 (MQ=255)
                                                                                                                                       aTCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGAt                                                                                                                                                            <  1:831553/198‑1 (MQ=255)
                                                                                                                                       aTCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATAc                                                                                                               <  1:433639/243‑1 (MQ=255)
                                                                                                                                         cTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTAc                                                                                                                                             >  1:644177/1‑211 (MQ=255)
                                                                                                                                                           gAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGa                                                                                                                                                                           >  1:61331/1‑163 (MQ=255)
                                                                                                                                                           gAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACg                                                                                    >  1:617692/1‑250 (MQ=255)
                                                                                                                                                             aTCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAg                                                                                                                                                                              >  1:478721/1‑158 (MQ=255)
                                                                                                                                                             aTCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGa                                                                                                                                                                           <  1:117299/161‑1 (MQ=255)
                                                                                                                                                             aTCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGAt                                                                                                                                                                 <  1:631883/171‑1 (MQ=255)
                                                                                                                                                             aTCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCtt                                                                                 >  1:256672/1‑251 (MQ=255)
                                                                                                                                                              tCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGaaa                                                                                                                     >  1:287681/1‑214 (MQ=255)
                                                                                                                                                                tataAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCtt                                                                                 >  1:372647/1‑248 (MQ=255)
                                                                                                                                                                     gTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGaaaa                                                                          >  1:125125/1‑250 (MQ=255)
                                                                                                                                                                              caGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGAt                                                                                                                                                                 >  1:367852/1‑154 (MQ=255)
                                                                                                                                                                              caGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGAt                                                                                                                                                                 >  1:707015/1‑154 (MQ=255)
                                                                                                                                                                                gTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAg                                                                                                                                                                              <  1:37228/139‑1 (MQ=255)
                                                                                                                                                                                gTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGAt                                                                                                                                                                 <  1:486500/152‑1 (MQ=255)
                                                                                                                                                                                  ttGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGa                                                                                                                                                                           >  1:597852/1‑140 (MQ=255)
                                                                                                                                                                                  ttGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCg                                                                                                                                                      >  1:867275/1‑161 (MQ=255)
                                                                                                                                                                                    gTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGa                                                                                                                                                                           >  1:84126/1‑138 (MQ=255)
                                                                                                                                                                                        gattgattGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAACTATTTCGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAATATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAAt                                                      >  1:269609/1‑251 (MQ=255)
                                                                                                                                                                                                gTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAg                                                                                                                                                                                                <  1:799341/105‑1 (MQ=255)
                                                                                                                                                                                                         cATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAg                                                                                                                                                                              >  1:826168/1‑114 (MQ=255)
                                                                                                                                                                                                              gCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGAt                                                                                                                                                                 >  1:120194/1‑122 (MQ=255)
                                                                                                                                                                                                              gCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCAc                                                                 >  1:197924/1‑218 (MQ=255)
                                                                                                                                                                                                              gCATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGt                                >  1:463494/1‑251 (MQ=255)
                                                                                                                                                                                                               cATTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTa                               <  1:32343/251‑1 (MQ=255)
                                                                                                                                                                                                                aTTAAGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTaa                              <  1:853011/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACg                                                                        >  1:12986/1‑206 (MQ=255)
                                                                                                                                                                                                                   aaGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATAtt                                             >  1:232654/1‑233 (MQ=255)
                                                                                                                                                                                                                   aaGCAAGGCCGCAAATTATAGATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGAt                           >  1:654153/1‑251 (MQ=255)
                                                                                                                                                                                                                                       gATGAAGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGATCCGTTTCTCCGCCTGTGATa       >  1:24903/1‑251 (MQ=255)
                                                                                                                                                                                                                                           aaGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGATCCGTTTCTCCGCCTGTGATACGc    <  1:645083/250‑1 (MQ=255)
                                                                                                                                                                                                                                            aGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAacttactt                                                                                                      >  1:832462/1‑151 (MQ=255)
                                                                                                                                                                                                                                            aGATATTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACATTTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGATCCGTTTCTCCGCCTGTGATACGCat  >  1:839663/1‑251 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
AAGCCTGCCGCCAGCGTTCAATCTGAGCCATGATCAAACTCTTCAATTCAAAAAGCTTAATCCCTCAATAGCTGACTTAATTAATCTATTATAAATAATATAAATAAACGAATCTTCCGGCACCTATTAAGTTCAATCTTAAAATCATTTAAAACGAATCTATAAGTGCCCACACAGTTTGTATGATTGATTGTTAAAGAACATTTGCATTAAGCAAGGCCGCAAATTATAGATGAAGATCTTATTTGCCGTCAAGAAAATATTTTGTTTATTTTCCTTTAGTTACTGCCCTTGCAGTGCGGATGCGTATTATAGGGATCTCTATGATAGGATCAAGCGGATTTTTACA‑TTTTTTGTCGATAGCTTAGAAAAAATACAACTTACTTATTTCAACTGCAAAAACGCTTCTGAAAACGACCGCACTTTTTGCCAATCCGTATATTCGATTTCTTTTGTAAGATCCGTTTCTCCGCCTGTGATACGCAT  >  1/1640260‑1640745

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑