Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 733,249 A→T 7.5% 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.5% 196.3 / 5.4 67intergenic (+224/‑139)clpB/PP_t07chaperone protein ClpB/tRNA‑Asn
Reads supporting (aligned to +/- strand):  ref base A (32/30);  new base T (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.63e-01

GGGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGCGATA  >  NC_002947/733102‑733396
                                                                                                                                                   |                                                                                                                                                   
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                                                                                                                                         <  6:169677/149‑5 (MQ=255)
           ccAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACTGATTTGTAAGTtgtttt                                                                                                                                         >  7:125791/1‑145 (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                                                                                                                              <  7:178667/143‑1 (MQ=255)
                            ccTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTgg                                                                                                                              >  8:178667/1‑143 (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                                                                                            >  5:230841/1‑131 (MQ=255)
                                                                          tAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATATCCgg                                                                                            <  6:230841/131‑1 (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                                                                                  <  7:75045/108‑1 (MQ=255)
                                                                                                           acacCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCagcagc                                                                                  >  8:75045/1‑108 (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                                   <  8:13820/140‑1 (MQ=255)
                                                                                                                          gTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTc                                   >  7:13820/1‑140 (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)
                                                                                                                                                  cAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGCGa    >  3:254745/1‑147 (MQ=255)
                                                                                                                                                   aGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGCGATa  <  2:162050/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GGGCCTTTGGGCCAACCCGGGGTAACCGCCTGAAATTTTTGAAATTTTTGCTTGCATCAAATTTAGAGTGCCCCTAATATACGCCGCGTTGTCAGGCACTAAGCGAAACACCAGCAACATCGGTGATCGCAAAACAACTTACAAATCAGTAAGTTGAATGAAAAAAAGTGGTTGACAAGCAAAACGAAGAATGTAGAATAGCCGGCCTCAGCAGCGACAACGCGAAAGCGAAGAAGCTAAAGAGTCCGAAGCGAACACAGTCTTCGGAGTTGTACCGAAGTTCAGTTCCGCGATA  >  NC_002947/733102‑733396

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: