Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,964,480 G→A P226S (CCG→TCG)  thyA ← thymidylate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,964,4800GA98.0% 325.2 / ‑2.5 102P226S (CCG→TCG) thyAthymidylate synthase
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (51/49);  total (52/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 103 total aligned reads displayed.

TGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACG  >  NC_000913/2964343‑2964611
                                                                                                                                         |                                                                                                                                   
tGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGGCGAAGATGAATTAGGGTTTACGTTTGATAATCAACTTCGGCAGCGAAc                                                                                                                                   >  2:587504/1‑140 (MQ=255)
       aTGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTc                                                                                                                            >  2:1112868/1‑140 (MQ=255)
        tGTTTCGTAATTAGATAGCCACGGGCGCTTTAATGCCCGAATGCGGATGGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTcg                                                                                                                           <  2:586437/140‑1 (MQ=255)
          tttCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTcgcg                                                                                                                         >  1:1515856/1‑140 (MQ=255)
          tttCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTcgcg                                                                                                                         <  1:36844/140‑1 (MQ=255)
            tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGc                                                                                                                       >  1:292065/1‑140 (MQ=255)
            tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGc                                                                                                                       <  2:292065/140‑1 (MQ=255)
            tCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAACCAATTTGGGCAGCCAACGCGGTTCGCGGc                                                                                                                       >  1:576450/1‑140 (MQ=255)
               tAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTa                                                                                                                    >  2:1152344/1‑140 (MQ=255)
                 aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTcg                                                                                                                           <  1:1253302/131‑1 (MQ=255)
                 aTTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTcg                                                                                                                           >  2:1253302/1‑131 (MQ=255)
                  ttAGATAGCCACCGGCGATTTAATGCCCGGATGCGGATCTTAGCCTTCTATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAAtt                                                                                                                 <  2:576450/140‑1 (MQ=255)
                    aGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAAGCAACTTCGGCATCGTACGCGGATCGCGGCTTAATTGc                                                                                                               >  1:720885/1‑140 (MQ=255)
                       tAGCCACCGGCGCTTTAGTGCCCGGATGCGGTTCGTAGCTTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGa                                                                                                            <  2:662849/140‑1 (MQ=255)
                         gCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCtt                                                                                                                     >  2:43259/1‑129 (MQ=255)
                         gCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCtt                                                                                                                     <  1:43259/129‑1 (MQ=255)
                          ccACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGa                                                                                                         >  2:1006293/1‑140 (MQ=255)
                          ccACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGa                                                                                                         >  2:874816/1‑140 (MQ=255)
                          ccACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTGCGGCAGCGAACTCGGTTCGCGTCTTAATTGCAGATGa                                                                                                         >  1:410429/1‑140 (MQ=255)
                              cGGCGCTTTAATGCCCGGATGCGGATCGTGGCCTTCATTCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGttt                                                                                                     <  1:587504/140‑1 (MQ=255)
                               ggCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCCGCGAACGAGGTACGAGGCTTACTTGCAGATGAGTTTg                                                                                                    >  1:722896/1‑140 (MQ=255)
                               ggCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTg                                                                                                    >  2:1181129/1‑140 (MQ=255)
                               ggCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTg                                                                                                    >  1:993945/1‑140 (MQ=255)
                               ggCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTg                                                                                                    <  1:1181129/140‑1 (MQ=255)
                                gcgcTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTg                                                                                                                >  1:1356453/1‑127 (MQ=255)
                                gcgcTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTg                                                                                                                <  2:1356453/127‑1 (MQ=255)
                                  gcTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATc                                                                                                 <  2:1443064/140‑1 (MQ=255)
                                    tttAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGtt                                                                                                      >  1:952620/1‑133 (MQ=255)
                                    tttAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGtt                                                                                                      <  2:952620/133‑1 (MQ=255)
                                       aaTGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCatat                                                                                            <  1:1152344/140‑1 (MQ=255)
                                          gCCCGGATGCGGATCGTAGCCTTCAATCTCACAGTCTTCGAAACGGTAATCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGt                                                                                         <  2:548443/140‑1 (MQ=255)
                                          gCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGt                                                                                         <  2:1038921/140‑1 (MQ=255)
                                          gCCCGGATGCGGATCGTAGCCTCCAATCTCATAGTTTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGt                                                                                         <  2:653287/140‑1 (MQ=255)
                                            ccGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGATGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTg                                                                                       <  1:121226/140‑1 (MQ=255)
                                            ccGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTg                                                                                       <  2:1127002/140‑1 (MQ=255)
                                            ccGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTg                                                                                       <  2:1515856/140‑1 (MQ=255)
                                             cGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGc                                                                                      >  2:1011309/1‑140 (MQ=255)
                                                 tGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTa                                                                                  <  1:1384548/140‑1 (MQ=255)
                                                      aTCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGt                                                                                         >  2:1343863/1‑128 (MQ=255)
                                                      aTCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGt                                                                                         <  1:1343863/128‑1 (MQ=255)
                                                       tCGTAGCCTTCAATCTTCAAGTCTTCGAAACGGGACTCGCAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATg                                                                            <  2:754594/140‑1 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCCGATGAGTTTGAGCCATATGGTGGCTGTACAGATGCgt                                                                         >  1:754594/1‑140 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATgg                                                                                          >  1:905327/1‑123 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATgg                                                                                          <  2:905327/123‑1 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATg                                                                            >  1:863040/1‑137 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATg                                                                            <  2:863040/137‑1 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCgt                                                                         >  1:1151759/1‑140 (MQ=255)
                                                          tAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCAGCTTAATTGCAGTTGACTTTGATCCATATGATTGCTGTACAGATGCgt                                                                         >  2:493946/1‑140 (MQ=255)
                                                             ccTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTc                                                                      >  1:247711/1‑140 (MQ=255)
                                                              cTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGACTGTCATCCATATGGTTGATGTACAGATGctgtctc                                                                     >  2:673091/1‑134 (MQ=255)
                                                                tCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGcc                                                                   >  1:88200/1‑140 (MQ=255)
                                                                  aaTCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCAc                                                                 >  2:837815/1‑140 (MQ=255)
                                                                  aaTCTCAAAGTCTTCGAAACGGTAGTCGAAGATGAATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCAc                                                                 >  2:771158/1‑140 (MQ=255)
                                                                   aTCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCAcc                                                                >  2:547171/1‑140 (MQ=255)
                                                                       cAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATTTGATTGCTGTGCAGATGCGTGTCGACATCGGTc                                                            >  2:737664/1‑140 (MQ=255)
                                                                         aaGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCa                                                          >  2:119336/1‑140 (MQ=255)
                                                                             cTTCGAAACGGTAGTTGAAGATGGATTCGCGTTTACGCTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCATCGGTCCAGACa                                                      <  1:737664/140‑1 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTcc                                                           <  2:1182475/134‑1 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTcc                                                           >  1:1182475/1‑134 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTc                                                                      >  1:291612/1‑123 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGt                                                             <  1:1361640/132‑1 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGt                                                             >  2:1361640/1‑132 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACaa                                                     >  1:1032478/1‑140 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACaa                                                     >  2:193732/1‑140 (MQ=255)
                                                                              ttCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGCTTGCTGTACAGATGAGTGTCGCCACCGGTCCAGACaa                                                     >  1:611317/1‑140 (MQ=255)
                                                                               tCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGGACAGATGCGTGTCGCCACCGGTCCAGACaaa                                                    >  2:551137/1‑140 (MQ=255)
                                                                                cGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACaaaa                                                   >  2:376615/1‑140 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  1:1047871/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  1:547171/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  2:88200/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  1:1015647/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  2:1032478/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  1:193732/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  1:119336/140‑1 (MQ=255)
                                                                                     cGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  2:247711/140‑1 (MQ=255)
                                                                                      ggTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAccc                                             >  1:930953/1‑140 (MQ=255)
                                                                                      ggTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              <  2:369948/139‑1 (MQ=255)
                                                                                      ggTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCAcc                                              >  1:369948/1‑139 (MQ=255)
                                                                                           tCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTc                                        >  2:1362319/1‑140 (MQ=255)
                                                                                                gATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGAt                                   >  2:364848/1‑140 (MQ=255)
                                                                                                 aTGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCCGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATc                                  >  1:528688/1‑140 (MQ=255)
                                                                                                  tGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACaaaa                                                   >  1:377591/1‑122 (MQ=255)
                                                                                                  tGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACaaaa                                                   <  2:377591/122‑1 (MQ=255)
                                                                                                    gATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCa                               >  1:217607/1‑140 (MQ=255)
                                                                                                       tCGGGTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCActg                            <  1:1112868/140‑1 (MQ=255)
                                                                                                         gggTTTACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCActgct                          <  1:1362319/140‑1 (MQ=255)
                                                                                                              tACGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGcc                     >  2:851077/1‑140 (MQ=255)
                                                                                                               aCGTTTGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCca                    <  2:239959/140‑1 (MQ=255)
                                                                                                                    tGATAATCAACTTCGGCCGCGATCGCGGTTCTCGGTTTATTTGCAGATTAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCata               <  2:477824/140‑1 (MQ=255)
                                                                                                                    tGATAATCAACTTCGGCCGCGAGCGCGGTTCTCGGGTTACTTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCata               <  2:627266/140‑1 (MQ=255)
                                                                                                                    tGATAATCAACTTCGGCAGCGATCGCGGTTCGCGGCTTATTTGCAGATTAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCata               <  2:648292/140‑1 (MQ=255)
                                                                                                                    tGATAATCAACTTCGGCAGCGACCGCGGTTCCCGGATTAATTGCAGATAAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCata               <  2:410429/140‑1 (MQ=255)
                                                                                                                    tGATAATCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCata               <  1:771158/140‑1 (MQ=255)
                                                                                                                      aTAATCAATTTCGGCAGCGAACACGGTTCGTAGTTGAATCGCAGATGACTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATg             <  2:722896/140‑1 (MQ=255)
                                                                                                                      aTAATCAACTTCGGCAGCGAACCCGGTTCGCTGCTTAATCGCAGATGAATTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATg             <  2:611317/140‑1 (MQ=255)
                                                                                                                          tCAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCCGATGGCCCTGCTGGGCCATCATATGCAcc         >  1:566531/1‑140 (MQ=255)
                                                                                                                           cAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCa        >  2:24286/1‑140 (MQ=255)
                                                                                                                           cAACTTCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAAATCGCACTGCTGCCCCATCATATGCACCa        >  1:627266/1‑140 (MQ=255)
                                                                                                                                tCGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAGTCACCCACCTCCAGATAGCACTGCTGCGCCCTCATATGCACCAATAAc   >  2:559568/1‑140 (MQ=255)
                                                                                                                                 cGGCAGCGAACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCCCTGCTGAGCCATCGTATGCACCAATAAcg  >  1:739451/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                   
TGGCAGGATGTTTCGTAATTAGATAGCCACCGGCGCTTTAATGCCCGGATGCGGATCGTAGCCTTCAATCTCAAAGTCTTCGAAACGGTAGTCGAAGATGGATTCGGGTTTACGTTTGATAATCAACTTCGGCAGCGGACGCGGTTCGCGGCTTAATTGCAGATGAGTTTGATCCATATGGTTGCTGTACAGATGCGTGTCGCCACCGGTCCAGACAAAATCACCCACTTCCAGATCGCACTGCTGCGCCATCATATGCACCAATAACG  >  NC_000913/2964343‑2964611

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: