Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,261 G→A 8.0% intergenic (+23/‑124) yeaH → / → cdgI DUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,870,2610GA 8.0% 370.3 / 12.1 123intergenic (+23/‑124)yeaH/cdgIDUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base G (68/44);  new base A (8/2);  total (77/46)
Fisher's exact test for biased strand distribution p-value = 3.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 125 total aligned reads displayed.

GCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCAT  >  NC_000913/1870112‑1870410
                                                                                                                                                     |                                                                                                                                                     
gCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCg                                                                                                                                                       >  1:1336756/1‑150 (MQ=255)
  gagaATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGcc                                                                                                                                                     >  1:210332/1‑150 (MQ=255)
    gaATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTg                                                                                                                                                   >  2:1006731/1‑150 (MQ=255)
      atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGc                                                                                                                                                 <  1:599890/150‑1 (MQ=255)
          gAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCTGAt                                                                                                                                             <  2:924654/150‑1 (MQ=255)
                 tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                      <  2:383373/150‑1 (MQ=255)
                 tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                      <  1:87800/150‑1 (MQ=255)
                  gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGTCTAAAACTATCAGCCAGGTCATTATCGCCTCGCTGATTTTTAGCt                                                                                                                                     >  1:11215/1‑150 (MQ=255)
                       cTACTTTCGACAACTTTGCTATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTg                                                                                                                                <  1:947890/150‑1 (MQ=255)
                          cTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTaa                                                                                                                             >  1:189279/1‑150 (MQ=255)
                             tCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                          <  2:1326916/150‑1 (MQ=255)
                              cGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTa                                                                                                                         >  1:1603184/1‑150 (MQ=255)
                                    cTTTGCGATGCAGCACATCAGCGACCAGGATGATATTTATCAGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAt                                                                                                                   >  2:415164/1‑150 (MQ=255)
                                     tttGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAtt                                                                                                                  <  2:410254/150‑1 (MQ=255)
                                       tGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATtta                                                                                                                <  1:752536/150‑1 (MQ=255)
                                          gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                             <  2:1603184/150‑1 (MQ=255)
                                            tGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtaca                                                                                                           <  2:1232753/150‑1 (MQ=255)
                                            tGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtaca                                                                                                           <  1:641218/150‑1 (MQ=255)
                                              cagcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTAGTGTAATTTATCTCATTTTTtagata                                                                                                         <  2:854429/150‑1 (MQ=255)
                                              cagcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacata                                                                                                         <  2:394740/150‑1 (MQ=255)
                                              cagcaCATCCGCGACCAGGATGATAATTATCCGGTGTTCCGTGAACTGTATCAAAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacata                                                                                                         <  2:1445970/150‑1 (MQ=255)
                                                   caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                    >  2:1575680/1‑150 (MQ=255)
                                                     tCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGc                                                                                                  >  1:1507862/1‑150 (MQ=255)
                                                      ccGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCt                                                                                                 >  1:510222/1‑150 (MQ=255)
                                                        gcgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTa                                                                                               <  2:1237877/150‑1 (MQ=255)
                                                            ccAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAAt                                                                                           >  2:666848/1‑150 (MQ=255)
                                                               ggATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGc                                                                                        >  2:1573960/1‑150 (MQ=255)
                                                                   gatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTa                                                                                    >  1:576922/1‑150 (MQ=255)
                                                                       tttATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAata                                                                                <  2:350061/150‑1 (MQ=255)
                                                                         tATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                              >  2:1396761/1‑150 (MQ=255)
                                                                         tATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                              <  1:765920/150‑1 (MQ=255)
                                                                          aTCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGa                                                                             >  1:1058727/1‑150 (MQ=255)
                                                                          aTCCGATCTTCCTTTAACTGTTTCATAAACAAAATTCAACATCTAAATGCTAAAACTATCATCAAGATAATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGa                                                                             >  2:1478903/1‑150 (MQ=255)
                                                                           tCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAc                                                                            <  1:935905/150‑1 (MQ=255)
                                                                             cGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAg                                                                          <  1:1265423/150‑1 (MQ=255)
                                                                              ggTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGc                                                                         >  1:1069833/1‑150 (MQ=255)
                                                                                   tCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAc                                                                    >  2:198571/1‑150 (MQ=255)
                                                                                   tCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAc                                                                    >  2:14626/1‑150 (MQ=255)
                                                                                    ccGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAcc                                                                   >  2:1360465/1‑150 (MQ=255)
                                                                                    ccGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCCCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAcc                                                                   >  2:562316/1‑150 (MQ=255)
                                                                                      gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                                 <  1:1124212/150‑1 (MQ=255)
                                                                                           cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAt                                                            >  1:1072094/1‑150 (MQ=255)
                                                                                           cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATAGCTTTAAATATGACAGCATAACCTTTACAt                                                            >  1:763603/1‑150 (MQ=255)
                                                                                             gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATaa                                                          >  2:196967/1‑150 (MQ=255)
                                                                                             gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATaa                                                          >  2:1403932/1‑150 (MQ=255)
                                                                                                tCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                       >  2:413430/1‑150 (MQ=255)
                                                                                                tCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                       >  1:932378/1‑150 (MQ=255)
                                                                                                tCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                       >  2:801181/1‑150 (MQ=255)
                                                                                                 cATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTa                                                      >  1:666533/1‑150 (MQ=255)
                                                                                                  aTAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAg                                                     >  1:260790/1‑150 (MQ=255)
                                                                                                   tAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGt                                                    >  1:136484/1‑150 (MQ=255)
                                                                                                   tAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGt                                                    >  2:1612812/1‑150 (MQ=255)
                                                                                                    aaaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                   >  2:1533434/1‑150 (MQ=255)
                                                                                                      aCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTcc                                                 <  2:786448/150‑1 (MQ=255)
                                                                                                      aCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAGAATTTAGTTcc                                                 >  2:592794/1‑150 (MQ=255)
                                                                                                       cAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTAAAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCa                                                >  1:482612/1‑150 (MQ=255)
                                                                                                        aaaaTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAg                                               >  1:964695/1‑150 (MQ=255)
                                                                                                           aTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaa                                            >  1:59947/1‑150 (MQ=255)
                                                                                                              cAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACa                                         <  1:863936/150‑1 (MQ=255)
                                                                                                              cAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACa                                         >  2:551312/1‑150 (MQ=255)
                                                                                                               aaCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACaa                                        >  2:181224/1‑150 (MQ=255)
                                                                                                                aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                       >  1:1562892/1‑150 (MQ=255)
                                                                                                                aCAGCTAAACGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                       <  2:1468747/150‑1 (MQ=255)
                                                                                                                 cAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATc                                      <  1:611986/150‑1 (MQ=255)
                                                                                                                   gCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                    >  2:450722/1‑150 (MQ=255)
                                                                                                                        atgCTAAAACTATCATCCAGGTCATTATAGCCTGGCTTATTTTTATCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGa                               <  1:1031517/148‑1 (MQ=255)
                                                                                                                        aGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGa                               <  1:1472457/150‑1 (MQ=255)
                                                                                                                          gCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaa                             >  1:89519/1‑150 (MQ=255)
                                                                                                                           cTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaa                            >  1:849928/1‑150 (MQ=255)
                                                                                                                             ctaaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                          >  2:1460762/7‑150 (MQ=255)
                                                                                                                             aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGATAAAt                          >  1:34621/1‑150 (MQ=255)
                                                                                                                                aCTATCATCCAGGTCATTATCGACTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAt                       <  1:453428/150‑1 (MQ=255)
                                                                                                                                aCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAt                       <  2:1098477/150‑1 (MQ=255)
                                                                                                                                 aaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGAtt                      >  1:1550561/3‑150 (MQ=255)
                                                                                                                                  aaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATCATTc                     >  2:779111/2‑150 (MQ=255)
                                                                                                                                   aTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCa                    >  1:1394064/1‑150 (MQ=255)
                                                                                                                                   aTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCa                    >  1:936506/1‑150 (MQ=255)
                                                                                                                                   aTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCa                    >  1:297771/1‑150 (MQ=255)
                                                                                                                                    tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAg                   >  2:1000573/1‑150 (MQ=255)
                                                                                                                                     cAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt                  >  2:392053/1‑150 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 <  2:1507862/150‑1 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 >  2:917344/1‑150 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 <  2:903089/150‑1 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 <  1:1575680/150‑1 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 >  2:586528/1‑150 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 <  1:305583/150‑1 (MQ=255)
                                                                                                                                      aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 >  1:548469/1‑150 (MQ=255)
                                                                                                                                      aGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc                 >  1:58794/1‑150 (MQ=255)
                                                                                                                                       gCCAGGTCATTATCGCCTTGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa                <  2:1214989/150‑1 (MQ=255)
                                                                                                                                       gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa                >  2:4211/1‑150 (MQ=255)
                                                                                                                                       gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa                >  2:626882/1‑150 (MQ=255)
                                                                                                                                        ccAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa               >  2:646202/1‑150 (MQ=255)
                                                                                                                                        ccAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCaa               <  2:1072094/150‑1 (MQ=255)
                                                                                                                                         cAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAac              >  1:674361/1‑150 (MQ=255)
                                                                                                                                         cAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAac              <  1:1360465/150‑1 (MQ=255)
                                                                                                                                           ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac            >  2:25321/1‑150 (MQ=255)
                                                                                                                                              cATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTa         <  2:1544374/150‑1 (MQ=255)
                                                                                                                                               aTTATCGCCTGTCTGATTTTTATCTTAATGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAt        <  1:760700/150‑1 (MQ=255)
                                                                                                                                                    cGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCa   >  2:751587/1‑150 (MQ=255)
                                                                                                                                                     gCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCAt  >  1:122121/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTCCAT  >  NC_000913/1870112‑1870410

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