Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,703,294:1 +T 14.1% intergenic (+101/‑66) lysS → / → PP_1497 lysine‑‑tRNA ligase/TetR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,703,2941.T14.1% 238.9 / 42.9 78intergenic (+101/‑66)lysS/PP_1497lysine‑‑tRNA ligase/TetR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (34/33);  new base T (5/6);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-02

CGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAA  >  NC_002947/1703147‑1703439
                                                                                                                                                    |                                                                                                                                                 
cGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCtt                                                                                                                                                 <  1:1545836/150‑1 (MQ=255)
cGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCtt                                                                                                                                                 <  1:827402/150‑1 (MQ=255)
   tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTa                                                                                                                                              <  1:434529/150‑1 (MQ=255)
   tCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTCGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttgta                                                                                                                                             >  2:1700977/1‑146 (MQ=37)
   ggatcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTa                                                                                                                                              >  2:1274529/5‑150 (MQ=255)
      gcgcCGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGa                                                                                                                                           <  2:389770/150‑1 (MQ=38)
      gcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGa                                                                                                                                           <  1:747361/150‑1 (MQ=255)
      gcgACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGa                                                                                                                                           <  2:1248668/150‑1 (MQ=255)
      gcgACGTCATCCTGTTCACGCACATCCGCCCGCAGGCCTGCGTGAGCTGCACGCGCCACATTTCGACGCTGCTTTTTCTTGCCTGTAGCTTTATGTTGTCAGTACCGGCCGCTTCGCGGGCAAACCCGCTTCTACAAGGCCTTTGTAGGa                                                                                                                                           <  2:1403378/150‑1 (MQ=25)
        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                         <  2:842175/150‑1 (MQ=255)
        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGc                                                                                                                                         <  1:1274529/150‑1 (MQ=255)
        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGCCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                         <  2:339329/150‑1 (MQ=255)
        gACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGCCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTACACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                         <  1:349381/150‑1 (MQ=38)
          cGTCATCCTGTTCCCGCGCATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCgg                                                                                                                                       <  2:263615/150‑1 (MQ=38)
          cGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                       <  2:921956/150‑1 (MQ=255)
          cGTCATCCGGTTCCCGCCCATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCCCCTGTTGGGGCGGCTTTTTCTTGCCTGAAGCTGTATGTTGTCAGTACCGGCCTCTTCGCGGGTCACCCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                       <  2:198188/150‑1 (MQ=255)
          cGTCATCCGGTTCCCGCACATGCGTCCGCAGGCCGGGGTGAGCTGAACGAGCCGCCTTTCGAGGCGCCTTTTTCTCGCCTGAAGCTTTATGTTGTCAGTACCGGCCGCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCgg                                                                                                                                       <  2:1409932/150‑1 (MQ=39)
           gTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCggg                                                                                                                                      <  2:792625/150‑1 (MQ=255)
            tCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGt                                                                                                                                     >  1:1771002/1‑150 (MQ=255)
              aTCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                         >  1:3533/1‑144 (MQ=255)
              aTCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGc                                                                                                                                         <  2:3533/144‑1 (MQ=255)
               tCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCCCATACAAAGGCCTTGTAGTAGCGGGTTTa                                                                                                                                  >  1:1744510/1‑150 (MQ=255)
                 cTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAcc                                                                                                                                <  2:1586655/150‑1 (MQ=255)
                  tGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAccc                                                                                                                               <  2:181519/150‑1 (MQ=255)
                    ttCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCcgc                                                                                                                             <  1:1642219/150‑1 (MQ=255)
                          cacaTGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAggg                                                                                                                       >  2:558173/1‑150 (MQ=255)
                             aTGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCtt                                                                                                                                                 >  1:389092/1‑121 (MQ=255)
                             aTGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCtt                                                                                                                                                 <  2:389092/121‑1 (MQ=255)
                                         gCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAcc                                                                                                                                >  1:782165/1‑126 (MQ=255)
                                         gCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTAcc                                                                                                                                <  2:782165/126‑1 (MQ=255)
                                           cTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAac                                                                                                      >  1:648008/1‑150 (MQ=255)
                                            tgagtgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTCGTCAGTGCTGGCCCCTTCGCGGGTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGAcaacata                                                                                                     <  1:1700977/150‑7 (MQ=14)
                                               gtgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACAccc                                                                                                  <  1:558173/150‑1 (MQ=255)
                                                tgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCa                                                                                                                   <  2:1237179/132‑1 (MQ=255)
                                                tgagCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCa                                                                                                                   >  1:1237179/1‑132 (MQ=255)
                                                    cTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGc                                                                                                                 <  1:606777/130‑1 (MQ=255)
                                                    cTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGc                                                                                                                 >  2:606777/1‑130 (MQ=255)
                                                      gAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAg                                                                                           <  2:1771002/150‑1 (MQ=255)
                                                        aCGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCCGCACTGACGAAACACCCCCAtaggga                                                                                         >  2:1644691/1‑145 (MQ=255)
                                                               gccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCAAGAGGTATccccc                                                                                  >  2:283315/1‑150 (MQ=255)
                                                                ccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGa                                                                                                            <  2:1693299/123‑1 (MQ=255)
                                                                ccTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGa                                                                                                            >  1:1693299/1‑123 (MQ=255)
                                                                        gggggttttttGTTGCCTGAGGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACAccc                                                                         <  2:1744510/144‑1 (MQ=255)
                                                                          cggcTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTGCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGc                                                                       <  1:1049178/150‑1 (MQ=255)
                                                                                    ttGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                             >  2:583925/1‑150 (MQ=255)
                                                                                    ttGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGAc                                                             >  2:955216/1‑150 (MQ=255)
                                                                                          gAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGc                                                                                                                                         <  2:1244835/68‑1 (MQ=255)
                                                                                          gAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGc                                                                                                                                         >  1:1244835/1‑68 (MQ=255)
                                                                                             gCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAg                                                         >  2:1408216/1‑145 (MQ=255)
                                                                                             gCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAg                                                         <  1:1408216/145‑1 (MQ=255)
                                                                                               tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGcgccgc                                                  >  1:750536/1‑150 (MQ=255)
                                                                                               tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttg                                                                                                                                               >  2:91312/1‑54 (MQ=17)
                                                                                               tttATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAGGCCtttg                                                                                                                                               <  1:91312/56‑3 (MQ=17)
                                                                                                  aTGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCgg                                               <  1:1217596/150‑1 (MQ=255)
                                                                                                                 gCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCggg                                                                                                                                      >  1:420055/1‑48 (MQ=255)
                                                                                                                 gCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCggg                                                                                                                                      <  2:420055/48‑1 (MQ=255)
                                                                                                                            ggTAAACCCGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCg                     >  1:1405379/1‑150 (MQ=255)
                                                                                                                             gTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACAccc                                                                                                  >  1:370844/1‑72 (MQ=255)
                                                                                                                             gTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACAccc                                                                                                  <  2:370844/72‑1 (MQ=255)
                                                                                                                             gTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGt                    <  1:283315/150‑1 (MQ=255)
                                                                                                                                 acacGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCGTGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                >  2:1034756/4‑150 (MQ=255)
                                                                                                                                 aCCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                >  1:1492839/1‑150 (MQ=255)
                                                                                                                                 aCCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCaa                <  1:955216/150‑1 (MQ=255)
                                                                                                                                  cccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAAt               >  1:175808/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa              >  2:189614/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa              >  2:1629118/1‑150 (MQ=255)
                                                                                                                                   ccGCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATa              >  2:1096575/1‑150 (MQ=255)
                                                                                                                                     gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                                                       <  1:340791/59‑1 (MQ=14)
                                                                                                                                     gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGaaa                                                                                                       >  2:340791/1‑59 (MQ=14)
                                                                                                                                     gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGa                                                                                            >  2:1456071/1‑70 (MQ=18)
                                                                                                                                     gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGa                                                                                            <  1:1456071/70‑1 (MQ=18)
                                                                                                                                     gCTCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATAcc            >  2:1771977/1‑150 (MQ=37)
                                                                                                                                       tCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCaa          >  2:799228/1‑150 (MQ=255)
                                                                                                                                       tCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCaa          >  1:1310675/1‑150 (MQ=255)
                                                                                                                                       tCCTACAAGGCCTTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCAGGAAAAGCAAACGCCGTGGCCGAAAGGGTGGAATAACaa          >  2:495154/1‑150 (MQ=25)
                                                                                                                                           aCAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGac                                                                             >  1:813787/1‑79 (MQ=255)
                                                                                                                                           aCAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGac                                                                             <  2:813787/79‑1 (MQ=255)
                                                                                                                                           aCAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGcc      >  1:817119/1‑150 (MQ=255)
                                                                                                                                             aaaGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGc       >  1:1174050/1‑147 (MQ=255)
                                                                                                                                             aaaGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGc       <  2:1174050/147‑1 (MQ=255)
                                                                                                                                             aaaGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGc    >  2:783035/1‑150 (MQ=255)
                                                                                                                                             aaaGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGc    >  2:1444944/1‑150 (MQ=255)
                                                                                                                                               aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGCTTAGTGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCaa  >  1:1261448/1‑150 (MQ=255)
                                                                                                                                               aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCaa  <  2:648008/150‑1 (MQ=255)
                                                                                                                                               aGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCaa  >  1:517686/1‑150 (MQ=255)
                                                                                                                                                  ccTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAg                                                         >  1:1582866/1‑92 (MQ=255)
                                                                                                                                                  ccTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAg                                                         <  2:1582866/92‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGATCCGCGACGTCATCCTGTTCCCGCACATGCGTCCGCAGGCCTGAGTGAGCTGAACGAGCCGCCTTTCGAGGCGGCTTTTTCTTGCCTGAAGCTTTATGTTGTCAGTACCGGCCTCTTCGCGGGTAAACCCGCTCCTACAAAGGCCTTGTAGGAGCGGGTTTACCCGCGAAGGGGCCAGCACTGACGAAACACCCCCATGAGGTATCCCCCGTGACACCCGCCATGCCTGACCAAGGCGCCGCCGGCATCGCCACCGCCGTCGCCGAAAGCGTGCAATACCAAGGCCGCAA  >  NC_002947/1703147‑1703439

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: