Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,244 A→T 8.6% intergenic (+219/‑144) 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,2440AT 8.6% 193.3 / 4.5 66intergenic (+219/‑144)clpB/PP_t07chaperone protein ClpB/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base A (30/30);  new base T (2/4);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGC  >  NC_002947/733102‑733392
                                                                                                                                              |                                                                                                                                                    
gggCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCa                                                                                                                                                 <  3:48137/148‑1 (MQ=255)
 ggCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAg                                                                                                                                                >  6:313411/1‑148 (MQ=255)
    cTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAg                                                                                                                                            >  3:299720/1‑149 (MQ=255)
         ggCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgtt                                                                                                                                       <  6:146742/149‑3 (MQ=255)
           ccAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgtttt                                                                                                                                     >  7:125791/1‑145 (MQ=255)
           ccAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgtttt                                                                                                                                     <  6:169677/149‑5 (MQ=255)
            cAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtgaatg                                                                                                                                     >  3:255975/1‑148 (MQ=255)
             aaCCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgttttgc                                                                                                                                   <  8:125791/149‑7 (MQ=255)
                  ggggTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGaaaaaaa                                                                                                                              <  6:182080/149‑1 (MQ=255)
                  ggggTAACCGCCTGAAATTTTTGAAATTTTAGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGaaaaaaa                                                                                                                              <  4:47370/149‑1 (MQ=255)
                     gcaaCCCCCTCAACTTTAAGAAATTTTTTCTTCCACAACATTCAGAGTCCCCCTAATATACGCCGCGCTGTCAGGCACTAACCCAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTg                                                                                                                           <  4:255975/147‑1 (MQ=255)
                        aCCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGt                                                                                                                         <  3:227272/148‑1 (MQ=255)
                         ccGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGtt                                                                                                                        >  6:63073/1‑148 (MQ=255)
                          cGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGa                                                                                                                      >  5:269607/1‑149 (MQ=255)
                            ccTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTgg                                                                                                                          >  8:178667/1‑143 (MQ=255)
                            ccTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTgg                                                                                                                          <  7:178667/143‑1 (MQ=255)
                              ttaaatttttgaaatttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAg                                                                                                                  <  5:159651/147‑1 (MQ=255)
                              tgaaatttttgaaattttTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtgaat                                                                                                                                      >  3:300352/1‑129 (MQ=255)
                              tgaaatttttgaaattttTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTtgaat                                                                                                                                      <  4:300352/129‑1 (MQ=255)
                                aaatttttgaaatttttgCTTTCATAAAATTTAAAGTGCCCATAATAAACGCCGCGTTGTCTGCAAATAAGCGAAAAACCAGAAACTTGGGTCATCGAATAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAg                                                                                                                  <  1:75527/147‑1 (MQ=255)
                                   tttttgaaatttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGaaaaaaa                                                                                                                              <  5:185248/132‑1 (MQ=255)
                                   tttttgaaatttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGaaaaaaa                                                                                                                              >  6:185248/1‑132 (MQ=255)
                                            tttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATg                                                                                                    <  5:313411/149‑1 (MQ=255)
                                            tttttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATg                                                                                                    >  2:207297/1‑149 (MQ=255)
                                              tttgCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGt                                                                                                   <  6:269607/148‑1 (MQ=255)
                                                      cATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGcc                                                                                          >  7:56839/1‑149 (MQ=255)
                                                           aaTTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTg                                                                                                                       >  8:69243/1‑115 (MQ=255)
                                                           aaTTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTg                                                                                                                       <  7:69243/115‑1 (MQ=255)
                                                                 gagTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTcagca                                                                                <  5:131572/148‑1 (MQ=255)
                                                                        ccTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGcc                                                                                      <  4:273040/135‑1 (MQ=255)
                                                                        ccTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGcc                                                                                      >  3:273040/1‑135 (MQ=255)
                                                                          tAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATATCCgg                                                                                        <  6:230841/131‑1 (MQ=255)
                                                                          tAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATATCCgg                                                                                        >  5:230841/1‑131 (MQ=255)
                                                                                 cgccgcGTTGTCAGGCACTAAGCGAAACAACAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCg                                                               <  5:63073/149‑1 (MQ=255)
                                                                                   ccgcGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAATGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaa                                                             <  1:250635/149‑1 (MQ=255)
                                                                                    ggcgTTGTCAGGAACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAAAAACTTACTAATCAGTAAGTTGAATGAAAATAATTTGTTGAAAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaa                                                             <  6:285307/147‑1 (MQ=255)
                                                                                     gcgTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaaga                                                           <  3:145321/149‑1 (MQ=255)
                                                                                     gcgTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCgaag                                                            >  1:146707/1‑148 (MQ=255)
                                                                                           tCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTa                                                      >  6:258515/1‑148 (MQ=255)
                                                                                             aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTaaa                                                    <  4:299720/148‑1 (MQ=255)
                                                                                             aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTaaa                                                    <  3:54631/148‑1 (MQ=255)
                                                                                             aGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTa                                                      >  4:54631/1‑146 (MQ=255)
                                                                                                   tAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAATGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAaga                                                  <  5:216473/144‑1 (MQ=255)
                                                                                                   tAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTCGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTcc                                              >  8:44514/1‑148 (MQ=255)
                                                                                                   tAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAaga                                                  >  6:216473/1‑144 (MQ=255)
                                                                                                      gCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAg                                          <  2:146707/149‑1 (MQ=255)
                                                                                                           acacCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCagcagc                                                                              >  8:75045/1‑108 (MQ=255)
                                                                                                           acacCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCagcagc                                                                              <  7:75045/108‑1 (MQ=255)
                                                                                                              ccAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGAGACAACGCGAAAGAGAAGAAGCTAAAGAGTCAGAAGCGAAcaca                                  >  5:262620/1‑149 (MQ=255)
                                                                                                               cAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCt                                                       >  4:111608/1‑127 (MQ=255)
                                                                                                               cAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCt                                                       <  3:111608/127‑1 (MQ=255)
                                                                                                                aGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGt                                >  1:162050/1‑149 (MQ=255)
                                                                                                                         ggTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAg                        <  5:258515/148‑1 (MQ=255)
                                                                                                                          gTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTc                               >  7:13820/1‑140 (MQ=255)
                                                                                                                          gTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTc                               <  8:13820/140‑1 (MQ=255)
                                                                                                                                gCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTAcc                 >  3:213900/1‑148 (MQ=255)
                                                                                                                                 cAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCg                >  6:278579/1‑148 (MQ=255)
                                                                                                                                  aaaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGaa              <  8:56839/149‑1 (MQ=255)
                                                                                                                                   attcaacTTACTGATCTGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAag             >  1:8732/4‑149 (MQ=255)
                                                                                                                                   aaaCAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAag             >  2:180447/1‑149 (MQ=255)
                                                                                                                                    aacaacTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAagt            >  2:3482/1‑149 (MQ=255)
                                                                                                                                      ctacTTAGAAATCAGTAAGTTGAATGAAAAATAGTGGTTGACAAGCAAAACGAATAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAg                        <  8:140198/133‑1 (MQ=255)
                                                                                                                                      caacTTACAAATCAGTAAGTTGAATCAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACTCGAAAGCGAAGAAGATAAAGAGTCCGAACCGAACACAGTCTTCGGAg                        >  7:140198/1‑135 (MQ=255)
                                                                                                                                       aacTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAAcaca                                  <  7:257656/124‑1 (MQ=255)
                                                                                                                                       aacTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAagtt           >  6:155944/1‑147 (MQ=255)
                                                                                                                                          ttACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAagttcagtt      >  2:32834/1‑149 (MQ=255)
                                                                                                                                              aaaTCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCcgc  >  3:180681/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                    
GGGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGC  >  NC_002947/733102‑733392

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: