Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,041 Δ1 bp 10.0% intergenic (+130/+59) PP_5440 → / ← PP_t72 hypothetical protein/tRNA‑Phe

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947344,0410A.10.0% 181.7 / 14.6 60intergenic (+130/+59)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base . (4/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTG  >  NC_002947/343895‑344184
                                                                                                                                                   |                                                                                                                                                
ccATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCAc                                                                                                                                                 >  5:184731/1‑148 (MQ=255)
     tCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCg                                                                                                                                          <  6:184731/149‑1 (MQ=255)
             gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGtttttt                                                                                                                                  >  3:136825/1‑149 (MQ=255)
             gTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                                  <  5:249808/149‑1 (MQ=255)
              tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttg                                                                                                                                 >  7:212655/1‑149 (MQ=255)
              tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                                 >  1:191030/1‑149 (MQ=255)
              tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttg                                                                                                                                 >  5:71907/1‑149 (MQ=255)
               aaGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTTCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttg                                                                                                                                 <  4:136825/148‑1 (MQ=255)
                aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttgtt                                                                                                                               <  8:212655/149‑1 (MQ=39)
                        cccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTc                                                                                                                       <  1:89299/149‑1 (MQ=255)
                         ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGTCGGGTTTCTTGTTGTCGCTTCg                                                                                                                      >  8:187877/1‑149 (MQ=255)
                         ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCg                                                                                                                      >  5:14505/1‑149 (MQ=255)
                                                     tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                           >  8:8669/1‑148 (MQ=255)
                                                       gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                        >  8:217455/1‑149 (MQ=255)
                                                       gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAg                                                                                        >  8:221190/1‑149 (MQ=255)
                                                           gCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtg                                                                                     <  5:261925/148‑1 (MQ=255)
                                                             aTGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGgtggtg                                                                                  >  1:36335/1‑149 (MQ=255)
                                                                    gCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCaga                                                                            <  7:217455/148‑1 (MQ=255)
                                                                    gCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCaga                                                                            <  7:221190/148‑1 (MQ=255)
                                                                     ccTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAgaga                                                                          <  6:28267/149‑1 (MQ=255)
                                                                     ccTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAgaga                                                                          >  7:62154/1‑149 (MQ=255)
                                                                       ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                        >  1:238369/1‑149 (MQ=255)
                                                                         cAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                        >  3:189449/1‑147 (MQ=255)
                                                                           ggTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACATACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGt                                                                    <  2:191030/149‑1 (MQ=255)
                                                                           ggTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGt                                                                    <  7:116722/149‑1 (MQ=255)
                                                                            gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                                <  8:205082/136‑1 (MQ=255)
                                                                            gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                                >  7:205082/1‑136 (MQ=255)
                                                                                cGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAc                                                               >  2:9041/1‑149 (MQ=255)
                                                                                     tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                           >  7:212894/1‑148 (MQ=255)
                                                                                      gAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCaga                                                                            >  3:174871/1‑130 (MQ=255)
                                                                                      gAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCaga                                                                            <  4:174871/130‑1 (MQ=255)
                                                                                          tctcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCAAGAGACGGAGTCGAACCGCCGacac                                                      >  7:215216/1‑148 (MQ=255)
                                                                                               tgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAt                                                <  8:212894/149‑1 (MQ=255)
                                                                                                  tcttcGACGCCGCCTGAAGTTAAAAATCACAGACAATCAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAttt                                              >  1:180528/1‑149 (MQ=255)
                                                                                                   cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGt                                                                                                                   <  2:218525/78‑1 (MQ=255)
                                                                                                   cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                             <  6:14505/148‑1 (MQ=255)
                                                                                                   cttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAtttt                                             <  7:8669/148‑1 (MQ=255)
                                                                                                     tcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                                >  6:180615/1‑111 (MQ=255)
                                                                                                     tcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                                <  5:180615/111‑1 (MQ=255)
                                                                                                      cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCaa                                          <  1:9041/148‑1 (MQ=255)
                                                                                                      cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCa                                           <  6:71907/147‑1 (MQ=255)
                                                                                                          gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                      <  7:44834/116‑1 (MQ=255)
                                                                                                          gccgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                      >  8:44834/1‑116 (MQ=255)
                                                                                                            cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctcgt                                   <  8:62154/149‑3 (MQ=255)
                                                                                                            cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctc                                     <  6:184467/147‑1 (MQ=255)
                                                                                                            cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTGAAAATCctc                                     >  5:184467/1‑147 (MQ=255)
                                                                                                                 gAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTa                              >  5:131050/1‑149 (MQ=255)
                                                                                                                  aaGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAc                             <  2:238369/149‑1 (MQ=255)
                                                                                                                  aaGTTAAAAATCACAGACAACAAAAAACACTCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTa                              <  7:187877/148‑1 (MQ=255)
                                                                                                                       aaaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACt                        >  1:250893/1‑149 (MQ=255)
                                                                                                                                   caacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGTTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg            <  2:250893/149‑1 (MQ=255)
                                                                                                                                   caacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTggg            <  2:180528/149‑1 (MQ=255)
                                                                                                                                    aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc           >  8:220207/1‑149 (MQ=255)
                                                                                                                                       caaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTg                       <  1:231712/135‑1 (MQ=255)
                                                                                                                                       caaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTg                       >  2:231712/1‑135 (MQ=255)
                                                                                                                                        aagaaaCCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACGAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac        >  2:209311/4‑149 (MQ=255)
                                                                                                                                        aagaaaCCCGCCTAGTGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACAACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTACAATCCTCTGATCTACCGACTGAGCTATCTGGGcgac        >  6:76130/4‑149 (MQ=255)
                                                                                                                                         aaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgacg       <  2:36335/149‑1 (MQ=255)
                                                                                                                                               ccGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCTCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTg  >  5:12312/1‑148 (MQ=255)
                                                                                                                                               ccGCACTAGGCGCGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTg  >  7:208431/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CCATGTCGGACGAGTAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGTG  >  NC_002947/343895‑344184

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: