Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,249 A→T 7.0% intergenic (+224/‑139) clpB → / → PP_t07 chaperone protein ClpB/tRNA‑Asn

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947733,2490AT 7.0% 215.0 / 4.6 71intergenic (+224/‑139)clpB/PP_t07chaperone protein ClpB/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base A (29/37);  new base T (3/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 6.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

GGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGC  >  NC_002947/733103‑733392
                                                                                                                                                    |                                                                                                                                               
ggCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGt                                                                                                                                               <  2:250916/149‑1 (MQ=255)
    tttGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGt                                                                                                                                           >  1:117882/1‑149 (MQ=255)
          ccAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtgaat                                                                                                                                      >  5:285618/1‑148 (MQ=255)
           cAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgttttg                                                                                                                                    <  7:79297/149‑6 (MQ=255)
           cAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtgaatga                                                                                                                                    <  7:74447/149‑1 (MQ=255)
                 ggggTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGaaaaaaa                                                                                                                              <  7:115976/149‑1 (MQ=255)
                       aCCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGtt                                                                                                                        <  2:135713/149‑1 (MQ=255)
                       aCCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGt                                                                                                                         <  4:145663/148‑1 (MQ=255)
                         cGCCTGAAATTTTTGAAATTTTTGCTTGCAACAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTg                                                                                                                       >  1:25614/1‑148 (MQ=255)
                         cGCCTGAAATTTTTGAAATTTTTGCTGGCATCTAATTTATAGTGCCCCTAATATACGCCGCGATGTAAGGCACTAAGCGAAAAACCACCAACATCTGTGATCTCAAAACAACTTACAAATCAGAAAGATGAATGAAAAAAAGTGGTTg                                                                                                                       >  5:250361/1‑148 (MQ=255)
                                 atttttgaaatttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCaaa                                                                                                              <  1:44975/149‑1 (MQ=255)
                                         aatttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACgaagaa                                                                                                      >  4:107156/1‑149 (MQ=255)
                                             tttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGt                                                                                                   <  3:107156/148‑1 (MQ=255)
                                                   tgcATCAAATTTAGAGTGCCCCTAATATAAGCCGCGTTGTCAGTCACTAAGCGAAACACCAGCAAAATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTTAATGAAAAAAAGTGGTTGACAAGAAAAACGATTAATTTAGAATa                                                                                             <  6:250361/148‑1 (MQ=255)
                                                     cATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtg                                                                                                                                         <  6:94508/102‑1 (MQ=255)
                                                     cATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtg                                                                                                                                         >  5:94508/1‑102 (MQ=255)
                                                     cATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGcc                                                                                          >  8:129554/1‑149 (MQ=255)
                                                      aTCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCg                                                                                         >  6:161799/1‑149 (MQ=255)
                                                       tCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCgg                                                                                        >  1:317571/1‑149 (MQ=255)
                                                       tCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCgg                                                                                        >  5:229074/1‑149 (MQ=255)
                                                         aaaTTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGc                                                                                       <  2:317571/148‑1 (MQ=255)
                                                            tttAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCt                                                                                     >  2:314869/1‑147 (MQ=255)
                                                              tAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTcag                                                                                  <  6:285618/148‑1 (MQ=255)
                                                                       ccTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAAc                                                                        <  2:25614/149‑1 (MQ=255)
                                                                       ccTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAAc                                                                        >  5:263238/1‑149 (MQ=255)
                                                                                 gccgcGTAGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCg                                                               >  2:288959/1‑148 (MQ=255)
                                                                                  ccgcGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaa                                                             <  5:161799/149‑1 (MQ=255)
                                                                                  ccgcGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaa                                                             <  6:263238/149‑1 (MQ=255)
                                                                                      gTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAGACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaagaa                                                          >  1:232705/1‑149 (MQ=255)
                                                                                            aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATa                                                                                             >  4:154152/1‑107 (MQ=255)
                                                                                            aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATa                                                                                             <  3:154152/107‑1 (MQ=255)
                                                                                            aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTaaa                                                    <  7:129554/148‑1 (MQ=255)
                                                                                            aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTa                                                      >  6:5828/1‑146 (MQ=255)
                                                                                             ggCACTAAGCGAAACACCATCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAATAATGTAGAATAGCCGGCCTCAGCAGCGACAGCGCGAAAGCGAAGAAGCTAAaga                                                  <  8:266320/149‑1 (MQ=255)
                                                                                                aCTAAGCGAAACACCAGCAACATCGGAGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTc                                               >  1:4582/1‑149 (MQ=255)
                                                                                                 cTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTcc                                              <  2:176765/149‑1 (MQ=255)
                                                                                                    aGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCg                                             >  4:85060/1‑147 (MQ=255)
                                                                                                       gAAACACCAGCAACATCGGTGATCGCAAAACAACTTGACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGc                                         >  1:52061/1‑149 (MQ=255)
                                                                                                       gAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGc                                         <  2:52061/148‑1 (MQ=255)
                                                                                                            acCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAacac                                   >  8:175879/1‑149 (MQ=255)
                                                                                                              cAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAg                                 >  2:145988/1‑149 (MQ=255)
                                                                                                                gCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAg                                          <  1:244456/138‑1 (MQ=255)
                                                                                                                gCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAg                                          >  2:244456/1‑138 (MQ=255)
                                                                                                                   aCATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGGCCGAAGCGAACACAGTCTTc                            <  2:117882/149‑1 (MQ=255)
                                                                                                                     aTCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGTTTGACAAGCAAAACGAAGAATGTAGAATAGcc                                                                                          <  2:216940/85‑1 (MQ=255)
                                                                                                                     aTCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCg                           >  7:12619/1‑148 (MQ=255)
                                                                                                                     aTCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTAAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGcc                                                                                          >  1:216940/1‑85 (MQ=255)
                                                                                                                      tCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGaa                                           <  5:72650/131‑1 (MQ=255)
                                                                                                                      tCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGaa                                           >  6:72650/1‑131 (MQ=255)
                                                                                                                             tCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTa                   <  1:314869/148‑1 (MQ=255)
                                                                                                                             tCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTAc                  <  8:12619/149‑1 (MQ=255)
                                                                                                                                 tattcaacTTACTGATTTGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  1:186439/145‑1 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCg                >  5:44579/1‑147 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  7:72156/149‑1 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  6:229074/149‑1 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  1:288959/149‑1 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  2:232705/149‑1 (MQ=255)
                                                                                                                                 aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  3:85060/149‑1 (MQ=255)
                                                                                                                                  attcaacTTACTGATTTGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAag             >  6:130729/4‑149 (MQ=255)
                                                                                                                                   ttcaacTTACTGATTTGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              >  2:186439/3‑147 (MQ=255)
                                                                                                                                   ttcaacTTACTGATTTGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAagt            >  6:39828/3‑149 (MQ=255)
                                                                                                                                     caacTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAagttc          <  6:44579/149‑1 (MQ=255)
                                                                                                                                          ttACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGa                         <  3:16507/130‑1 (MQ=255)
                                                                                                                                          ttACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGa                         >  4:16507/1‑130 (MQ=255)
                                                                                                                                             aCAAATCAGTAAGTTTAATGAAAAAAAGTGGTTGACAATCAAAACGAAGAATTTAGAATAGCCTTCATCATAAGCGACAACGCGAAATCGAAGAAGCTAAAGAGTCCGAAg                                          <  6:217344/111‑1 (MQ=255)
                                                                                                                                              cAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCcg   <  5:144542/149‑1 (MQ=255)
                                                                                                                                               aaaTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaagaa                                                          <  1:52750/93‑1 (MQ=255)
                                                                                                                                               aaaTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaagaa                                                          >  2:52750/1‑93 (MQ=255)
                                                                                                                                               aaaTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCcgc  <  5:5828/149‑1 (MQ=255)
                                                                                                                                                 aTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTc                                               <  1:55693/102‑1 (MQ=255)
                                                                                                                                                 aTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTc                                               >  2:55693/1‑102 (MQ=255)
                                                                                                                                                    aGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAg                        >  6:166447/1‑122 (MQ=255)
                                                                                                                                                    aGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAg                        <  5:166447/122‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGC  >  NC_002947/733103‑733392

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