Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,871,770 T→A 5.6% R462R (CGT→CGA yeaI → putative membrane‑anchored diguanylate cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,871,7700TA 5.6% 312.5 / 5.9 107R462R (CGT→CGAyeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (51/50);  new base A (3/3);  total (54/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 109 total aligned reads displayed.

TTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTAT  >  NC_000913/1871621‑1871919
                                                                                                                                                     |                                                                                                                                                     
tttACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGt                                                                                                                                                       <  2:1581221/150‑1 (MQ=255)
  tACCGAACTCAATAGTGCCCAGGCAAAAATCATTTCGGAACTTATTCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCATCGGGTGACACGCGTaa                                                                                                                                                     <  2:1159895/150‑1 (MQ=255)
  tACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTCGTAACAGATATGATGATCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTaa                                                                                                                                                     <  2:721571/150‑1 (MQ=255)
    ccGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAata                                                                                                                                                   <  2:1529864/150‑1 (MQ=255)
        aCTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATAtctc                                                                                                                                               >  2:1630664/1‑150 (MQ=255)
             aTAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTg                                                                                                                                          >  1:1683595/1‑150 (MQ=255)
              tAGTGCTCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGt                                                                                                                                         >  2:1531679/1‑150 (MQ=255)
               aGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTc                                                                                                                                        <  1:1092537/150‑1 (MQ=255)
               aGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACAAGCGAGATATTACGCGTGTc                                                                                                                                        <  2:461289/150‑1 (MQ=255)
               aGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAAAAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTc                                                                                                                                        >  1:870524/1‑150 (MQ=255)
                gTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCa                                                                                                                                       >  1:474291/1‑150 (MQ=255)
                gTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCa                                                                                                                                       >  2:1067942/1‑150 (MQ=255)
                     cAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACAAGCGAGATATTACGCGTGTCacccgt                                                                                                                                  <  2:1140854/150‑6 (MQ=255)
                       ggCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGa                                                                                                                                <  2:638411/150‑1 (MQ=255)
                       ggCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACAAGCGAGATATTACGCGTGTCacccgttg                                                                                                                                <  1:1646913/150‑8 (MQ=255)
                       ggCAAAAATCATTGCGGAACGGATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGa                                                                                                                                <  1:432542/150‑1 (MQ=255)
                        gCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGaa                                                                                                                               >  1:647197/1‑150 (MQ=255)
                         cAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAg                                                                                                                              >  1:822385/1‑150 (MQ=255)
                           aaaaTCATTGCGGAACGAATGAGTAAAAATGTCGAACTCCAGACCGGCTTTATTAAAATATAAGATGTTCATGAACAAATGACCATCAGTAATGGCACTGTTTTTTCAACGGGTGACACGCGTAAAAACTCGCTTCTCATGACGGAAGCa                                                                                                                            >  2:1538161/1‑150 (MQ=255)
                                     cGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCt                                                                                                                  <  2:81263/150‑1 (MQ=255)
                                            tATGAGTAAAAATGTCGAACTCCTGCCCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcga                                                                                                           <  2:1604703/150‑1 (MQ=255)
                                                cGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcg                                                                                                       >  2:436641/1‑150 (MQ=255)
                                                cGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcg                                                                                                       >  1:1704945/1‑150 (MQ=255)
                                                 gTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcga                                                                                                      >  2:1285195/1‑150 (MQ=255)
                                                   aaaaaTGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTAcg                                                                                                              <  2:1292783/140‑1 (MQ=255)
                                                   aaaaaTGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTAcg                                                                                                              >  1:1292783/1‑140 (MQ=255)
                                                    aaaaTGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGaaaa                                                                                                   >  1:1598374/1‑150 (MQ=255)
                                                        tGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGa                                                                                               >  1:1037331/1‑150 (MQ=255)
                                                          tCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgg                                                                                             <  2:952082/150‑1 (MQ=255)
                                                            gAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggg                                                                                           >  1:223748/1‑150 (MQ=255)
                                                             aaCTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGCCCATCAGTAATGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAggggg                                                                                          <  1:338028/150‑1 (MQ=255)
                                                              aCTCCTGACAGGATTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTAATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggtgg                                                                                         <  2:335307/150‑1 (MQ=255)
                                                               cTCCTGACCGGCTTTAGTATCAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGc                                                                                        >  2:241996/1‑150 (MQ=255)
                                                               cTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGc                                                                                        <  1:436641/150‑1 (MQ=255)
                                                                tCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGca                                                                                       <  1:368069/150‑1 (MQ=255)
                                                                tCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGca                                                                                       <  1:616231/150‑1 (MQ=255)
                                                                      ccGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAg                                                                                 <  2:566338/150‑1 (MQ=255)
                                                                      ccGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAg                                                                                 >  2:1250277/1‑150 (MQ=255)
                                                                       cGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGt                                                                                <  1:1336632/150‑1 (MQ=255)
                                                                         gCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACCCGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGa                                                                              >  1:833344/1‑150 (MQ=255)
                                                                         ccTTTAGAAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGACGGGGGCAACAAAGTGa                                                                              >  2:242839/2‑150 (MQ=255)
                                                                          cTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGat                                                                             >  2:422235/1‑150 (MQ=255)
                                                                          cTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGat                                                                             <  1:1193731/150‑1 (MQ=255)
                                                                          cTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGat                                                                             <  1:1391554/150‑1 (MQ=255)
                                                                           tttAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGatt                                                                            >  1:1587686/1‑150 (MQ=255)
                                                                             tAGTTACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGattat                                                                          <  2:1683595/150‑1 (MQ=255)
                                                                              aGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGattatt                                                                         >  2:455226/1‑150 (MQ=255)
                                                                              aGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTCTTCAACGGGTGACACGCGTAATACCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGattatt                                                                         <  1:1121142/150‑1 (MQ=255)
                                                                              aGTAACAGATATGATGTTCCTGAAAAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACAAGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGAGAGGGGGGCAACAAAGAGattatt                                                                         >  1:89308/1‑150 (MQ=255)
                                                                                tAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATtca                                                                       <  2:1603262/150‑1 (MQ=255)
                                                                                   cAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACTCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATtcagca                                                                    >  1:1153231/1‑150 (MQ=255)
                                                                                    aGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTcatcat                                                                   <  2:870524/150‑1 (MQ=255)
                                                                                     gATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCata                                                                  >  2:1428707/1‑150 (MQ=255)
                                                                                         tgatgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTa                                                              >  1:975683/1‑150 (MQ=255)
                                                                                          gatgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTaa                                                             <  2:1037331/150‑1 (MQ=255)
                                                                                            tgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACATGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGt                                                           >  1:482156/1‑150 (MQ=255)
                                                                                            tgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGt                                                           >  2:262980/1‑150 (MQ=255)
                                                                                                  tGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaa                                                     >  2:1190411/1‑150 (MQ=255)
                                                                                                   gAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaaa                                                    >  1:751250/1‑150 (MQ=255)
                                                                                                     aCAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAata                                                  <  1:1630664/150‑1 (MQ=255)
                                                                                                     aCAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTAAAGTGCAAAAata                                                  >  2:807807/1‑150 (MQ=255)
                                                                                                      cAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGCGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAatat                                                 >  2:150597/1‑150 (MQ=255)
                                                                                                        aaTGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATAtta                                               >  2:181528/1‑149 (MQ=255)
                                                                                                             ccATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGGGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCTGAGc                                          <  2:822385/150‑1 (MQ=255)
                                                                                                             ccATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGc                                          >  2:1653385/1‑150 (MQ=255)
                                                                                                              cATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGcc                                         >  1:761025/1‑150 (MQ=255)
                                                                                                              cATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCACGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGcc                                         >  2:671719/1‑150 (MQ=255)
                                                                                                                 cAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATg                                      >  1:816428/1‑150 (MQ=255)
                                                                                                                  aGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGAAAAAATATTCAGAGCCATGc                                     >  2:869461/1‑150 (MQ=255)
                                                                                                                   gTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCt                                    >  2:885055/1‑150 (MQ=255)
                                                                                                                       tGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCttttt                                >  1:976637/1‑150 (MQ=255)
                                                                                                                        ggCACTGTTTTTTCAACTGGTTACACGCTTAATATCTCGCTTTTCATGACGGAATCATATAAAGCCATCCGCGAAGCGAAAAGCGAGGGGGGCTACAAAGTGATTATTCATAATATTTAAGTGAAAAAATTTTCAGATCCATGCTTTTTg                               <  1:1538161/150‑1 (MQ=255)
                                                                                                                        ggCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTg                               >  2:1546880/1‑150 (MQ=255)
                                                                                                                         gCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCGTATTTAAGTGCAATAATATTCAGAGCCATGCTTTTTGc                              >  2:1448473/1‑150 (MQ=255)
                                                                                                                          cACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCg                             <  1:1428707/150‑1 (MQ=255)
                                                                                                                            cGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTg                           <  2:1345043/150‑1 (MQ=255)
                                                                                                                            cGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTg                           <  1:1225294/150‑1 (MQ=255)
                                                                                                                             ggTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTg                               <  1:44092/145‑1 (MQ=255)
                                                                                                                             ggTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTg                               >  2:44092/1‑145 (MQ=255)
                                                                                                                             ggTTTTTACAACGGGTGACACGCGTAATATCTCTCTTGTAATGACGGAATCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTgg                          <  2:89308/150‑1 (MQ=255)
                                                                                                                              gTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGc                         <  2:474291/150‑1 (MQ=255)
                                                                                                                                tttttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCtt                       >  1:140574/1‑150 (MQ=255)
                                                                                                                                  tttCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCtttt                     >  1:396589/1‑150 (MQ=255)
                                                                                                                                   ttCTACGGGTGTCACGCGTAATATCTCGCTTGTCATGACTGAAGCAGATAAAGGCTTTCGCGAAGCGAAAAGCGATGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGTTTTTTGCGTGTCTTTTg                    <  1:1322047/150‑1 (MQ=255)
                                                                                                                                     cAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                  >  2:843876/1‑150 (MQ=255)
                                                                                                                                     cAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                  <  2:833344/150‑1 (MQ=255)
                                                                                                                                      aaCGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGACAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCAt                 <  1:1285195/150‑1 (MQ=255)
                                                                                                                                         tcatGACAAGCGAGATATTACGCGTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACa              >  1:855081/4‑150 (MQ=255)
                                                                                                                                         gtgtgtCACGCGTAATATCTCGTTTGTCATGACGGAAGCTTATAAAGCCTTACGGGAATCGAAAAGCGTGGGGGGCATCAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGGGTGGCTTTTGCATAct              <  1:242839/148‑2 (MQ=255)
                                                                                                                                         gggTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATAct              <  1:671719/150‑2 (MQ=255)
                                                                                                                                         gggTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACa              >  2:17054/1‑150 (MQ=255)
                                                                                                                                          catGACAAGCGAGATATTACGCGTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACaa             >  2:212910/3‑150 (MQ=255)
                                                                                                                                             gACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAATTGATTATTAATAATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAAttt          >  2:956462/1‑150 (MQ=255)
                                                                                                                                               cacGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttat        <  1:410050/150‑1 (MQ=255)
                                                                                                                                                acGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATTTTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatt       <  1:1190411/150‑1 (MQ=255)
                                                                                                                                                acGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatt       <  2:976637/150‑1 (MQ=255)
                                                                                                                                                 cgcgTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATttatta      <  1:1250277/150‑1 (MQ=255)
                                                                                                                                                   cgTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattact    <  1:262980/150‑1 (MQ=255)
                                                                                                                                                    gTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattacta   <  2:223748/150‑1 (MQ=255)
                                                                                                                                                     tAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtataactat  <  2:1587686/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTTACCGAACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTAT  >  NC_000913/1871621‑1871919

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