Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,781,049 C→T intergenic (+32/‑166) lldD → / → trmL L‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0490CT97.2% 230.3 / NA 74intergenic (+32/‑166)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  major base T (43/27);  minor base G (0/2);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 1.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCG  >  NC_000913/3780803‑3781281
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
tACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGAGGGTGTAGCTCACCTGGTCAATATGATAGGAAAAGAGATGAAAGCGGGGATGACGCTGAATGGCGCGCAATAGATAAGCGAAATTACGCAAAATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGGCATGGCTCCCATGGCGAAAGGGAATGCGGCCTAGTCGTTTGCGCCCCCTCACCCCAACCCTCTACCTc                                                                                                                                                                                                                                          >  1:338513/1‑249 (MQ=255)
 aCTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAACAGAGATGAAAGTGGCGATGGCGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGCTTCGCTGGTGCAGGGGCTGGGTAAAAAGTTTCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAgg                                                                                                                                                                                                                                      >  1:288723/1‑251 (MQ=255)
        ggTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTCGCTAACCTGCTCAATCTGATCGAAAAAGCGATGAAAGGGGCGATGACGCTGACTGGCGCGAAAACGATTAGCGAAATTACGCAAGATTTGCTGGGGGAGGGGCTGGGTAAAGAGGTGCCTGCGGAACTGGGTCCCATGGCGAAAGGGGATTCGGCATAGTCGTTTGGGGCCCCTCACCCTTACCCTCTCCTTCAGGGGGAggg                                                                                                                                                                                                                               >  2:63396/1‑251 (MQ=255)
         gTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGGGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGGAAGTGGCGATGAAGCTGGCTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGGGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCAACGGCCCCCATGGCGGAAGGGAAAGCGGCATAGTAGTTTGCGCCCCCCCCCCCTAACCCTCTCCTTCCGGGGGAgggg                                                                                                                                                                                                                              >  2:450059/1‑251 (MQ=255)
            gTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGAcc                                                                                                                                                                                                                           <  2:337689/251‑1 (MQ=255)
            gTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGGGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAAATACGCAAGATTCGCTGGGGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCCCTGGCTCCCATGGCGGAAGGGAAAGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGAcc                                                                                                                                                                                                                           >  1:247234/1‑251 (MQ=255)
                tttCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGAGGAGCTAAACTGCTAAATCTGGTCGGAAAAGAGATGAAAGGGGCGATGACGCTGACTGGCGCGGAATCGATCAACGAAAATACGCAAGATTCGCCGGTGCAGGGGCTGGGGAAAGAGGTGCCTGCGGCACTGGCTCCCATGGCGGAAAGGGATGCGGCATAGTCGGTTTCGGCCCCCCACCCCAACCCCCCCCTTCCGGGAGAGGGGGCCggtc                                                                                                                                                                                                                       >  2:165471/1‑248 (MQ=255)
                      gTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                                                                                                 <  2:414419/251‑1 (MQ=255)
                      gTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCt                                                                                                                                                                                                                 <  1:447646/251‑1 (MQ=255)
                       tATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAggg                                                                                                                                                                                                                               >  1:128939/1‑236 (MQ=255)
                       tATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAggg                                                                                                                                                                                                                               <  2:128939/236‑1 (MQ=255)
                          gcgcTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGTCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCTCCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCtgtat                                                                                                                                                                                                             <  2:299432/251‑1 (MQ=255)
                             cTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACATGCTAAATCTGGTCGCAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGGAATCGATTAGCGAAATTACGCAAGCTTCGCTGGTGAAGGGGCTGGGGAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGGCTGCGGCATAGTCGTTTGCGCCCCCTCCACCTAACCCTCTCCTTCAGGGAGAGGGGGACGTTCGGCGCtggatgt                                                                                                                                                                                                           >  1:488778/1‑250 (MQ=255)
                                    cAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGGACTCCCTc                                                                                                                                                                                                   <  1:470093/251‑1 (MQ=255)
                                    cAGCGGGCCAGGCGGGTGTAGCTAACCTGCCTCATCTGATCGAAAAAGCGGTGAAAGTGGCGATGAAGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCCAGATTCGCTGGGGCAGGGGCTGGGTAAAAAGTTGCCTGCGGCCCTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCCAACCCTCTCCTTCCGGGGGAGGGGACCCTTCGGCGCTGTTTGTACTccccc                                                                                                                                                                                                   >  2:405093/1‑249 (MQ=255)
                                     aGCGGGCCAGGTGGGAGTAGCTAACCTGCTAAATCGGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGAGCAGCGAAATTTCGCAAGATTCGCTGGTGCAGGGGTTGGGTAAAGAGTTGCCTGCGTTATTGGCTCCCATGGCGAATGGGAATTCGGCAGAGTCGTTTGCGCCCCCTCACCCTAACCCGCGCCGTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                                                                                                  <  1:339522/251‑1 (MQ=255)
                                       ggggCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                                                                                                                <  1:115239/250‑1 (MQ=255)
                                       cGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGTACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                                                                                                                <  1:159089/251‑1 (MQ=255)
                                       cGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAct                                                                                                                                                                                                <  1:251825/251‑1 (MQ=255)
                                         ggCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCActct                                                                                                                                                                                              <  2:338513/251‑1 (MQ=255)
                                          gCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGTGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGCGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTg                                                                                                                                                                                             >  1:480683/1‑251 (MQ=255)
                                               gCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                                                                                                       <  1:468683/251‑1 (MQ=255)
                                                cGGGTGTAGCTAACCTGCTAACTCTGATCGAAAAAGAGATGAAAGTGGCGCTGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                                                                         >  2:269517/1‑248 (MQ=255)
                                                cGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGaaa                                                                                                                                                                                         <  1:269517/248‑1 (MQ=255)
                                                 gggTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTAGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGAAAGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGACAGGGAAAGCGGCATAGTCGGTTGCGCCCCCTCACACTCGCCCTCTCCATCATGGAGAGGGGACCGATAAGCGCTGCATGTCTTACCTCAATCTGAAAcca                                                                                                                                                                                      >  2:330474/1‑248 (MQ=255)
                                                  ggagtAGCTAACATGCTAAATCTGCTCGAAATAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATAGATCAGCGAAATTACGCAAGATTCGCTGGTGCCGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCATCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTACCTCACTCTGAAACGac                                                                                                                                                                                     >  2:468620/4‑250 (MQ=255)
                                                   gtgtAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                                                                                                   <  1:545140/251‑1 (MQ=255)
                                                       aGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCa                                                                                                                                                                               <  1:241124/251‑1 (MQ=255)
                                                          tAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                                                                                                            >  1:558038/1‑251 (MQ=255)
                                                           aaCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                                                                                                            <  1:405093/250‑1 (MQ=255)
                                                            aCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAg                                                                                                                                                                                                                                       >  2:76366/1‑191 (MQ=255)
                                                            aCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAg                                                                                                                                                                                                                                       <  1:76366/191‑1 (MQ=255)
                                                            aCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtt                                                                                                                                                                          <  1:408557/251‑1 (MQ=255)
                                                             ccTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                                                                                                         <  1:508521/251‑1 (MQ=255)
                                                                gCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTCtttt                                                                                                                                                                        <  2:31970/249‑1 (MQ=255)
                                                                 cTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                                                                                                     >  1:466863/1‑251 (MQ=255)
                                                                   aaaTCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTGACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                                                                                    <  1:456310/250‑1 (MQ=255)
                                                                         gATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                                                                                                               <  1:453804/249‑1 (MQ=255)
                                                                             gAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                                                                                                                                         >  2:190909/1‑251 (MQ=255)
                                                                              aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCt                                                                                                                                                        <  1:393793/251‑1 (MQ=255)
                                                                              aaaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCt                                                                                                                                                        >  1:400326/1‑251 (MQ=255)
                                                                                 aaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCg                                                                                                                                                     <  2:177323/251‑1 (MQ=255)
                                                                                              ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                        >  2:26728/1‑251 (MQ=255)
                                                                                                 gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggggg                                                                                                                                     >  2:420866/1‑250 (MQ=255)
                                                                                                 gATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGCTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCggggg                                                                                                                                     >  1:520094/1‑250 (MQ=255)
                                                                                                  tcgACGCCGACTGGGGCGAAATCGATCAGCGAAATTACGCAAGATTCGGTGGCGCAGGGGCTGGGTAAAGAGCTGCCTGCGGCACTGGCACCCATGCCGCACGGGAATGCGGCATAGCCGCTTGCGCCCCCTCACCCTAGCCCCCTCCTTCAGGGAGAGGGGACTGCTCGGCGCTGTATGTACTCCCTCACGTTGAACCGACAACGCATTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCGGGGGTg                                                                                                                                    >  2:262911/3‑251 (MQ=255)
                                                                                                       gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                 >  1:414880/1‑249 (MQ=255)
                                                                                                       gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCGCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                                                                 <  2:414880/249‑1 (MQ=255)
                                                                                                       gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCCCCCTTCAAGGGAGAGGGGCCGGTTCGCCCCTGTTTGTACTCCTTCACTTTGAAACGACACCGAACTCTTTTTTTCTCCCCCGCCCCTCCGGGGAGGGGGCCGGGGTgggggg                                                                                                                               >  1:527052/1‑251 (MQ=255)
                                                                                                            cgggCGGGAAAGCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAg                                                                                                                          <  2:505682/249‑1 (MQ=255)
                                                                                                                        cGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                                                              >  1:352784/1‑251 (MQ=255)
                                                                                                                        agATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTc                                                                                                              >  2:328829/2‑251 (MQ=255)
                                                                                                                                  aaTTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGTGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAAACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGAACTGCTCACCAATTTTc                                                                                                    >  1:438690/1‑251 (MQ=255)
                                                                                                                                          aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAg                                                                                                                                                                                                                                       <  1:6323/113‑1 (MQ=255)
                                                                                                                                          aaGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAg                                                                                                                                                                                                                                       >  2:6323/1‑113 (MQ=255)
                                                                                                                                           aGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCt                                                                                           >  2:444177/1‑251 (MQ=255)
                                                                                                                                              ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGCAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCCATTTTCTGCTATCCTGcc                                                                                        >  2:443220/1‑251 (MQ=255)
                                                                                                                                                       gCAGGGGCTGGGACAAGAATTGCATGCGGCAATGGCTCCCATGGAGAAAAGAAAAGCATCATAGTCGTTTGCGCCCCCTCACACTAACACGCTCCTTCAGGGAGAGGGGACCGTTAAACGCTGTATGTACTCCCTCACTATGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGCGAAGAGGGCCGGGGTGAGGGAAAAAGGCCGCACTACTCACCAAATTTCTGCTATCCTGCCCCCCCACt                                                                                >  1:517779/1‑250 (MQ=255)
                                                                                                                                                        ccggggCTGGGAAAAGAGTTGCCTTCGGCACTGGCTCCCATGGTGAAAGGGAATGCGGCATAGTCGTTTGAGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGCCCGTTCGGCGCTGTGTGCACTACCTCACTCTGAAAAGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCGGCTCACCAATTTGCTGCCATCCTGCCCCCGCaata                                                                               >  2:317765/3‑247 (MQ=255)
                                                                                                                                                                gggTAAAGAGTTGCCGGCGGCCCGGGCCCCAAGGGCGAAAGGTAAGGCGGCAGAGTGTTTTGCGCCCCCCCCCCCTACCCCTCCCCGACGGGGAGGGGGGCCCGTGCGGCGTTGTATGTCCTCCCCCACCTGGAAACGCCCCCGCACTTTTTTTTTCCCCCTCGCCCCCCCGGGGAGGGGGCGGGGGGGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                                                      <  2:400326/251‑1 (MQ=255)
                                                                                                                                                                    aaaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                                                                  >  1:408397/1‑251 (MQ=255)
                                                                                                                                                                      agagTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                                                                >  2:417125/1‑251 (MQ=255)
                                                                                                                                                                          ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCCGTATGCTAAAc                                                            >  2:66441/1‑251 (MQ=255)
                                                                                                                                                                                         ggCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg                                               >  1:141015/1‑250 (MQ=255)
                                                                                                                                                                                           cTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc                                            >  1:78049/1‑251 (MQ=255)
                                                                                                                                                                                            tCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCAACCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGACGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAACCATCGTACTTTACGAAcc                                           >  1:464771/1‑251 (MQ=255)
                                                                                                                                                                                                aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGa                                        >  2:445134/1‑250 (MQ=255)
                                                                                                                                                                                                 tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGaa                                       >  1:337350/1‑250 (MQ=255)
                                                                                                                                                                                                       aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTccc                                 >  1:41250/1‑250 (MQ=255)
                                                                                                                                                                                                                 cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAATTACTg                       >  1:282737/1‑251 (MQ=255)
                                                                                                                                                                                                                   gCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTCCTTTACGAACCAGAAATTCCGCCAAATACTGGc                     >  1:289687/1‑251 (MQ=255)
                                                                                                                                                                                                                                   ccccTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTccc                                                                                                                                                                >  1:239476/1‑96 (MQ=255)
                                                                                                                                                                                                                                   ccccTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACGCCCGCACTCTGAAACGACACCGCACTCCGTTGTTCTccc                                                                                                                                                                <  2:239476/96‑1 (MQ=255)
                                                                                                                                                                                                                                      cTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcg  >  1:472377/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                          
TACTGCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCG  >  NC_000913/3780803‑3781281

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