Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,260 C→G 5.5% intergenic (+22/‑125) 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,2600CG 5.5% 503.5 / 11.9 150intergenic (+22/‑125)yeaH/cdgIDUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base C (73/68);  new base G (5/3);  total (78/72)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.06e-01

Only 100 of 151 total aligned reads displayed.

GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870111‑1870407
                                                                                                                                                     |                                                                                                                                                   
ggCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATc                                                                                                                                                     >  2:188194/1‑150 (MQ=255)
  cGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGc                                                                                                                                                   <  2:2314324/150‑1 (MQ=255)
  cGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGc                                                                                                                                                   >  1:2474274/1‑150 (MQ=255)
   gagaATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGcc                                                                                                                                                  <  2:1420787/150‑1 (MQ=255)
      aaTATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCATgg                                                                                                                                               >  1:1207088/1‑150 (MQ=255)
       atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTggg                                                                                                                                              >  1:1584851/1‑149 (MQ=255)
       atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGc                                                                                                                                              >  2:87925/1‑150 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCt                                                                                                                                             >  1:1843210/1‑150 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCt                                                                                                                                             <  2:1822124/150‑1 (MQ=255)
         atGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTg                                                                                                                                            >  1:1474146/1‑150 (MQ=255)
         atGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCTg                                                                                                                                            <  1:2109516/150‑1 (MQ=255)
          tGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGa                                                                                                                                           >  1:1241711/1‑150 (MQ=255)
             gCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCTGAtag                                                                                                                                        <  2:1213995/150‑3 (MQ=255)
              cATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGAtttt                                                                                                                                       <  2:554774/150‑1 (MQ=255)
                  tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                   <  1:1225407/150‑1 (MQ=255)
                  tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                   <  2:1838343/150‑1 (MQ=255)
                     aaTCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTa                                                                                                                                <  1:1313719/150‑1 (MQ=255)
                      aTCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTAc                                                                                                                               <  2:162463/150‑1 (MQ=255)
                          aCTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTTGCTGATTTTTAGCTTACTGTa                                                                                                                           >  1:1418658/1‑150 (MQ=255)
                          aCTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTa                                                                                                                           <  1:251952/150‑1 (MQ=255)
                           cTTTCGACAACTTTGCGATGCAGCACATCCGAGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAATCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCAATATCGACTGGCTGATTTTTAGCTTACTTTaa                                                                                                                          >  1:2276639/1‑150 (MQ=255)
                               cGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTa                                                                                                                      >  1:1720470/1‑150 (MQ=255)
                                gACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAt                                                                                                                     >  2:2259537/1‑150 (MQ=255)
                                 aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtc                                                                                                                    <  1:1778199/150‑1 (MQ=255)
                                 aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtc                                                                                                                    <  2:1029586/150‑1 (MQ=255)
                                  cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtct                                                                                                                   >  2:1448234/1‑150 (MQ=255)
                                  cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtct                                                                                                                   >  2:1338771/1‑150 (MQ=255)
                                  cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtct                                                                                                                   <  2:2262645/150‑1 (MQ=255)
                                  cAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtct                                                                                                                   >  2:394014/1‑150 (MQ=255)
                                   aaCTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtctc                                                                                                                  <  2:2417646/150‑1 (MQ=255)
                                     cTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTCTTAGCTTACTGTAAATTATCTCAt                                                                                                                >  1:458549/1‑150 (MQ=255)
                                           gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                          >  2:1927884/1‑150 (MQ=255)
                                           gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTGTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                          <  1:1694376/150‑1 (MQ=255)
                                                agcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacatac                                                                                                     <  2:650572/150‑1 (MQ=255)
                                                 gcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                    <  2:2107707/150‑1 (MQ=255)
                                                   acaTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACAtt                                                                                                  <  2:1235680/150‑1 (MQ=255)
                                                    caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                 <  1:2376547/150‑1 (MQ=255)
                                                     aTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCATCGTTGTCATTTTTGCGTGGATGATTTTTAGCTTACTGTTAATTATCTCATTTATTACATTAATTAg                                                                                                <  2:693966/150‑1 (MQ=255)
                                                        cgcgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCtt                                                                                             >  1:674860/1‑150 (MQ=255)
                                                          cgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTAc                                                                                           <  2:660014/150‑1 (MQ=255)
                                                               aGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCg                                                                                      >  1:2242238/1‑150 (MQ=255)
                                                               aGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCg                                                                                      >  1:1492221/1‑150 (MQ=255)
                                                                ggATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGc                                                                                     <  2:1847204/150‑1 (MQ=255)
                                                                 gatgaaATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACTTTAGCTTACAATCGCt                                                                                    <  2:1702831/150‑1 (MQ=255)
                                                                  atgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCtt                                                                                   <  1:1927884/150‑1 (MQ=255)
                                                                  atgatATTTATCCGGTGTTCCGTGAACCGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCtt                                                                                   <  2:763660/150‑1 (MQ=255)
                                                                   tgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCttt                                                                                  >  2:169243/1‑150 (MQ=255)
                                                                    gatATTTATCCGGTGGTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTa                                                                                 <  2:593113/150‑1 (MQ=255)
                                                                      tatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTaaa                                                                               <  2:645838/150‑1 (MQ=255)
                                                                       atTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAat                                                                              >  1:378116/1‑150 (MQ=255)
                                                                         ttATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAatat                                                                            >  1:604519/1‑150 (MQ=255)
                                                                          tATCCGGTGTTCCGTTAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                           <  1:556873/150‑1 (MQ=255)
                                                                           aTCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGa                                                                          >  2:2363113/1‑150 (MQ=255)
                                                                             ccGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACa                                                                        >  1:1329901/1‑150 (MQ=255)
                                                                             ccGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACa                                                                        >  1:148035/1‑150 (MQ=255)
                                                                               ggTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGc                                                                      <  2:2075177/150‑1 (MQ=255)
                                                                                  gtTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATa                                                                   >  1:1452067/1‑150 (MQ=255)
                                                                                   ttCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATaa                                                                  <  1:800056/150‑1 (MQ=255)
                                                                                       gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACTTTAGCTTACAATCGATTTAAATATGAAAGCATAACAtt                                                              >  2:2136054/1‑150 (MQ=255)
                                                                                       gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                              >  2:533436/1‑150 (MQ=255)
                                                                                           aCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACa                                                          >  1:1562019/1‑150 (MQ=255)
                                                                                            cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAt                                                         <  2:2474274/150‑1 (MQ=255)
                                                                                              gTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATaa                                                       >  1:1986017/1‑150 (MQ=255)
                                                                                                 tCATAAACAAAATGCAACAGATAAAGGCTAAAACTATCAGCCAGGTCATTATCTCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                    <  1:2163985/150‑1 (MQ=255)
                                                                                                  cATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTa                                                   >  1:1663654/1‑150 (MQ=255)
                                                                                                     aaaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGtt                                                >  2:849260/1‑150 (MQ=255)
                                                                                                      aaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  1:1979509/1‑150 (MQ=255)
                                                                                                      aaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                               >  2:60662/1‑150 (MQ=255)
                                                                                                       aCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTcc                                              >  1:2077285/1‑150 (MQ=255)
                                                                                                             tGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtac                                                                                                         <  1:1042437/85‑1 (MQ=255)
                                                                                                             tGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtac                                                                                                         >  2:1042437/1‑85 (MQ=255)
                                                                                                             tGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaaa                                        >  2:1767410/1‑150 (MQ=255)
                                                                                                                 aCAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAAt                                    <  1:135910/150‑1 (MQ=255)
                                                                                                                  cAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATc                                   >  1:119042/1‑150 (MQ=255)
                                                                                                                    gCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                 >  2:2049372/1‑150 (MQ=255)
                                                                                                                    gCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                 >  1:2330437/1‑150 (MQ=255)
                                                                                                                     cTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAtt                                >  2:2348993/1‑150 (MQ=255)
                                                                                                                       aaaGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCg                              <  1:1143516/150‑1 (MQ=255)
                                                                                                                            cTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaa                         >  1:1718607/1‑150 (MQ=255)
                                                                                                                             tAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                        >  2:50190/1‑150 (MQ=255)
                                                                                                                             gctaaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTTCTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                        >  2:525259/8‑150 (MQ=38)
                                                                                                                              aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                       >  2:1207659/1‑150 (MQ=255)
                                                                                                                              aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                       <  2:1492221/150‑1 (MQ=255)
                                                                                                                               aaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATg                      <  1:615665/150‑1 (MQ=255)
                                                                                                                                aaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGa                     >  1:109562/5‑150 (MQ=255)
                                                                                                                                  cTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACATTTACATAATTTAGTTACAGAAAACAATCATTCGGAAAAATGAtt                   >  2:1538316/1‑150 (MQ=255)
                                                                                                                                   aaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTc                  >  1:1458673/2‑150 (MQ=255)
                                                                                                                                    aTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCa                 >  2:667995/1‑150 (MQ=255)
                                                                                                                                     tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAg                <  1:1912074/150‑1 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc              <  2:1979509/150‑1 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc              <  1:169243/150‑1 (MQ=255)
                                                                                                                                        gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa             >  2:2049283/1‑150 (MQ=255)
                                                                                                                                        gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTCTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa             <  2:1584851/150‑1 (MQ=255)
                                                                                                                                          cAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAac           >  1:2385259/1‑150 (MQ=255)
                                                                                                                                           aGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAaca          <  2:818349/150‑1 (MQ=255)
                                                                                                                                            ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac         >  1:1176075/1‑150 (MQ=255)
                                                                                                                                                  tATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAttt   >  1:538304/1‑150 (MQ=255)
                                                                                                                                                   aTCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc  >  2:1305680/1‑150 (MQ=255)
                                                                                                                                                   aTCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc  <  2:1678564/150‑1 (MQ=255)
                                                                                                                                                   aTCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTc  >  2:2207791/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTTC  >  NC_000913/1870111‑1870407

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: