Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 344,041 Δ1 bp 7.7% 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. 7.7% 201.5 / 10.7 65intergenic (+130/+59)PP_5440/PP_t72hypothetical protein/tRNA‑Phe
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base . (2/3);  total (30/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

TAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGT  >  NC_002947/343909‑344183
                                                                                                                                     |                                                                                                                                               
tAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCtt                                                                                                                                 <  4:80332/148‑1 (MQ=255)
 aaGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCCTAGTGCGGGTTTTttg                                                                                                                                <  5:155355/148‑1 (MQ=255)
  aGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgt                                                                                                                               <  8:289959/148‑1 (MQ=255)
   ggTTTGGCCTGGTGTCTGATGATGTTTGCGGGTTTCTCTATTGCATGGGTTGCCTTCAGGTTGCGTTGTGAGGTGTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATAACAGACAACAAAATACCAGCTCTAGGGGGGTTTCttgtt                                                                                                                              <  4:14604/148‑1 (MQ=255)
    gTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGTTTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTtgttgttgt                                                                                                                            >  8:199263/1‑149 (MQ=255)
         gCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTTTGTTGTctgt                                                                                                                        <  5:288608/148‑4 (MQ=38)
         gCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCtt                                                                                                                       <  4:149834/149‑1 (MQ=255)
         gCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCtt                                                                                                                       >  5:144343/1‑149 (MQ=255)
           ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttgtt                                                                                                                              >  6:146431/1‑140 (MQ=39)
           ccGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAGAAACCCGCCTAGTGCGGGTTTTttgtt                                                                                                                              <  5:146431/140‑1 (MQ=39)
              cAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGt                                                                                                                  <  7:259815/149‑1 (MQ=255)
                        aTGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                          <  6:187093/147‑1 (MQ=255)
                         tGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAac                                                                                                       <  1:287548/149‑1 (MQ=255)
                            cTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                    <  2:287035/149‑1 (MQ=255)
                             tGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACttt                                                                                                            <  5:43000/140‑1 (MQ=255)
                             tGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACttt                                                                                                            >  6:43000/1‑140 (MQ=255)
                                      ctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                    >  6:278936/1‑139 (MQ=255)
                                      ctcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCtt                                                                                           >  3:279662/1‑148 (MQ=255)
                                      ctcGATTGCATGGCTAGCCTTCAGGTTGCGTAGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAaccac                                                                                                    <  5:278936/139‑1 (MQ=255)
                                       tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                          >  6:13211/1‑148 (MQ=255)
                                       tcGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTg                                                                                          >  3:255682/1‑148 (MQ=255)
                                         gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                          >  8:930/1‑130 (MQ=255)
                                         gATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGc                                                                                                          <  7:930/130‑1 (MQ=255)
                                           ttGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggc                                                                                     <  5:13211/149‑2 (MQ=255)
                                              cATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAggtggt                                                                                  >  7:228736/1‑149 (MQ=255)
                                                tGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGc                                                                                <  4:279662/149‑1 (MQ=255)
                                                   cTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCa                                                                             <  8:228736/149‑1 (MQ=255)
                                                    tAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGcc                                                                               >  2:40388/1‑146 (MQ=255)
                                                         ttCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACg                                                                       >  7:130986/1‑149 (MQ=255)
                                                            aGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGa                                                                     <  3:77981/148‑1 (MQ=255)
                                                              gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTc                                                                  >  2:153020/1‑149 (MQ=255)
                                                              gTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTc                                                                  <  4:255682/149‑1 (MQ=255)
                                                               ttGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTc                                                                  >  8:184443/1‑148 (MQ=255)
                                                                       tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgccg                                                         >  2:290361/1‑149 (MQ=255)
                                                                       tgAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAAccgcc                                                          <  3:104020/148‑1 (MQ=255)
                                                                         aGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGa                                                        >  8:47949/1‑148 (MQ=255)
                                                                           gTCTCGTGTTCTTCGACGCCGCCTGACGTTAAAACTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCttgtt                                                                                                                              <  3:808/76‑1 (MQ=255)
                                                                           gTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACCGACAACAAAAAACCCGCACTAGGCGGGTGTCttgtg                                                                                                                              >  4:808/1‑75 (MQ=255)
                                                                              tcGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAg                                                  >  3:217463/1‑149 (MQ=255)
                                                                                 tgtTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGAt                                               >  6:174382/1‑149 (MQ=255)
                                                                                      ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGAAACCACCAGAAGATTGAAGTTGGTGCCCATAGACGGAGTCGAACCGCCGacac                                                     >  1:228257/1‑139 (MQ=255)
                                                                                      ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGacac                                                     <  2:228257/138‑1 (MQ=255)
                                                                                      ttcGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTc                                           >  7:85220/1‑148 (MQ=255)
                                                                                        cGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCaa                                         <  1:153020/148‑1 (MQ=255)
                                                                                           cgccgcCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCct                                     <  7:199263/149‑1 (MQ=255)
                                                                                              cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTg                                  <  6:144343/149‑1 (MQ=255)
                                                                                              cgccTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTg                                  <  5:174382/149‑1 (MQ=255)
                                                                                                      gTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTAcc                           >  1:291030/1‑148 (MQ=255)
                                                                                                       ttAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGa                         >  3:102493/1‑149 (MQ=255)
                                                                                                         aaaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGAc                        <  7:184443/148‑1 (MQ=255)
                                                                                                          aaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCaa                                         >  3:136298/1‑130 (MQ=255)
                                                                                                          aaaaTCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCaa                                         <  4:136298/130‑1 (MQ=255)
                                                                                                              tCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGc                   >  8:92611/1‑148 (MQ=255)
                                                                                                                  aGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCt              >  4:46200/1‑149 (MQ=255)
                                                                                                                    acaacaAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTCCTCTACCGACTGAGCTATCTgg            <  7:47949/149‑1 (MQ=255)
                                                                                                                      aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc          <  2:175706/149‑1 (MQ=255)
                                                                                                                      aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc          <  8:130986/149‑1 (MQ=255)
                                                                                                                      aacaaAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGc          <  4:102493/149‑1 (MQ=255)
                                                                                                                        caaGAAACCCGCCTAGTGCGGGTTTTTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcga        >  7:301912/1‑149 (MQ=255)
                                                                                                                         aaaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTAt                <  5:159561/140‑1 (MQ=255)
                                                                                                                         aaaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTAt                >  6:159561/1‑140 (MQ=255)
                                                                                                                         aaaaaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGcgac       >  6:129574/1‑149 (MQ=255)
                                                                                                                            aaaCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGTGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAg    >  4:293020/1‑149 (MQ=255)
                                                                                                                              aCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGt  <  5:129574/149‑1 (MQ=255)
                                                                                                                              aCCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGt  <  8:85220/149‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                               
TAAGGTTTGGCCCGCAGTCTGCTCATGGCTGCGGGTTTCTCGATTGCATGGCTAGCCTTCAGGTTGCGTTGTGAGGTCTCGTGTTCTTCGACGCCGCCTGAAGTTAAAAATCACAGACAACAAAAAACCCGCACTAGGCGGGTTTCTTGTTGTCGCTTCGGGTAACTTTGCAACCACCAGAAGCTTGAAGGTGGTGCCCAGAGACGGAGTCGAACCGCCGACACGAGGATTTTCAATCCTCTGCTCTACCGACTGAGCTATCTGGGCGACGAGGT  >  NC_002947/343909‑344183

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: