Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,871,776 C→T 5.6% I464I (ATC→ATT 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,7760CT 5.6% 342.4 / 5.4 107I464I (ATC→ATTyeaIputative membrane‑anchored diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base C (51/50);  new base T (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 = 3.44e-01

Only 100 of 108 total aligned reads displayed.

ACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTATTA  >  NC_000913/1871629‑1871921
                                                                                                                                                   |                                                                                                                                                 
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                                                                                                                                         >  2:1067942/1‑150 (MQ=255)
        gTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCa                                                                                                                                         >  1:474291/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)
                        cATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATaa                                                                                                                         <  2:1607878/150‑1 (MQ=255)
                             cGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCt                                                                                                                    <  2:81263/150‑1 (MQ=255)
                                    tATGAGTAAAAATGTCGAACTCCTGCCCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcga                                                                                                             <  2:1604703/150‑1 (MQ=255)
                                        cGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcg                                                                                                         >  1:1704945/1‑150 (MQ=255)
                                        cGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcg                                                                                                         <  1:589188/150‑1 (MQ=255)
                                        cGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACgcgaagcg                                                                                                         >  2:436641/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)
                                                   cGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAggg                                                                                              >  2:410050/1‑150 (MQ=255)
                                                    gAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggg                                                                                             >  1:223748/1‑150 (MQ=255)
                                                     aaCTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGCCCATCAGTAATGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAggggg                                                                                            <  1:338028/150‑1 (MQ=255)
                                                      aCTCCTGACAGGATTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTAATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAgggtgg                                                                                           <  2:335307/150‑1 (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:1250277/1‑150 (MQ=255)
                                                              ccGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAg                                                                                   <  2:566338/150‑1 (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)
                                                                   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)
                                                                                   atgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAg                                                              <  2:81050/150‑1 (MQ=255)
                                                                                    tgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACATGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGt                                                             >  1:482156/1‑150 (MQ=255)
                                                                                    tgTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGt                                                             >  2:262980/1‑150 (MQ=255)
                                                                                     gTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTg                                                            <  1:241996/150‑1 (MQ=255)
                                                                                          tGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCaaaa                                                       >  2:1190411/1‑150 (MQ=255)
                                                                                             aCAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAata                                                    >  1:1345043/1‑150 (MQ=255)
                                                                                             aCAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAata                                                    <  1:1630664/150‑1 (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                             <  1:1225294/150‑1 (MQ=255)
                                                                                                                    cGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTg                             <  2:1345043/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)
                                                                                                                             tccgtcatGACAAGCGAGATATTACGCGTGTCATGACGGAAGCAGATAAAGCATTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCa                    >  1:965414/8‑150 (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)
                                                                                                                                            gTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattacta     <  2:140574/150‑1 (MQ=255)
                                                                                                                                            gTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattacta     <  2:223748/150‑1 (MQ=255)
                                                                                                                                             tAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtataactat    <  2:1587686/150‑1 (MQ=255)
                                                                                                                                              aaTATCTTGCTTGTCATTACGGCATCATTTAAAGCCTTACGCGTATCGTAAAGCGTTGGGGGCAACAAAGTGTTTAATCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactatt   <  1:807807/150‑1 (MQ=255)
                                                                                                                                              aaTATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactatt   <  1:885055/150‑1 (MQ=255)
                                                                                                                                               atatCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTtattactatta  >  2:1128441/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACTCAATAGTGCCCAGGCAAAAATCATTGCGGAACGTATGCGTAAAAATGTCGAACTCCTGACCGGCTTTAGTAACAGATATGATGTTCCTGAACAAATGACCATCAGTATTGGCACGGTTTTTTCAACGGGTGACACGCGTAATATCTCGCTTGTCATGACGGAAGCAGATAAAGCCTTACGCGAAGCGAAAAGCGAGGGGGGCAACAAAGTGATTATTCATCATATTTAAGTGCAAAAATATTCAGAGCCATGCTTTTTGCGTGGCTTTTGCATACAATTTATTACTATTA  >  NC_000913/1871629‑1871921

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: