Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,870,257 T→A 5.5% intergenic (+19/‑128) 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,2570TA 5.5% 463.0 / 8.1 147intergenic (+19/‑128)yeaH/cdgIDUF444 domain‑containing protein YeaH/putative c‑di‑GMP binding protein CdgI
Reads supporting (aligned to +/- strand):  ref base T (70/68);  new base A (6/3);  total (76/71)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 148 total aligned reads displayed.

GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTT  >  NC_000913/1870111‑1870406
                                                                                                                                                  |                                                                                                                                                     
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)
       atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGc                                                                                                                                             >  2:87925/1‑150 (MQ=255)
       atatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGc                                                                                                                                             <  1:1112739/150‑1 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCt                                                                                                                                            <  1:1690708/150‑1 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCt                                                                                                                                            >  1:1843210/1‑150 (MQ=255)
        tatGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGCGATAATGACCTGGCt                                                                                                                                            <  2:1822124/150‑1 (MQ=255)
          tGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGa                                                                                                                                          >  1:1241711/1‑150 (MQ=255)
          tGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATATCGCCTGGCTGAt                                                                                                                                         >  1:1483520/1‑150 (MQ=255)
              cATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGAtttt                                                                                                                                      >  2:235842/1‑150 (MQ=255)
                  tGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGc                                                                                                                                  <  2:1838343/150‑1 (MQ=255)
                   gCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCt                                                                                                                                 >  2:135910/1‑150 (MQ=255)
                        cTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTg                                                                                                                            >  2:615665/1‑150 (MQ=255)
                        cTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATAATCGCCTGGCTGATTTTTAGCTTACTg                                                                                                                            >  1:1579832/1‑150 (MQ=255)
                          aCTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTTGCTGATTTTTAGCTTACTGTa                                                                                                                          >  1:1418658/1‑150 (MQ=255)
                           cTTTCGACAACTTTGCGATGCAGCACATCCGAGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAATCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCAATATCGACTGGCTGATTTTTAGCTTACTTTaa                                                                                                                         >  1:2276639/1‑150 (MQ=255)
                              tCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAATGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAAtt                                                                                                                      <  2:814003/150‑1 (MQ=255)
                               cGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTa                                                                                                                     >  1:1720470/1‑150 (MQ=255)
                               cGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTa                                                                                                                     >  1:1728920/1‑150 (MQ=255)
                                gACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAt                                                                                                                    >  2:2259537/1‑150 (MQ=255)
                                 aCAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTAtc                                                                                                                   <  2:1029586/150‑1 (MQ=255)
                                     cTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAt                                                                                                               >  2:528548/1‑150 (MQ=255)
                                      tttGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCAtt                                                                                                              <  2:1572955/150‑1 (MQ=255)
                                           gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                         >  2:1927884/1‑150 (MQ=255)
                                           gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                         <  1:930344/150‑1 (MQ=255)
                                           gATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTGTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATttatta                                                                                                         <  1:1694376/150‑1 (MQ=255)
                                               cagcaCATCCGCGACCAGGATGAAATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacata                                                                                                     >  1:1702831/1‑150 (MQ=255)
                                                 gcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                   <  1:1963944/150‑1 (MQ=255)
                                                 gcaCATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtacataca                                                                                                   <  2:2107707/150‑1 (MQ=255)
                                                   acaTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACAtt                                                                                                 <  2:1235680/150‑1 (MQ=255)
                                                    caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                >  1:1180904/1‑150 (MQ=255)
                                                    caTCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTa                                                                                                >  1:2270398/1‑150 (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)
                                                          cgACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTAc                                                                                          <  2:651614/150‑1 (MQ=255)
                                                              cAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATc                                                                                      >  1:2046046/1‑150 (MQ=255)
                                                                ggATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGc                                                                                    <  2:1847204/150‑1 (MQ=255)
                                                                 gatgaaATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACTTTAGCTTACAATCGCt                                                                                   <  2:1702831/150‑1 (MQ=255)
                                                                  atgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGcgt                                                                                  >  2:854058/1‑148 (MQ=255)
                                                                   tgatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCttt                                                                                 >  2:169243/1‑150 (MQ=255)
                                                                    gatATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTCGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATTCATTAGCTTACAATCGCTTTa                                                                                >  1:1223458/1‑150 (MQ=255)
                                                                    gatATTTATCCGGTGGTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTa                                                                                <  2:593113/150‑1 (MQ=255)
                                                                      tatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTaaa                                                                              <  2:1649480/150‑1 (MQ=255)
                                                                      tatTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTaaa                                                                              <  2:645838/150‑1 (MQ=255)
                                                                       atTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAat                                                                             >  1:378116/1‑150 (MQ=255)
                                                                        tttATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAata                                                                            <  1:576284/150‑1 (MQ=255)
                                                                         ttATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAatat                                                                           >  1:604519/1‑150 (MQ=255)
                                                                          tATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATg                                                                          <  2:241815/150‑1 (MQ=255)
                                                                           aTCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGa                                                                         >  2:2363113/1‑150 (MQ=255)
                                                                            tCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGAc                                                                        <  1:235842/150‑1 (MQ=255)
                                                                             ccGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACa                                                                       >  1:1329901/1‑150 (MQ=255)
                                                                             ccGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACa                                                                       >  1:148035/1‑150 (MQ=255)
                                                                               ggTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGc                                                                     <  2:2075177/150‑1 (MQ=255)
                                                                                 tgtTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCAt                                                                   <  2:588384/150‑1 (MQ=255)
                                                                                      cGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCt                                                              <  1:1261258/150‑1 (MQ=255)
                                                                                       gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACTTTAGCTTACAATCGATTTAAATATGAAAGCATAACAtt                                                             >  2:2136054/1‑150 (MQ=255)
                                                                                       gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                             <  1:1466124/150‑1 (MQ=255)
                                                                                       gTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCtt                                                             >  2:533436/1‑150 (MQ=255)
                                                                                         taaCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTa                                                           <  2:449272/149‑1 (MQ=255)
                                                                                          aaCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTAc                                                          <  2:1241711/150‑1 (MQ=255)
                                                                                           aCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACa                                                         >  1:1562019/1‑150 (MQ=255)
                                                                                           aCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACa                                                         <  1:1441389/150‑1 (MQ=255)
                                                                                           aCTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACa                                                         <  1:637869/150‑1 (MQ=255)
                                                                                            cTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACAt                                                        <  2:2474274/150‑1 (MQ=255)
                                                                                                 tCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                   <  2:1180904/150‑1 (MQ=255)
                                                                                                 tCATAAACAAAATGCAACAGATAAAGGCTAAAACTATCAGCCAGGTCATTATCTCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAAttt                                                   <  1:2163985/150‑1 (MQ=255)
                                                                                                  cATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTa                                                  >  1:1663654/1‑150 (MQ=255)
                                                                                                      aaCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTc                                              >  1:1979509/1‑150 (MQ=255)
                                                                                                       aCAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTcc                                             >  1:2077285/1‑150 (MQ=255)
                                                                                                           aaTGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaa                                         >  2:681096/1‑150 (MQ=255)
                                                                                                             tGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATtac                                                                                                        <  1:1042437/85‑1 (MQ=255)
                                                                                                             tGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGaaaa                                       >  2:1767410/1‑150 (MQ=255)
                                                                                                                    gCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                >  2:2049372/1‑150 (MQ=255)
                                                                                                                    gCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCAt                                >  1:287560/1‑150 (MQ=255)
                                                                                                                       aaaGGCTAAAACTATCAGCCAGGTCATATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCgg                            <  2:1483520/150‑1 (MQ=255)
                                                                                                                             tAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTATTTCCAGAAAACAATCATTCGGaaaaa                       <  2:2303430/150‑1 (MQ=255)
                                                                                                                             tAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGaaaaa                       >  2:50190/1‑150 (MQ=255)
                                                                                                                              ctaaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGAATAAACTTTACATAATTTAGATCCAGAAAACAATCATTCGGAAAAAt                      >  2:2227438/7‑150 (MQ=255)
                                                                                                                              aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                      <  2:1492221/150‑1 (MQ=255)
                                                                                                                              aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                      >  2:1211799/1‑150 (MQ=255)
                                                                                                                              aaaaCTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAAt                      >  2:1207659/1‑150 (MQ=255)
                                                                                                                                aaaaaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGa                    >  1:109562/5‑150 (MQ=255)
                                                                                                                                   aaTCAGCCAGGCGATAATGACCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTc                 >  1:1458673/2‑150 (MQ=255)
                                                                                                                                     tCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAg               <  1:1912074/150‑1 (MQ=255)
                                                                                                                                      cAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGt              >  1:2119815/1‑150 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc             >  1:2457819/1‑150 (MQ=255)
                                                                                                                                       aGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTc             <  2:1979509/150‑1 (MQ=255)
                                                                                                                                        gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa            >  2:2049283/1‑150 (MQ=255)
                                                                                                                                        gCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTCTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCa            <  2:1584851/150‑1 (MQ=255)
                                                                                                                                           aGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAaca         <  2:818349/150‑1 (MQ=255)
                                                                                                                                            ggTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAacac        >  1:1176075/1‑150 (MQ=255)
                                                                                                                                             gTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACg       <  1:1338771/150‑1 (MQ=255)
                                                                                                                                               cATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTa     <  1:533436/150‑1 (MQ=255)
                                                                                                                                                aTTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAAGCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAt    >  1:189648/1‑150 (MQ=255)
                                                                                                                                                  tATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTAttt  >  1:538304/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GGCGAGAATATGAGCATCTGCAATCTACTTTCGACAACTTTGCGATGCAGCACATCCGCGACCAGGATGATATTTATCCGGTGTTCCGTGAACTGTTTCATAAACAAAATGCAACAGCTAAAGGCTAAAACTATCAGCCAGGTCATTATCGCCTGGCTGATTTTTAGCTTACTGTAAATTATCTCATTTATTACATACATTAGCTTACAATCGCTTTAAATATGACAGCATAACCTTTACATAATTTAGTTCCAGAAAACAATCATTCGGAAAAATGATTCAGTCAACACGTATTT  >  NC_000913/1870111‑1870406

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: