Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,777 T→A 6.5% intergenic (‑17/+51) 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,7770TA 6.5% 453.1 / 16.4 154intergenic (‑17/+51)ycjD/sapFDUF559 family endonuclease‑related protein/antimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (70/73);  new base A (4/6);  total (75/79)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.49e-01

Only 100 of 158 total aligned reads displayed.

CCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351629‑1351926
                                                                                                                                                    |                                                                                                                                                     
ccccACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGAttt                                                                                                                                                      >  2:737491/1‑150 (MQ=255)
  ccACTGGATGTTGACGCCGAAATTTGAAAACACTAAAACGTCGGCAGCGAAGATATCGCCAGAGCTTTCGTACCTGCAAAGTGAGATTGCGTCTTAAATAACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTaa                                                                                                                                                    >  2:143006/1‑150 (MQ=255)
   cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                                   >  1:861076/1‑150 (MQ=255)
   cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                                   <  2:232877/150‑1 (MQ=255)
    aCTTGATGTTGACGGCGAAATTTTAAATCACTAAAACGTCGGCTGCGAAGATATCGACAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTAGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATc                                                                                                                                                  <  1:856696/150‑1 (MQ=255)
     cTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCt                                                                                                                                                 >  2:2037013/1‑150 (MQ=255)
         aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCTGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGAATGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                             <  1:346820/150‑1 (MQ=255)
         aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                             >  1:1114469/1‑150 (MQ=255)
         aTGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTg                                                                                                                                             <  2:1941621/150‑1 (MQ=255)
          tGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                                                                                                                            >  2:1734400/1‑150 (MQ=255)
          tGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGt                                                                                                                                            <  1:682061/150‑1 (MQ=255)
            ttGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCa                                                                                                                                          <  2:1551047/150‑1 (MQ=255)
             tGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAg                                                                                                                                         <  1:1625610/150‑1 (MQ=255)
             tGACGGCGAAATTTGAAATCACTAAAACGTCGGCAGCGAAGATATCGCCAGAGATTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCAGTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATGTAATCTTTTGTCat                                                                                                                                         >  2:258702/1‑149 (MQ=255)
              gACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGa                                                                                                                                        >  1:152651/1‑150 (MQ=255)
                cGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAAt                                                                                                                                      >  2:162234/1‑150 (MQ=255)
                 ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATg                                                                                                                                     <  1:632232/150‑1 (MQ=38)
                 ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATACCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATg                                                                                                                                     >  2:1889746/1‑150 (MQ=255)
                  gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGa                                                                                                                                    <  2:601074/150‑1 (MQ=38)
                   cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                   >  1:1465339/1‑150 (MQ=255)
                      aaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGATTGAGGTg                                                                                                                                >  1:1858481/1‑150 (MQ=255)
                              aTCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCt                                                                                                                        <  2:1780680/150‑1 (MQ=255)
                               tCACTAAAACGTCGGCTGAGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTg                                                                                                                       <  2:1021991/150‑1 (MQ=255)
                                cACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTgg                                                                                                                      >  2:193995/1‑150 (MQ=255)
                                cACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTgg                                                                                                                      >  2:1095446/1‑150 (MQ=255)
                                 aCTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGa                                                                                                                     >  1:23560/1‑150 (MQ=255)
                                    aaaaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAagg                                                                                                                  >  1:780772/1‑150 (MQ=255)
                                        cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                              >  2:1708590/1‑150 (MQ=255)
                                        cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                              <  2:842091/150‑1 (MQ=255)
                                           cGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGCCAAGATTTAATCTTTTGTCAGAAGGAGGTGAATAATCTGGAAGGAGGATTc                                                                                                           >  2:1272552/1‑150 (MQ=255)
                                               tGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTAGTCAGAATGAGGTGAATAAACTGGAAGTAGGATTCAGaa                                                                                                       >  2:1434735/1‑150 (MQ=255)
                                                 cGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                     >  1:226542/1‑150 (MQ=255)
                                                 cGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGAGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                     <  1:2038260/150‑1 (MQ=255)
                                                  gAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCCTTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAAt                                                                                                    <  2:1844173/150‑1 (MQ=255)
                                                     gATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAg                                                                                                 <  2:1051477/150‑1 (MQ=255)
                                                      atatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGc                                                                                                >  1:596455/1‑150 (MQ=255)
                                                       tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                               <  1:1435563/150‑1 (MQ=255)
                                                        atCGCGAGAGCTTTCGTTCCTGCAATGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTTGAAGGAGGATTCAGAAAATTAGCGa                                                                                              >  1:140988/1‑150 (MQ=255)
                                                         tCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGAt                                                                                             >  1:565767/1‑150 (MQ=255)
                                                          cGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAATTCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCCGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGATGATTCAGAAAATTAGCGAta                                                                                            >  2:406323/1‑149 (MQ=255)
                                                          cGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATc                                                                                            >  2:372054/1‑150 (MQ=255)
                                                            ccAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCtt                                                                                          >  1:1307521/1‑150 (MQ=255)
                                                             cAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAACTCACGCGCATTTGAGTTAATTTTATCTATCAACTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTCGAATGAGGATTCAGAAAAATAGCGATCttt                                                                                         >  1:178802/1‑150 (MQ=255)
                                                             cAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCttt                                                                                         <  1:1480449/150‑1 (MQ=255)
                                                                     tCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACg                                                                                 <  2:1314705/150‑1 (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:1891495/150‑1 (MQ=255)
                                                                                      gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                                <  2:1212368/150‑1 (MQ=255)
                                                                                         tGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCg                                                             >  2:1479682/1‑150 (MQ=255)
                                                                                          gCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGc                                                            >  1:916782/1‑150 (MQ=255)
                                                                                           cGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGcc                                                           <  1:1034447/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)
                                                                                                 aaaTCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATg                                                     <  1:436148/150‑1 (MQ=255)
                                                                                                   aTCACGCTCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCATAATGATGCGATTAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACTCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAc                                                   <  2:1026306/150‑1 (MQ=255)
                                                                                                   aTCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAATATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGAc                                                   <  1:411059/150‑1 (MQ=255)
                                                                                                      actcgcATTTGATTTAATATTATCTATCACCTCATTCTGACAAGATTTAATCTTTTTTCAGAATGAGGTGAATAATCTTGACGGAGGATTCAGAAAATTAGCGATCTTTCCGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCgg                                                <  2:161884/147‑1 (MQ=255)
                                                                                                      aCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCgg                                                <  2:494214/150‑1 (MQ=255)
                                                                                                       cgcgGATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  1:606244/150‑1 (MQ=255)
                                                                                                       cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  1:759748/150‑1 (MQ=255)
                                                                                                       cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               >  2:1630004/1‑150 (MQ=255)
                                                                                                       cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  2:1639765/150‑1 (MQ=255)
                                                                                                       cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                               <  2:1647837/150‑1 (MQ=255)
                                                                                                                gATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCg                                      <  2:1821604/150‑1 (MQ=255)
                                                                                                                 aTTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGt                                     >  1:487488/1‑150 (MQ=255)
                                                                                                                  tttAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTAATAATGAGTTGAATAATATGGAAGGAGGATTCAGAAAATTAGCGATCTTTACTCCACGCATACGCCGTCAATGCCTCGCCAAAATGACCGGCAATCATCCGtt                                    >  1:312460/1‑150 (MQ=255)
                                                                                                                    tAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCTGGCATCCGCCGTCAATGGCTCGCCAAAATGACCGGCAATCAGCCGTTTg                                  <  1:1019390/150‑1 (MQ=255)
                                                                                                                      aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                                <  2:95815/150‑1 (MQ=255)
                                                                                                                      aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                                >  2:1726699/1‑150 (MQ=255)
                                                                                                                       ttttATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTg                               <  2:1790547/150‑1 (MQ=255)
                                                                                                                            tctatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTc                          <  2:1401044/150‑1 (MQ=255)
                                                                                                                             ctatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCa                         >  1:440924/1‑150 (MQ=255)
                                                                                                                              tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTTTAAGGAGGATTCAGAAAATTAGCTATCTTTACGCCACGCATCCTACGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTctt                        >  1:1075196/1‑148 (MQ=255)
                                                                                                                              tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                        <  1:193995/150‑1 (MQ=255)
                                                                                                                               atcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGGGAGTTCATg                       <  2:1292594/150‑1 (MQ=255)
                                                                                                                                tcACCTCATTCTGACAAGATTTAATCTTTTGTGAGAATGAGGTGAATAATCTGGATGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCATTGCCTCGCCAAAATGACCGTCAATCAGCCGTTTGGTGAGTTCATGc                      <  1:1740949/150‑1 (MQ=255)
                                                                                                                                tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCCCGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                      >  2:529271/1‑150 (MQ=255)
                                                                                                                                tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                      <  1:1488476/150‑1 (MQ=255)
                                                                                                                                 cACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                     <  2:905858/150‑1 (MQ=255)
                                                                                                                                  aCCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAg                    >  2:709972/1‑150 (MQ=38)
                                                                                                                                   ccTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgc                   >  1:161912/1‑150 (MQ=38)
                                                                                                                                     tCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcgg                 <  1:1218626/150‑1 (MQ=255)
                                                                                                                                     tCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCAGAAATCAGCCGTTTGGTGAGTTCATCCAgcgg                 >  1:401176/1‑150 (MQ=38)
                                                                                                                                       aTTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgcggcg               <  1:1889747/150‑1 (MQ=255)
                                                                                                                                           tGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGc           <  1:1906192/150‑1 (MQ=255)
                                                                                                                                           tGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGc           >  1:1576713/1‑150 (MQ=255)
                                                                                                                                             aCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCa         <  2:54442/150‑1 (MQ=255)
                                                                                                                                             aCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCa         <  2:487488/150‑1 (MQ=255)
                                                                                                                                              cAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCGCGCATCCGCCGTCAATGCCTCGCCAAAATGAACGGCAATCAGCCGTTTGGTGAGTTCAAGCAGCGGCGATGCCAg        >  2:1254955/1‑150 (MQ=255)
                                                                                                                                              cAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAg        >  2:1683289/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)
                                                                                                                                                aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca      <  1:780105/150‑1 (MQ=255)
                                                                                                                                                    ttAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCAGCAGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCCTGCAGCGGCTATGACAGCACATc  >  2:1844592/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCCACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATC  >  NC_000913/1351629‑1351926

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: