Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,351,781 T→A 6.4% intergenic (‑21/+47) 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,7810TA 6.4% 466.3 / 15.4 155intergenic (‑21/+47)ycjD/sapFDUF559 family endonuclease‑related protein/antimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (74/71);  new base A (4/6);  total (78/77)
Fisher's exact test for biased strand distribution p-value = 5.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 159 total aligned reads displayed.

CACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCG  >  NC_000913/1351632‑1351930
                                                                                                                                                     |                                                                                                                                                     
cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                                       <  2:232877/150‑1 (MQ=255)
cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                                       >  1:861076/1‑150 (MQ=255)
cACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTGCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAAt                                                                                                                                                       <  2:1814349/150‑1 (MQ=255)
 aCTTGATGTTGACGGCGAAATTTTAAATCACTAAAACGTCGGCTGCGAAGATATCGACAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTAGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATc                                                                                                                                                      <  1:856696/150‑1 (MQ=255)
  cTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCt                                                                                                                                                     >  2:2037013/1‑150 (MQ=255)
      aTGTTGTCGGCGAAATTTGAAATCACTAACACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGCCAAGATTTAATCTTTTg                                                                                                                                                 <  2:1307610/150‑1 (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)
        gTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTc                                                                                                                                               <  2:1607766/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                                                                                                                                            >  2:785753/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)
              ggCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATACCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATg                                                                                                                                         >  2:1889746/1‑150 (MQ=255)
               gCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGa                                                                                                                                        <  2:601074/150‑1 (MQ=38)
                cGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAg                                                                                                                                       >  1:1465339/1‑150 (MQ=255)
                  aaaTTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGt                                                                                                                                     <  1:2023128/150‑1 (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)
                                 aaaaCGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAagg                                                                                                                      >  1:780772/1‑150 (MQ=255)
                                     cGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAggagga                                                                                                                  >  2:1708590/1‑150 (MQ=255)
                                        cGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGCCAAGATTTAATCTTTTGTCAGAAGGAGGTGAATAATCTGGAAGGAGGATTc                                                                                                               >  2:1272552/1‑150 (MQ=255)
                                              cGAAGATTTCGCCAGAGGTTTCGTTCCTGCATTGTGAGTTTGCGTCTTAAATCACGCGCATTTGATTTAAATTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                         <  1:187858/150‑1 (MQ=255)
                                              cGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGaaaa                                                                                                         >  1:226542/1‑150 (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)
                                                    tatCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCg                                                                                                   >  2:1488475/1‑150 (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)
                                                        gCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTCAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCt                                                                                               >  1:154699/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)
                                                                   cGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGc                                                                                    >  2:414764/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)
                                                                                  agaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGc                                                                     <  2:1629377/150‑1 (MQ=255)
                                                                                   gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                                    <  2:1212368/150‑1 (MQ=255)
                                                                                   gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                                    <  1:902794/150‑1 (MQ=255)
                                                                                   gaTTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGcc                                                                    <  2:219800/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)
                                                                                           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)
                                                                                                   aCGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCgg                                                    <  2:494214/150‑1 (MQ=255)
                                                                                                    cgcgGATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                                   <  1:606244/150‑1 (MQ=255)
                                                                                                    cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                                   <  1:759748/150‑1 (MQ=255)
                                                                                                    cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                                   <  2:1639765/150‑1 (MQ=255)
                                                                                                    cgcgCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGc                                                   <  2:1647837/150‑1 (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)
                                                                                                                   aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                                    <  2:95815/150‑1 (MQ=255)
                                                                                                                   aTTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGt                                    >  2:1434532/1‑150 (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)
                                                                                                                           tatcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATAAGCGATCTTTACGCCACGCATGCGCCGTCAAAGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCAt                            >  1:1655950/1‑150 (MQ=255)
                                                                                                                             tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCCCGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                          >  2:529271/1‑150 (MQ=255)
                                                                                                                             tcACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGc                          <  1:1488476/150‑1 (MQ=255)
                                                                                                                              cACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                         <  2:905858/150‑1 (MQ=255)
                                                                                                                              cACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCa                         >  1:234693/1‑150 (MQ=255)
                                                                                                                               aCCTCATTCTGACAAAAGATTAAATCTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAg                        >  2:709972/1‑150 (MQ=38)
                                                                                                                                ccTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAgc                       >  2:1715308/1‑150 (MQ=255)
                                                                                                                                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:1576713/1‑150 (MQ=255)
                                                                                                                                          aCAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCa             <  2:54442/150‑1 (MQ=255)
                                                                                                                                           cAAGATTTAATCTTTTGTCAGAATGAGTTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCGCGCATCCGCCGTCAATGCCTCGCCAAAATGAACGGCAATCAGCCGTTTGGTGAGTTCAAGCAGCGGCGATGCCAg            >  2:1254955/1‑150 (MQ=255)
                                                                                                                                            aaGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGACGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGc           <  2:1117051/150‑1 (MQ=255)
                                                                                                                                             aGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGca          <  1:780105/150‑1 (MQ=255)
                                                                                                                                                tttAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTTAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATAAGCCGTTTGGTGTGTTCATGAAGCGGCGATGCCAGaacaa       >  2:166722/1‑149 (MQ=255)
                                                                                                                                                  tAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATcc     >  2:885230/1‑150 (MQ=255)
                                                                                                                                                    aTCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGAGAGTTCATGCAGCGGCGATGCCAGCACATCcgc   >  1:237715/1‑150 (MQ=255)
                                                                                                                                                    aTCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTGTGGGGAGTTCATGCAGCGGCGATGCCAGCACATCcgc   >  2:1754294/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CACTGGATGTTGACGGCGAAATTTGAAATCACTAAAACGTCGGCTGCGAAGATATCGCCAGAGCTTTCGTTCCTGCAAAGTGAGATTGCGTCTTAAATCACGCGCATTTGATTTAATTTTATCTATCACCTCATTCTGACAAGATTTAATCTTTTGTCAGAATGAGGTGAATAATCTGGAAGGAGGATTCAGAAAATTAGCGATCTTTACGCCACGCATCCGCCGTCAATGCCTCGCCAAAATGACCGGCAATCAGCCGTTTGGTGAGTTCATGCAGCGGCGATGCCAGCACATCCGCG  >  NC_000913/1351632‑1351930

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: