Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,142 A→T 5.4% intergenic (+27/+96) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,1420AT 5.4% 551.9 / 15.7 166intergenic (+27/+96)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (48/109);  new base T (4/5);  total (52/114)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

Only 100 of 175 total aligned reads displayed.

ACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAT  >  NC_002947/4890994‑4891291
                                                                                                                                                    |                                                                                                                                                     
aCTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGTCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCACGGGCaa                                                                                                                                                      <  2:99324/150‑1 (MQ=255)
aCTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCaa                                                                                                                                                      <  2:4832518/150‑1 (MQ=255)
aCTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCaa                                                                                                                                                      <  1:6870085/150‑1 (MQ=255)
aCTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCaa                                                                                                                                                      <  1:2052547/150‑1 (MQ=255)
     gACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCg                                                                                                                                                 <  2:511801/150‑1 (MQ=255)
      aCTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGc                                                                                                                                                <  1:1459325/150‑1 (MQ=255)
      aCTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGc                                                                                                                                                <  2:5439611/150‑1 (MQ=255)
      aCTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGc                                                                                                                                                <  1:3592567/150‑1 (MQ=255)
       cTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCt                                                                                                                                               <  1:2957927/150‑1 (MQ=255)
        tGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGcga                                                                                                                                              <  1:59369/150‑3 (MQ=255)
        tGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTc                                                                                                                                              >  2:6112116/1‑150 (MQ=255)
        tGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTc                                                                                                                                              <  2:3400687/150‑1 (MQ=255)
        tGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTc                                                                                                                                              >  2:3252259/1‑150 (MQ=255)
          gCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGcgaat                                                                                                                                            <  1:6605755/150‑5 (MQ=255)
           cAGATTTCTGCGGTGCCCGGTAAGGCCAAGCTGGCCGACAGTGATGATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGACGGGCCATTCGCGGGCAAGCCCGcgaatg                                                                                                                                           <  2:1166301/150‑6 (MQ=255)
            aGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCcac                                                                                                                                          <  2:6044695/150‑1 (MQ=255)
             gATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCcaca                                                                                                                                         <  1:1637689/150‑1 (MQ=255)
              aTTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAg                                                                                                                                        <  1:88898/150‑1 (MQ=255)
              aTTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAg                                                                                                                                        <  1:2462312/150‑1 (MQ=255)
                ttCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGa                                                                                                                                      <  2:6649918/150‑1 (MQ=255)
                ttCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGa                                                                                                                                      <  2:6649822/150‑1 (MQ=255)
                ttCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGa                                                                                                                                      <  1:515578/150‑1 (MQ=255)
                  cTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATc                                                                                                                                    <  2:1805557/150‑1 (MQ=255)
                  cTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATc                                                                                                                                    <  2:2549288/150‑1 (MQ=255)
                   tGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATca                                                                                                                                   <  2:4495165/150‑1 (MQ=255)
                   tGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATca                                                                                                                                   <  2:4495248/150‑1 (MQ=255)
                    gCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcac                                                                                                                                  <  1:28676/150‑1 (MQ=255)
                     cGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacc                                                                                                                                 <  1:4879111/150‑1 (MQ=255)
                     cGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacc                                                                                                                                 <  1:4092565/150‑1 (MQ=255)
                       gTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcaccac                                                                                                                               >  2:16226/1‑150 (MQ=255)
                        tGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACcaca                                                                                                                              <  2:5996017/150‑1 (MQ=255)
                          aCCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATg                                                                                                                            >  1:1638431/1‑150 (MQ=255)
                            cGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCTTCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCATGCCCGCTCCCTCAGGCTCACCACATGcc                                                                                                                          >  1:1972039/1‑150 (MQ=255)
                            cGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCCCATGcc                                                                                                                          >  1:1971976/1‑150 (MQ=255)
                              gTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGccag                                                                                                                        >  1:1022754/1‑148 (MQ=255)
                                aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                      <  1:253447/150‑1 (MQ=255)
                                aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                      <  2:2926746/150‑1 (MQ=255)
                                aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                      <  2:47159/150‑1 (MQ=255)
                                  ggCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAgg                                                                                                                    <  2:3430653/150‑1 (MQ=255)
                                   gCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCt                                                                                                                                               >  1:1359308/1‑122 (MQ=255)
                                   gCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCt                                                                                                                                               <  2:1359308/122‑1 (MQ=255)
                                     cAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtg                                                                                                                 <  1:2771851/150‑1 (MQ=255)
                                     cAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtg                                                                                                                 <  2:6866269/150‑1 (MQ=255)
                                     cAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtg                                                                                                                 <  2:5086686/150‑1 (MQ=255)
                                      aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                <  1:4819225/150‑1 (MQ=255)
                                       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                               <  2:3636135/150‑1 (MQ=255)
                                       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                               <  2:550791/150‑1 (MQ=255)
                                       aGCTGGCCGACAGTGATTATGACGCAGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                               <  1:6249378/150‑1 (MQ=255)
                                        ccTGGCCGACAGTCATTATGACGCGGCCGCCGAGACCATCAATCGTATCGACTAGAAAGGTTGTGAGCCTGAGGCACCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTgg                                                                                                              <  2:5609715/149‑1 (MQ=255)
                                              cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTccc                                                                                                        <  2:4809762/150‑1 (MQ=255)
                                                aCAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGcgaatg                                                                                                                                           <  1:1246331/113‑6 (MQ=255)
                                                aCAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCTTGCCCGcgaatg                                                                                                                                           >  2:1246331/1‑108 (MQ=255)
                                                 cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCTCAGGATCACCACATGGCTGTAGGTTGTGGTGTCGCtgt                                                                                                     <  2:4300191/150‑1 (MQ=255)
                                                  aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                    <  2:5939765/150‑1 (MQ=255)
                                                  aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                    <  2:4758531/150‑1 (MQ=255)
                                                   gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                   <  2:1284303/150‑1 (MQ=255)
                                                   gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                   <  2:5869882/150‑1 (MQ=255)
                                                   gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                   <  2:4737704/150‑1 (MQ=255)
                                                   gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                   <  2:5022575/150‑1 (MQ=255)
                                                      aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACcacaaccttc                                                                                                                        <  2:4242889/126‑7 (MQ=255)
                                                      aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACcacaaccttc                                                                                                                        >  1:4242889/1‑120 (MQ=255)
                                                      aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAg                                                                                                >  2:6768763/1‑150 (MQ=255)
                                                         aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                             <  1:2887440/150‑1 (MQ=255)
                                                                       gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGTGGCAGCCTGAACCATCTCGAAGGGCCATCCGCGGGCACCCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                               <  1:5922993/150‑1 (MQ=255)
                                                                        ataCCATAATTCTTTTGGTCTTGAATGGTTGTGTGCGTGTGGCAGCCTGAACCATCTCGATGGGCCTTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACTTGCCTGATGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATg                                                                              <  1:4931250/148‑1 (MQ=37)
                                                                        agaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCCCGCCCACAGGATCACCACCAGCCGGAAGGTTGTGGAGTCCCGGTGGGAGCGGGCGTGCCcgcgcagg                                                                              >  1:7053989/1‑146 (MQ=255)
                                                                            catcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGcc                                                                          >  2:5593225/1‑150 (MQ=255)
                                                                                    ttATGGTCTTGAATGGTTGTGTGCCTGTGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGTATCACCACATGCCTGACGGTTGTGGAGTCCCTGTGGGTGCGGGCTTGCCCGCGAATGGGCCAGAACAgg                                                                  <  1:3602001/149‑1 (MQ=38)
                                                                                       tGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATAACCACATGCCGGAAGGGTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCaa                                                               >  1:1661058/1‑150 (MQ=255)
                                                                                         gTCTTGAACGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATa                                                             <  2:2286474/150‑1 (MQ=255)
                                                                                         gTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATa                                                             <  2:5376436/150‑1 (MQ=255)
                                                                                         gTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATa                                                             <  1:7330280/150‑1 (MQ=255)
                                                                                            ttGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGACGGGCCTTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGACGGTTGTGGAGTCCCTGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATagag                                                          <  1:3045672/150‑1 (MQ=255)
                                                                                            ttGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCTGAAGGGTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCACTagag                                                          >  2:6218236/1‑150 (MQ=255)
                                                                                             tGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAgaga                                                         <  2:7284143/150‑1 (MQ=255)
                                                                                              gAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAt                                                        >  2:7178596/1‑150 (MQ=255)
                                                                                                 agtttGTGAGCCGGAGGCAGCCTGAACCTTCTCGCTGGGCCATTCGCGGGCAAGCCCGCGCCCACTGGATCACCACATGCCTGACGGTTGTGGAGTCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                     <  1:2296133/147‑1 (MQ=38)
                                                                                                       tgAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGc                                               >  2:197578/1‑150 (MQ=255)
                                                                                                          gCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTAc                                            >  1:7433992/1‑150 (MQ=255)
                                                                                                                     ccTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCGGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTg                                 <  2:3871931/150‑1 (MQ=255)
                                                                                                                       tGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                      >  1:46863/1‑63 (MQ=255)
                                                                                                                        gAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGccc                              <  2:2273182/150‑1 (MQ=255)
                                                                                                                            cATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAg                          <  2:1638431/150‑1 (MQ=255)
                                                                                                                            cATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAg                          <  1:6663224/150‑1 (MQ=255)
                                                                                                                             aTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGc                         <  2:7433992/150‑1 (MQ=255)
                                                                                                                               ctcGACGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGGGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGccc                       <  1:6112116/150‑1 (MQ=255)
                                                                                                                                tcGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                      >  2:4138471/1‑150 (MQ=255)
                                                                                                                                tcGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                      >  2:3998401/1‑150 (MQ=255)
                                                                                                                                  gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                    >  2:4141638/1‑150 (MQ=255)
                                                                                                                                        cggttCGCGGGCAAGCCCGCTCCCACAGGATCACCACCTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc              <  2:2653257/147‑1 (MQ=255)
                                                                                                                                        ccATCCGCGGGCAAGCCCGCTCCCACAGTATCACCACATGCCTGAAGGTTGTGGAGTCCCCGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc              <  1:1093810/150‑1 (MQ=255)
                                                                                                                                          aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                               <  2:5699929/51‑1 (MQ=255)
                                                                                                                                          aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                               >  1:5699929/1‑51 (MQ=255)
                                                                                                                                          aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGa            <  2:6103095/150‑1 (MQ=255)
                                                                                                                                           ggagcgGGCTTGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                     <  2:2902867/57‑1 (MQ=25)
                                                                                                                                           ggagcgGGCTTGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                     >  1:2902867/4‑60 (MQ=25)
                                                                                                                                                  gCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTg    >  2:1117070/1‑150 (MQ=17)
                                                                                                                                                    aaGCCCGCTCCCCCAGTATCACCAGATGCCTTATGGTTGTGGTGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt  <  1:4360879/150‑1 (MQ=255)
                                                                                                                                                    aaGCCCGCTCCCACTGCACCACCACCGGCCTGAAGTTTGTGGAGTCCCTGTGGGTGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt  <  1:4360909/150‑1 (MQ=255)
                                                                                                                                                    aaGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt  >  2:6555453/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACTTCGACTGGCAGATTTCTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAT  >  NC_002947/4890994‑4891291

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