Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,242,983 G→A 8.3% intergenic (‑99/‑126) PP_1083 ← / → tsaA (2Fe‑2S)‑binding protein/peroxiredoxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,242,9830GA 8.3% 199.5 / 4.8 71intergenic (‑99/‑126)PP_1083/tsaA(2Fe‑2S)‑binding protein/peroxiredoxin
Reads supporting (aligned to +/- strand):  ref base G (28/37);  new base A (4/2);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 3.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

GCATCGCGGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCGTT  >  NC_002947/1242835‑1243123
                                                                                                                                                     |                                                                                                                                            
gcatcgcGGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCg                                                                                                                                              >  1:256396/1‑149 (MQ=255)
   tcgcGGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCggtg                                                                                                                                           >  7:163394/1‑149 (MQ=255)
   tcgcGGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGGTGGGGCGAGGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCggtg                                                                                                                                           >  8:95175/1‑149 (MQ=255)
     gcgGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCggtggt                                                                                                                                         <  8:79635/149‑1 (MQ=255)
      cgGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTc                                                                                                                                        <  5:196295/149‑1 (MQ=255)
              tCCGTCGGTAACGCCGACACAAAGACACTCATACATAAGGGCTGACGATCGCTGGTTGGGTCGATGGTTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATGTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCttt                                                                                                                                   <  7:48806/146‑1 (MQ=255)
              tCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGGTGGGGCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCttt                                                                                                                                   <  7:95175/146‑1 (MQ=255)
                     gTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTACTGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAg                                                                                                                         >  5:285925/1‑149 (MQ=255)
                     gTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGc                                                                                                                                      <  1:307739/136‑1 (MQ=255)
                     gTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGc                                                                                                                                      >  2:307739/1‑136 (MQ=255)
                        aCGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCg                                                                                                                      <  7:286633/149‑1 (MQ=255)
                                 caAAGACACACATACATCAGGGCTGACCATCGCTGGTTGGGGCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCAt                                                                                                             <  8:333582/149‑1 (MQ=255)
                                 caAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTctt                                                                                                             <  5:351564/149‑3 (MQ=255)
                                  aaaGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCAt                                                                                                             >  4:178694/1‑148 (MQ=255)
                                      acacacATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGaa                                                                                                        <  2:256396/149‑1 (MQ=255)
                                        acacaTACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGc                                                                                                      >  3:380072/1‑149 (MQ=255)
                                        acacaTACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCAGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGc                                                                                                      <  8:50571/149‑1 (MQ=255)
                                         cacaTACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGcc                                                                                                     <  2:127544/149‑1 (MQ=255)
                                          acatacatAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGc                                                                                                      >  2:370677/1‑147 (MQ=255)
                                          acatacatAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGc                                                                                                      <  1:370677/147‑1 (MQ=255)
                                          acatacatAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCg                                                                                                    >  1:352880/1‑149 (MQ=255)
                                           catacataAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCgt                                                                                                   >  3:104496/1‑149 (MQ=255)
                                           catacataAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCgt                                                                                                   <  4:104496/149‑1 (MQ=255)
                                             tacataAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACAGAGTCGGCCTGTCATGAGAAGCCgtgt                                                                                                 <  6:357270/149‑1 (MQ=255)
                                              acataAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTa                                                                                                >  5:153068/1‑149 (MQ=255)
                                                     ggCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGt                                                                                         <  6:379495/149‑1 (MQ=255)
                                                     ggCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGt                                                                                         >  1:226315/1‑149 (MQ=255)
                                                       cTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAGTGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCg                                                                                       <  2:10693/149‑1 (MQ=255)
                                                          aCCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTATGTCGACGGAGTCGGCCTGTGATGAGAAGCAGTGTATGATGGTCGGc                                                                                     >  8:284326/1‑148 (MQ=255)
                                                           ccATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCt                                                                                   <  3:178694/149‑1 (MQ=255)
                                                            cATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCtt                                                                                  >  8:167979/1‑149 (MQ=255)
                                                             aTCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCttt                                                                                 <  6:392501/149‑1 (MQ=255)
                                                              tCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTttg                                                                                >  6:311563/1‑149 (MQ=255)
                                                                 cTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTttgttg                                                                             <  4:380072/149‑1 (MQ=255)
                                                                 cTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTttgttg                                                                             <  8:363789/149‑1 (MQ=255)
                                                                 cTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTttgttg                                                                             <  6:162816/149‑1 (MQ=255)
                                                                    gTTTGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTGATGAGAAGCCGTGTATGGTGGTCGGCCTTTGTTgg                                                                            <  5:11113/147‑1 (MQ=255)
                                                                    gTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTc                                                                                                                                >  4:267213/1‑95 (MQ=255)
                                                                    gTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTc                                                                                                                                <  3:267213/95‑1 (MQ=255)
                                                                    gTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGAt                                                                          <  2:352880/149‑1 (MQ=255)
                                                                    gTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGAAGGAGTCGGCCTGTCATGATAAGCCGTGTatgat                                                                                            >  5:404809/1‑131 (MQ=255)
                                                                         gTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCgg                                                                     <  3:126317/149‑1 (MQ=255)
                                                                         gTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCgg                                                                     <  2:226315/149‑1 (MQ=255)
                                                                         gTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGATTGAGCTAACATTTGTCTATTTGTCGAAGTAGTAGGCCTGTCATGAGAAGCGGTGTATGATGGTCGGACTTTGTTGGATgaagg                                                                     >  7:151887/1‑145 (MQ=255)
                                                                          tCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATg                                                                                                            <  4:377508/109‑1 (MQ=255)
                                                                                  ttGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGccc                                                             >  8:139225/1‑148 (MQ=255)
                                                                                     gAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCt                                                         >  1:81906/1‑149 (MQ=255)
                                                                                       gCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTgg                                                       <  1:201048/149‑1 (MQ=255)
                                                                                            atgttaatgAGAATGCTTGTCAAAGATCTTTATGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTCTTGGATGTCGGGCAAGCCCGTCTGGACCgg                                                  >  6:59538/1‑149 (MQ=255)
                                                                                              gttaatgAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGc                                                 >  8:336171/1‑148 (MQ=255)
                                                                                              gttaatgAGAATGCTTGTCAAAGATCTTTCAGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGc                                                 >  7:261701/1‑148 (MQ=255)
                                                                                               ttaatgAGAATGCTTGTCAAAGATCTTTCTGAGGGTGGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGc                                                 >  1:277096/1‑148 (MQ=255)
                                                                                                  atgAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCttt                                                                                 >  8:359663/1‑112 (MQ=255)
                                                                                                  atgAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCttt                                                                                 <  7:359663/112‑1 (MQ=255)
                                                                                                           cTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCa                                    <  2:81906/148‑1 (MQ=255)
                                                                                                           cTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCAc                                   <  7:167979/149‑1 (MQ=255)
                                                                                                            ttGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCAcc                                  >  8:371408/1‑149 (MQ=255)
                                                                                                            ttGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCAc                                   >  5:77969/1‑148 (MQ=255)
                                                                                                               tCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCtc                               <  7:371408/149‑1 (MQ=255)
                                                                                                                 aaaGATCTTTCTGAGGGTGCTGTGCTTGAGCGACGGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCtca                              <  8:151887/148‑1 (MQ=255)
                                                                                                                    gATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAg                          <  7:336171/149‑1 (MQ=255)
                                                                                                                                ggTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGga                        >  4:7358/1‑139 (MQ=255)
                                                                                                                                ggTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGga                        <  3:7358/139‑1 (MQ=255)
                                                                                                                                ggTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAAt               <  2:277096/148‑1 (MQ=255)
                                                                                                                                 gtgctgtgctTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGa             <  7:300129/149‑1 (MQ=255)
                                                                                                                                 gtgctgtgctTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGa             <  6:77969/149‑1 (MQ=255)
                                                                                                                                    gctgtgctTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGc           >  4:244020/1‑149 (MQ=255)
                                                                                                                                    gctgtgctTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTACACCCTCATCAGGAGATTACAATGAGc           >  8:353509/1‑149 (MQ=255)
                                                                                                                                           ttGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCg    >  1:8542/1‑149 (MQ=255)
                                                                                                                                           ttGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCg    >  2:187301/1‑149 (MQ=255)
                                                                                                                                           caaaGCGACCACCGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCg    >  7:207236/4‑149 (MQ=255)
                                                                                                                                            aaaGCGACCACCGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGc           >  5:277822/3‑141 (MQ=255)
                                                                                                                                            aaaGCGACCACCGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGc           <  6:277822/139‑1 (MQ=255)
                                                                                                                                             aaGCGACCACCGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCGtt  >  2:141779/2‑149 (MQ=255)
                                                                                                                                              aGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCGtt  <  7:284326/148‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
GCATCGCGGATCTGTCCGTCGGTAACGCCGACACAAAGACACACATACATAAGGGCTGACCATCGCTGGTTGGGTCGATGGGTTGGAGCTTAATGTTAATGAGAATGCTTGTCAAAGATCTTTCTGAGGGTGCTGTGCTTGAGCGACCGGTGGTCGCTTTGTCGACGGAGTCGGCCTGTCATGAGAAGCCGTGTATGATGGTCGGCCTTTGTTGGATGTCGGGCAAGCCCGTCTGGCCCGGCAAACGGTTTTTCACCCTCATCAGGAGATTACAATGAGCGTACTCGTT  >  NC_002947/1242835‑1243123

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: