Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,774 G→A 6.3% intergenic (‑14/+54) ycjD ← / ← sapF DUF559 family endonuclease‑related protein/antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,351,7740GA 6.3% 494.5 / 14.1 158intergenic (‑14/+54)ycjD/sapFDUF559 family endonuclease‑related protein/antimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (74/73);  new base A (4/6);  total (78/80)
Fisher's exact test for biased strand distribution p-value = 7.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 160 total aligned reads displayed.

AGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCA  >  NC_000913/1351625‑1351922
                                                                                                                                                     |                                                                                                                                                    
agatCCCCTCTGGATGTTGCCGGCGTAATTTGAAATCACTAAAACGAGTGCTTCTAAGATATCGCCAGAGATTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGAATTAATTTAATCTATCACCTCATTCTGACAAg                                                                                                                                                      <  1:321542/147‑1 (MQ=255)
 gCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGa                                                                                                                                                     >  2:436148/1‑150 (MQ=255)
 gCTCCCCACTGGATGTTGACGGCGAAATTTGAAATAACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTAATGCAAAGTGAGATTGCGTCTTAAATAAAGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGa                                                                                                                                                     >  1:1026307/1‑150 (MQ=255)
   tctccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGAtt                                                                                                                                                   <  2:1915214/147‑1 (MQ=255)
    ccccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGAttt                                                                                                                                                  >  2:737491/1‑150 (MQ=255)
      ccACTGGATGTTGACGCCGAAATTTGAAAACACTAAAACGTCGGCAGCGAAGATATCGCCAGAGCTTTCGTACCTGCAAAGTGAGATTGCGTCTTAAATAACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTaa                                                                                                                                                >  2:143006/1‑150 (MQ=255)
       cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                               >  1:861076/1‑150 (MQ=255)
           ggATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCttt                                                                                                                                           >  1:1391495/1‑150 (MQ=255)
             aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCTGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGAATGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                         <  1:346820/150‑1 (MQ=255)
             aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                         <  2:1941621/150‑1 (MQ=255)
              tGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                                                                                                                        <  1:682061/150‑1 (MQ=255)
              tGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                                                                                                                        >  2:1734400/1‑150 (MQ=255)
                ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTTCAAAGTGCGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                      <  2:2012741/150‑1 (MQ=255)
                ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                      <  2:1750084/150‑1 (MQ=255)
                ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                      <  1:1127244/150‑1 (MQ=255)
                ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                      <  2:1551047/150‑1 (MQ=255)
                 tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAg                                                                                                                                     <  1:1625610/150‑1 (MQ=255)
                  gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGa                                                                                                                                    >  1:152651/1‑150 (MQ=255)
                    cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAAt                                                                                                                                  >  2:162234/1‑150 (MQ=255)
                     ggCGAAATTTGAACTCACTAAAACGTCGGCTGGGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATg                                                                                                                                 <  1:309238/150‑1 (MQ=38)
                     ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATg                                                                                                                                 <  1:632232/150‑1 (MQ=38)
                     ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATg                                                                                                                                 <  1:612419/150‑1 (MQ=38)
                     ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATACCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATg                                                                                                                                 >  2:1889746/1‑150 (MQ=255)
                      gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGa                                                                                                                                <  2:601074/150‑1 (MQ=38)
                       cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGAGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAg                                                                                                                               <  2:1832220/150‑1 (MQ=38)
                        gAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCTCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAgg                                                                                                                              <  1:909469/150‑1 (MQ=38)
                        gAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAgg                                                                                                                              <  1:1263872/150‑1 (MQ=255)
                         aaaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGt                                                                                                                             <  1:2023128/150‑1 (MQ=255)
                            tttGAAATCACTAAAACGTCGGCTGCGAAGATATAGCCCGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATATAATTTTATCTATCACCTCATTCTGACAAGATTTAATCGTTAGGCAGAATGATGTGaa                                                                                                                          >  2:1019389/1‑150 (MQ=255)
                                     aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                 >  1:23560/1‑150 (MQ=255)
                                        aaaaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAagg                                                                                                              >  1:780772/1‑150 (MQ=255)
                                            cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                          <  2:842091/150‑1 (MQ=255)
                                            cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                          >  2:1708590/1‑150 (MQ=255)
                                               cGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGCCAAGATTTAATCTTTTGTCAGAAGGAGGTGAATAATCTGGAAGGAGGATTc                                                                                                       >  2:1272552/1‑150 (MQ=255)
                                                  cTGAGAAGATATCGCCAGATCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGa                                                                                                    <  2:1071107/150‑1 (MQ=255)
                                                   tGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTAGTCAGAATGAGGTGAATAAACTGGAAGTAGGATTCAGaa                                                                                                   >  2:1434735/1‑150 (MQ=255)
                                                     cGAAGATTTCGCCAGAGGTTTCGTTCCTGCATTGTGAGTTTGCGTCTTAAATCACGCGCATTTGATTTAAATTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                 <  1:187858/150‑1 (MQ=255)
                                                     cGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                 >  1:226542/1‑150 (MQ=255)
                                                         gATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAg                                                                                             <  2:1051477/150‑1 (MQ=255)
                                                          atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                            >  1:596455/1‑150 (MQ=255)
                                                            atCGCGAGAGCTTTCGTTCCTGCAATGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTTGAAGGAGGATTCAGAAAATTAGCGa                                                                                          >  1:140988/1‑150 (MQ=255)
                                                              cGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATc                                                                                        >  2:372054/1‑150 (MQ=255)
                                                               gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                       >  1:587528/1‑150 (MQ=255)
                                                               gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTCAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                       >  1:154699/1‑150 (MQ=255)
                                                                 cAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAACTCACGCGCATTTGAGTTAATTTTATCTATCAACTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTCGAATGAGGATTCAGAAAAATAGCGATCttt                                                                                     >  1:178802/1‑150 (MQ=255)
                                                                 cAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCttt                                                                                     <  1:767887/150‑1 (MQ=255)
                                                                 cAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCttt                                                                                     <  1:1480449/150‑1 (MQ=255)
                                                                         tCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACg                                                                             <  2:1314705/150‑1 (MQ=255)
                                                                          cGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                                                                            >  1:419010/1‑150 (MQ=255)
                                                                          cGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                                                                            >  2:414764/1‑150 (MQ=255)
                                                                              ccTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATcc                                                                        >  1:451036/1‑150 (MQ=255)
                                                                               cTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACTCCACGCATccg                                                                       >  1:1596711/1‑150 (MQ=255)
                                                                                  cAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATccgccg                                                                    >  2:1888364/1‑150 (MQ=255)
                                                                                   aaaGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGt                                                                   >  1:1820304/1‑150 (MQ=255)
                                                                                     aGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCa                                                                 >  2:88440/1‑150 (MQ=255)
                                                                                        gagaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATg                                                              >  2:1602567/1‑150 (MQ=255)
                                                                                         agaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCTCCGTCAATGc                                                             >  1:979693/1‑150 (MQ=255)
                                                                                          gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                            <  1:902794/150‑1 (MQ=255)
                                                                                          gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                            <  2:219800/150‑1 (MQ=255)
                                                                                             tGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCg                                                         >  2:1479682/1‑150 (MQ=255)
                                                                                               cGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                       <  2:399371/150‑1 (MQ=255)
                                                                                               cGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                       <  1:1244283/150‑1 (MQ=255)
                                                                                                gTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCa                                                      <  1:556926/150‑1 (MQ=255)
                                                                                                 tCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaa                                                     >  2:1530837/1‑150 (MQ=255)
                                                                                                  cTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGTCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCaaa                                                    >  2:914309/1‑150 (MQ=255)
                                                                                                    tAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAAt                                                  >  2:1181477/1‑150 (MQ=255)
                                                                                                     aaaTCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATg                                                 <  1:436148/150‑1 (MQ=255)
                                                                                                       aTCACGCTCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCATAATGATGCGATTAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACTCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAc                                               <  2:1026306/150‑1 (MQ=255)
                                                                                                          actcgcATTTGATTTAATATTATCTATCACCTCATTCTGACAAGATTTAATCTTTTTTCAGAATGAGGTGAATAATCTTGACGGAGGATTCAGAAAATTAGCGATCTTTCCGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCgg                                            <  2:161884/147‑1 (MQ=255)
                                                                                                          aCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCgg                                            <  2:494214/150‑1 (MQ=255)
                                                                                                           cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  1:1485080/150‑1 (MQ=255)
                                                                                                           cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  1:759748/150‑1 (MQ=255)
                                                                                                           cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  2:1647837/150‑1 (MQ=255)
                                                                                                           cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           <  2:1639765/150‑1 (MQ=255)
                                                                                                           cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                           >  2:1630004/1‑150 (MQ=255)
                                                                                                                   tGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGcc                                   <  1:1919009/150‑1 (MQ=255)
                                                                                                                    gATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCg                                  <  2:1821604/150‑1 (MQ=255)
                                                                                                                     aTTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                 >  1:487488/1‑150 (MQ=255)
                                                                                                                        tAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCTGGCATCCGCCGTCAATGGCTCGCCAAAATGACCGGCAATCAGCCGTTTg                              <  1:1019390/150‑1 (MQ=255)
                                                                                                                         aaTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTgg                             <  2:23560/150‑1 (MQ=255)
                                                                                                                          aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                            >  2:1370511/1‑150 (MQ=255)
                                                                                                                          aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                            <  2:95815/150‑1 (MQ=255)
                                                                                                                           ttttATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                           >  1:781048/1‑150 (MQ=255)
                                                                                                                           ttttATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                           <  2:1790547/150‑1 (MQ=255)
                                                                                                                                 ctatcACCTCATTCTGACAAGATTTAATCTTATGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCa                     <  1:1071396/150‑1 (MQ=255)
                                                                                                                                  tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTTTAAGGAGGATTCAGAAAATTAGCTATCTTTACGCCACGCATCCTACGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTctt                    >  1:1075196/1‑148 (MQ=255)
                                                                                                                                  tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATAAGCGATCTTTACGCCACGCATGCGCCGTCAAAGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                    >  1:1655950/1‑150 (MQ=255)
                                                                                                                                   atcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATg                   >  1:1438902/1‑150 (MQ=255)
                                                                                                                                   atcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGGGAGTTCATg                   <  2:1292594/150‑1 (MQ=255)
                                                                                                                                    tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCCCGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                  >  2:529271/1‑150 (MQ=255)
                                                                                                                                     cACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                 <  2:905858/150‑1 (MQ=255)
                                                                                                                                      aCCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAg                >  2:709972/1‑150 (MQ=38)
                                                                                                                                       ccTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgc               >  2:1715308/1‑150 (MQ=255)
                                                                                                                                         tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCAGAAATCAGCCGTTTGGTGAGTTCATCCAgcgg             >  1:401176/1‑150 (MQ=38)
                                                                                                                                               tGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGc       <  1:1906192/150‑1 (MQ=255)
                                                                                                                                                 aCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCa     <  2:54442/150‑1 (MQ=255)
                                                                                                                                                  cAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCGCGCATCCGCCGTCAATGCCTCGCCAAAATGAACGGCAATCAGCCGTTTGGTGAGTTCAAGCAGCGGCGATGCCAg    >  2:1254955/1‑150 (MQ=255)
                                                                                                                                                   aaGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGc   <  2:226542/150‑1 (MQ=255)
                                                                                                                                                   aaGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGACGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGc   <  2:1117051/150‑1 (MQ=255)
                                                                                                                                                   aaGATTTAATCTTTTGGCAGAATGAGGTGAATAATCTGGAAGGAGGAGTCAGAAAATTAGCGATCTTTACGCCACGCAGCCGCCGTCAATGCCTCGCCAAAAAGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCTGc   <  1:2003161/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTCCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCA  >  NC_000913/1351625‑1351922

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: