Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,897,910 G→A 6.9% intergenic (‑40/+54) ygcS ← / ← ygcU putative MFS sugar transporter; membrane protein/putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,897,9100GA 6.9% 353.5 / 8.8 115intergenic (‑40/+54)ygcS/ygcUputative MFS sugar transporter; membrane protein/putative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (53/54);  new base A (2/6);  total (55/60)
Fisher's exact test for biased strand distribution p-value = 2.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 116 total aligned reads displayed.

CCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGTG  >  NC_000913/2897764‑2898059
                                                                                                                                                  |                                                                                                                                                     
ccGAGAACATAACCGTCGATCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCt                                                                                                                                                    >  2:195604/1‑150 (MQ=255)
 cGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTg                                                                                                                                                   >  2:1147959/1‑150 (MQ=255)
  gagaACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTga                                                                                                                                                  <  2:921499/150‑1 (MQ=255)
     aaCATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGaa                                                                                                                                               >  2:518521/1‑150 (MQ=255)
      aCATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAg                                                                                                                                              >  1:1396234/1‑150 (MQ=255)
       cATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGc                                                                                                                                             >  2:1034796/1‑150 (MQ=255)
        aTAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGcc                                                                                                                                            <  1:1086471/150‑1 (MQ=255)
        aTAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGcc                                                                                                                                            <  1:360523/150‑1 (MQ=255)
         tAACCGTCGGTCCGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCTGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCg                                                                                                                                           <  1:1006337/150‑1 (MQ=255)
         tAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCg                                                                                                                                           <  1:1565608/150‑1 (MQ=255)
          aaCCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGa                                                                                                                                          >  2:789380/1‑150 (MQ=255)
          aaCCGTCGGGCATGTGTGCGCCGAAAGTGAGCGCAGCAACGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgt                                                                                                                                          <  1:304762/150‑2 (MQ=255)
            ncNTCGGTCAGGTGTGCGCCGAAAGTNAGCGCAGCANTNCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGATGCCGaaa                                                                                                                                        >  2:1011281/2‑150 (MQ=255)
            ccGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttc                                                                                                                                        <  1:836854/150‑4 (MQ=255)
             cGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttct                                                                                                                                       <  1:818699/150‑5 (MQ=255)
              gTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAg                                                                                                                                      <  1:528857/150‑1 (MQ=255)
              gTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCAGGCGAGAAGCCgttctt                                                                                                                                      <  2:425465/150‑6 (MQ=255)
              gTCGGGCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAg                                                                                                                                      <  2:41415/150‑1 (MQ=255)
                  gTCATGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAACCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAATTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCaa                                                                                                                                  <  1:1045038/150‑1 (MQ=255)
                  gTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGAGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACACCTTCATTGTCCGGTATTAAGAACGGCTTCTCGACTGAGAAGCCGAAAAGGCaa                                                                                                                                  >  2:1271984/1‑150 (MQ=255)
                   gcAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAAc                                                                                                                                 <  1:124957/149‑1 (MQ=255)
                      ggTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCTTTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGcc                                                                                                                              <  2:484344/150‑1 (MQ=255)
                      ggTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGcc                                                                                                                              >  1:76044/1‑150 (MQ=255)
                      ggTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGACTTCATCTTTACCCCTTCATTGTACGGTAATAAGAACGGCTTATCGCCAGAGAAGCCGAAAACGCAACGcc                                                                                                                              >  2:1622563/1‑150 (MQ=255)
                             gcagAAAGTGAGCTCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCAt                                                                                                                       <  1:537375/147‑1 (MQ=255)
                             gcCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCAt                                                                                                                       >  2:101332/1‑150 (MQ=255)
                              ccGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATc                                                                                                                      >  1:1189022/1‑150 (MQ=255)
                                  aaGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAAGCACCGCCCCGTCATCGCgg                                                                                                                  >  2:1022750/1‑150 (MQ=255)
                                   aGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGTAAAGGCAACGCCCCGTCATCGCGGa                                                                                                                 >  2:1176548/1‑150 (MQ=255)
                                          gcagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGc                                                                                                          <  2:246540/150‑1 (MQ=255)
                                           cagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGcc                                                                                                         <  1:789380/150‑1 (MQ=255)
                                           cagcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGcc                                                                                                         >  1:1052348/1‑150 (MQ=255)
                                            agcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCt                                                                                                        <  1:232814/150‑1 (MQ=255)
                                             gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCAGCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                       <  2:1040936/150‑1 (MQ=255)
                                             gcaATGCGGCAGTGAAAACGGTTAAGCGGTAAAGCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTGCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                       <  2:207123/150‑1 (MQ=255)
                                             gcaATGCGGCAGTGAAAACGGCTAAGCGGTAAATCCTCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTg                                                                                                       <  2:1719897/150‑1 (MQ=255)
                                               aaTGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGac                                                                                                     <  2:49629/150‑1 (MQ=255)
                                                 tGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGacac                                                                                                   >  2:1385064/1‑150 (MQ=255)
                                                  gCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGAcaca                                                                                                  >  2:391284/1‑150 (MQ=255)
                                                     gCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACatta                                                                                               <  1:1209651/150‑1 (MQ=255)
                                                      cAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACattat                                                                                              >  1:1167454/1‑150 (MQ=255)
                                                        gTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGATCGCTGACTGACACATTAttt                                                                                            >  2:215586/1‑150 (MQ=255)
                                                           aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                         <  2:82576/150‑1 (MQ=255)
                                                           aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTc                                                                                         >  2:909636/1‑150 (MQ=255)
                                                           aaaaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAAAAAGAACGGCTTCTCGCCTGACAAGCCGAAAAGGCAACGCCCCGTCATAGCGGAACGCTGCCTGACACATTATTTTTc                                                                                         >  2:192982/1‑150 (MQ=255)
                                                           aaaaCGGTAAAGCGGTAAATCATACATTCGAACCGGTGAAGTGTTAATATTTACCCCTTCATTGTCCGGTAATAAAAACGCCTTCTCGCATGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCATGACACattattatta                                                                                         >  2:495155/1‑149 (MQ=255)
                                                             aaCGGTTAAGCGGTAAATCATCCATTCGCACCGGTGTAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGTACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCaa                                                                                       <  2:1641427/150‑1 (MQ=255)
                                                             aaCGGTTAAGCGGGAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCaa                                                                                       <  1:851155/150‑1 (MQ=255)
                                                                 gTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCgg                                                                                   <  1:761350/150‑1 (MQ=255)
                                                                      gCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgataga                                                                              >  1:391171/1‑150 (MQ=255)
                                                                        ggTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGTCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGgatagata                                                                            >  2:366874/1‑150 (MQ=255)
                                                                           aaaTCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTa                                                                         <  1:436147/150‑1 (MQ=255)
                                                                             atcatcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGACGCCGAAAAGGCACCGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAcc                                                                       >  2:873681/1‑150 (MQ=255)
                                                                             atcatcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAcc                                                                       <  1:1147959/150‑1 (MQ=255)
                                                                             atcatcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAcc                                                                       <  2:1444009/150‑1 (MQ=255)
                                                                             atcatcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTAcc                                                                       <  2:390002/150‑1 (MQ=255)
                                                                               catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCTGTTATAAGAACGGCTTCTCGCCTGAGAAGCAGAAAAGGGAACGCCCCGTCATCGCGGAACGCTGCCTGAGACATTATTTTTCTATCGGATAGATAGTAcgcg                                                                     <  2:217896/150‑1 (MQ=255)
                                                                               catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                     <  1:1415252/150‑1 (MQ=255)
                                                                               catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                     <  2:448838/150‑1 (MQ=255)
                                                                               catcCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCg                                                                     >  2:1690604/1‑150 (MQ=255)
                                                                                   cATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGCTAGTACCCGTAtt                                                                 >  2:1262928/1‑150 (MQ=255)
                                                                                   cATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTAtt                                                                 >  2:221005/1‑150 (MQ=255)
                                                                                     ttCGCACCGGTGAAGAGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCa                                                               <  1:1483741/150‑1 (MQ=255)
                                                                                      tCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCAt                                                              <  1:114687/150‑1 (MQ=255)
                                                                                       cGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATa                                                             <  1:1264566/150‑1 (MQ=255)
                                                                                           ccGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATTAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCAGTGATCGCGGATCGCTGGCTGACACATTTTTTTTCATTCGGATAGTTAGTTCCCGTATTctttttgg                                                         >  2:1497039/1‑143 (MQ=255)
                                                                                               tGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCAtt                                                     <  2:76044/150‑1 (MQ=255)
                                                                                               ggAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCAtt                                                     <  2:729301/149‑1 (MQ=255)
                                                                                                     gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                               >  1:1253439/1‑150 (MQ=255)
                                                                                                     gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                               >  1:700588/1‑150 (MQ=255)
                                                                                                     gtTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATc                                               >  1:694514/1‑150 (MQ=255)
                                                                                                      ttCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                              <  1:656898/150‑1 (MQ=255)
                                                                                                      ttCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCg                                              <  1:1710378/150‑1 (MQ=255)
                                                                                                       tCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATATAGCCATTAGGATCGa                                             >  2:797235/1‑150 (MQ=255)
                                                                                                        cATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGaa                                            >  1:1393635/1‑150 (MQ=255)
                                                                                                         aTCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAAc                                           <  1:704210/150‑1 (MQ=255)
                                                                                                            tttACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTAGCTTGATAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCGTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGc                                        <  1:88642/150‑1 (MQ=255)
                                                                                                                  ccTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACACGTATTCATAATGCCATTAGGATCGAACTGCTTTTTc                                  >  1:304567/1‑150 (MQ=255)
                                                                                                                  ccTTCATTGACCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACTCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATACATAATGCCATTAGGATCGAACTGCTTTTTc                                  >  1:846809/1‑150 (MQ=255)
                                                                                                                     tCATTGTCCTGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGc                               <  2:583204/150‑1 (MQ=255)
                                                                                                                     tCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGc                               <  2:887782/150‑1 (MQ=255)
                                                                                                                      cATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGAGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGcc                              <  1:126957/150‑1 (MQ=255)
                                                                                                                           tCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTcc                         <  2:1676772/150‑1 (MQ=255)
                                                                                                                           tCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTcc                         <  1:391284/150‑1 (MQ=255)
                                                                                                                             cGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAg                       <  1:195604/150‑1 (MQ=255)
                                                                                                                             cGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAg                       >  1:17878/1‑150 (MQ=255)
                                                                                                                              ggTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGc                      >  1:965481/1‑150 (MQ=255)
                                                                                                                                 aataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAAc                   >  2:1043611/1‑150 (MQ=255)
                                                                                                                                  ataaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACAATATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                  >  2:1189581/1‑150 (MQ=255)
                                                                                                                                  aaaaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTAAAAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACg                  <  1:1497039/148‑1 (MQ=255)
                                                                                                                                   ccttttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGc                 >  1:1218562/7‑150 (MQ=255)
                                                                                                                                    aaGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGcc                <  1:1532547/150‑1 (MQ=255)
                                                                                                                                      tttcGGCTTCTCAGGCGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGAGAGATAGTACCCGTATTCATTATGCCATTAGGATCGAACTGCTTTTTCAGACCTTCCCGCAACGcccc              >  1:1549116/4‑149 (MQ=255)
                                                                                                                                      gAACGGCTTCTCGCCTGAGAAGCCGAAAAGTATACGCCCCGTCATCGCGGAACGCTTCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCa              <  1:192982/150‑1 (MQ=255)
                                                                                                                                       aaCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAc             <  2:420411/150‑1 (MQ=255)
                                                                                                                                       aaCGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAc             >  1:1110053/1‑150 (MQ=255)
                                                                                                                                          ggCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCAcgcg          >  2:1465869/1‑150 (MQ=255)
                                                                                                                                           gCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACgcgc         <  1:1622563/150‑1 (MQ=255)
                                                                                                                                           gCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGGAGTCATATTGCCATTAGGATCGCACTGCTTTTTCAGCCCTTCCAGCAACGCCCACtcgc         >  1:1422887/1‑150 (MQ=255)
                                                                                                                                                  gCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCgtg  <  2:1263972/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CCGAGAACATAACCGTCGGTCAGGTGTGCGCCGAAAGTGAGCGCAGCAATGCGGCAGTGAAAACGGTTAAGCGGTAAATCATCCATTCGCACCGGTGAAGTGTTCATCTTTACCCCTTCATTGTCCGGTAATAAGAACGGCTTCTCGCCTGAGAAGCCGAAAAGGCAACGCCCCGTCATCGCGGAACGCTGCCTGACACATTATTTTTCAATCGGATAGATAGTACCCGTATTCATAATGCCATTAGGATCGAACTGCTTTTTCAGCCCTTCCAGCAACGCCCACGCGCTGCCGTG  >  NC_000913/2897764‑2898059

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