Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,558,481 Δ1 bp 5.3% coding (420/732 nt) PP_1369 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,558,4810C. 5.3% 246.5 / 6.4 75coding (420/732 nt)PP_1369hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (38/33);  new base . (2/2);  total (40/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 = 1.00e+00

AGTCGCACAGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCTTCATGG  >  NC_002947/1558342‑1558628
                                                                                                                                           |                                                                                                                                                    
aGTCGCACAGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGc                                                                                                                                            >  4:20324/1‑149 (MQ=255)
 gTCGCACAGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCACGCGGCAAGCGCtgtgcc                                                                                                                                           >  2:343076/1‑144 (MQ=255)
  tCGCACAGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCTGGCGCTCCAGACCAGGCCTGGCGCATCAATCGACGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCcg                                                                                                                                          >  6:295334/1‑149 (MQ=255)
      acaGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTg                                                                                                                                      >  1:209036/1‑149 (MQ=255)
       caGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCACGCGGCAAGCGCTGTGCCGCCTGc                                                                                                                                     <  1:343076/149‑1 (MQ=255)
                 gACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCAt                                                                                                                           >  1:160696/1‑149 (MQ=255)
                      tttGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGc                                                                                                                       >  3:77010/1‑148 (MQ=255)
                      tttGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGc                                                                                                                       >  5:340590/1‑148 (MQ=255)
                           gACGCCATGCGCTCCGGGCGCTCCTGACCAGGGCTGGCGCATCAATCGCCGTAAAGTCGGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGccc                                                                                                                                   <  3:318238/131‑1 (MQ=255)
                            aCGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCg                                                                                                                >  2:253491/1‑149 (MQ=255)
                                cATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGc                                                                                                            >  5:77246/1‑149 (MQ=255)
                                  tGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGa                                                                                                          >  1:47652/1‑149 (MQ=255)
                                     gcTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGc                                                                                                                                            <  2:293608/112‑1 (MQ=255)
                                     gcTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGc                                                                                                                                            >  1:293608/1‑112 (MQ=255)
                                       tCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAAc                                                                                                     >  5:313344/1‑149 (MQ=255)
                                           ggCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTg                                                                                                 >  5:312202/1‑149 (MQ=255)
                                                tCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCca                                                                                              >  8:288646/1‑147 (MQ=255)
                                                tCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCca                                                                                              >  1:62580/1‑147 (MQ=255)
                                                tCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCca                                                                                              >  5:216036/1‑147 (MQ=255)
                                                  cAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAgg                                                                                          >  2:138810/1‑149 (MQ=255)
                                                  cAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAgg                                                                                          <  2:209036/149‑1 (MQ=255)
                                                           ccTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTg                                                                                 <  3:299215/149‑1 (MQ=255)
                                                            cTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTgg                                                                                >  7:116794/1‑149 (MQ=255)
                                                               gcgcATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTc                                                                              <  6:312202/148‑1 (MQ=255)
                                                                 gcATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTc                                                                           <  6:340590/149‑1 (MQ=255)
                                                                  cATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTg                                                                                 >  2:201816/1‑142 (MQ=255)
                                                                  cATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTg                                                                                 <  1:201816/142‑1 (MQ=255)
                                                                  cATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGACGTGTACTGTGGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCg                                                                          <  5:295334/149‑1 (MQ=255)
                                                                      aaTCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTcc                                                                             <  2:17844/142‑1 (MQ=255)
                                                                      aaTCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTcc                                                                             >  1:17844/1‑142 (MQ=255)
                                                                          gCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTg                                                                                 <  7:183417/134‑1 (MQ=255)
                                                                          gCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTg                                                                                 >  8:183417/1‑134 (MQ=255)
                                                                             gTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAg                                                               <  3:329550/149‑1 (MQ=255)
                                                                              tAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCCCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAgg                                                              >  5:166458/1‑149 (MQ=255)
                                                                                  gTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcct                                                          >  8:152343/1‑149 (MQ=255)
                                                                                  gTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcct                                                          >  8:282179/1‑149 (MQ=255)
                                                                                   tCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcctc                                                         >  4:87689/1‑149 (MQ=255)
                                                                                   tCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcctc                                                         >  6:163265/1‑149 (MQ=255)
                                                                                   tCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcctc                                                         <  5:200699/149‑1 (MQ=255)
                                                                                       gCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGt                                                     <  1:138810/149‑1 (MQ=255)
                                                                                         aaGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAg                                                   >  8:304124/1‑149 (MQ=255)
                                                                                             ggTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTTGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCgg                                               >  5:146302/1‑149 (MQ=255)
                                                                                              gTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCt                                                                                   <  2:200658/112‑1 (MQ=255)
                                                                                              gTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCt                                                                                   >  1:200658/1‑112 (MQ=255)
                                                                                               tACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGt                                             <  1:253491/149‑1 (MQ=255)
                                                                                                 cAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTaa                                           <  5:163265/149‑1 (MQ=255)
                                                                                                  aGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGa                                                                                                                         >  7:8491/1‑70 (MQ=255)
                                                                                                  aGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGa                                                                                                                         <  8:8491/70‑1 (MQ=255)
                                                                                                  aGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAg                                          >  5:231492/1‑149 (MQ=255)
                                                                                                    aTGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcctc                                                         <  5:216008/132‑1 (MQ=255)
                                                                                                    aTGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCcctc                                                         >  6:216008/1‑132 (MQ=255)
                                                                                                     tGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAgcgc                                       >  2:149360/1‑149 (MQ=255)
                                                                                                       cTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTTGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCtt                                     <  6:146302/149‑1 (MQ=255)
                                                                                                        tGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGACAGTCCATAATGTAGCAGGGTAAgatattc                                    >  8:169568/1‑143 (MQ=255)
                                                                                                          aGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACATCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCt                                                                                   <  3:88625/100‑1 (MQ=255)
                                                                                                          aGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCt                                                                                   >  4:88625/1‑100 (MQ=255)
                                                                                                           ggTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGtt                                 <  4:309435/149‑1 (MQ=255)
                                                                                                           ggTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGtt                                 <  6:216036/149‑1 (MQ=255)
                                                                                                            gtgtACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTtc                                <  2:160696/149‑1 (MQ=255)
                                                                                                               tACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTtctc                              >  7:106846/1‑148 (MQ=255)
                                                                                                                aCTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCtca                             <  8:106846/148‑1 (MQ=255)
                                                                                                                    gtgtATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAATCCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCaa                         <  2:47652/148‑1 (MQ=255)
                                                                                                                     tgtATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAgg                       <  7:152343/149‑1 (MQ=255)
                                                                                                                         tttCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGt                   <  3:20324/149‑1 (MQ=255)
                                                                                                                            cAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGCACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAg                 >  8:287526/1‑148 (MQ=255)
                                                                                                                             aGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGaa               <  1:149360/149‑1 (MQ=255)
                                                                                                                              gTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAAc              <  2:62580/149‑1 (MQ=255)
                                                                                                                                  cAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCCGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGacggcgc          >  5:178674/1‑143 (MQ=255)
                                                                                                                                  cAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGc          <  3:34520/149‑1 (MQ=255)
                                                                                                                                  cAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGc          <  4:77010/149‑1 (MQ=255)
                                                                                                                                   acgcggcAAGCGCTGTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAgg                       >  2:180868/3‑135 (MQ=255)
                                                                                                                                   acgcggcAAGCGCTGTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAgg                       <  1:180868/133‑1 (MQ=255)
                                                                                                                                    ggcggcACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCt        >  6:194546/1‑149 (MQ=255)
                                                                                                                                        gcACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCTTCAt    <  7:304124/149‑1 (MQ=255)
                                                                                                                                        gcACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCTTCAt    <  7:169568/149‑1 (MQ=255)
                                                                                                                                          acaGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCTTCATgg  >  8:65858/1‑149 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
AGTCGCACAGCGTGCGGGACAGTTTGCGACGCCATGCGCTCCGGGCGCTCCAGACCAGGCCTGGCGCATCAATCGCCGTAAAGTCTGGCAAGCGGTACAGATGCTGAGGTGTACTGGTGTATTTCAGTGGCAGGCGGCACAGCGCTTGCCGCGTGCCCTGGCTCATGAGCCAGGCCGGTGCGACAAACCCTGCCACAGGCCACCCTTGGTCCTCGAACAGTGCCAGGCCCTCCTGTAGCCGGGTAAGCGCTTCGTTCTCATCAAGGCTGTAGAACTCGCCTTCATGG  >  NC_002947/1558342‑1558628

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: