Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 557,332 Δ1 bp 10.2% intergenic (+75/‑134) rpsN → / → rpsH 30S ribosomal protein S14/30S ribosomal protein S8

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947557,3320T.10.2% 259.7 / 25.0 88intergenic (+75/‑134)rpsN/rpsH30S ribosomal protein S14/30S ribosomal protein S8
Reads supporting (aligned to +/- strand):  ref base T (42/37);  new base . (5/4);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATG  >  NC_002947/557195‑557474
                                                                                                                                            |                                                                                                                                               
aaCATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCtt                                                                                                                                     >  3:298850/1‑149 (MQ=255)
       tgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgt                                                                                                                              >  8:38386/1‑149 (MQ=255)
         cGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTgtttgtt                                                                                                                              <  1:181694/147‑4 (MQ=255)
             aaGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTg                                                                                                                         <  5:358411/148‑1 (MQ=255)
               gCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCt                                                                                                                      >  5:41041/1‑149 (MQ=255)
                 tGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtg                                                                                                                    >  3:300684/1‑149 (MQ=255)
                  gCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                   >  1:385261/1‑149 (MQ=255)
                   cAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgt                                                                                                                              >  4:48699/1‑137 (MQ=255)
                   cAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgt                                                                                                                              <  3:48699/137‑1 (MQ=255)
                       gcgcGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgtttc                                                                                                                           >  2:49192/1‑136 (MQ=255)
                             tGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGa                                                                                                        >  1:340086/1‑149 (MQ=255)
                               aCGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       <  2:385261/148‑1 (MQ=255)
                               aCGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       <  2:340086/148‑1 (MQ=255)
                                    ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCt                                                                                                  <  4:300684/148‑1 (MQ=255)
                                              tCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTgttt                                                                                                                                 >  7:204659/1‑107 (MQ=39)
                                              tCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTgttt                                                                                                                                 <  8:204659/107‑1 (MQ=39)
                                              tCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGATGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       <  4:394070/133‑1 (MQ=255)
                                              tCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGATGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       >  3:394070/1‑133 (MQ=255)
                                                aaGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATgg                                                                                                                                          >  6:409912/1‑96 (MQ=255)
                                                aaGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATgg                                                                                                                                          <  5:409912/96‑1 (MQ=255)
                                                 aGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTggg                                                                                    >  8:179994/1‑149 (MQ=255)
                                                  ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                   >  2:26667/1‑117 (MQ=255)
                                                  ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                   <  1:26667/117‑1 (MQ=255)
                                                   gCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGtt                                                                                  <  2:290786/149‑1 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACCAAACAAGCCCCTATATGGGGCTTgtttc                                                                                                                                <  1:30709/102‑1 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttc                                                                                                                                >  2:30709/1‑101 (MQ=255)
                                                      tcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGTAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                               >  1:388904/1‑149 (MQ=255)
                                                        ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGTAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc                                                                             <  2:388904/149‑1 (MQ=255)
                                                        ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCTGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc                                                                             <  7:177816/149‑1 (MQ=255)
                                                        ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGc                                                                             >  1:305939/1‑149 (MQ=255)
                                                          ggTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGcc                                                                            >  6:128995/1‑148 (MQ=255)
                                                           gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                   >  4:312199/1‑108 (MQ=255)
                                                           gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTtgt                                                                                                                   <  3:312199/108‑1 (MQ=255)
                                                            tAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAACTCATTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGccc                                                                           <  1:135934/149‑1 (MQ=255)
                                                            tAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCa                                                                          >  3:328942/1‑148 (MQ=255)
                                                            tAATCGCTGCAGGCAAGGCGCGGGTGGAACGGGACCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTAGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCa                                                                          >  4:298333/1‑148 (MQ=255)
                                                              aTCGCTGCAGGCAAGGCGCGGGTGGAAGGGGATCAATCTTACGAACGCCGGTTACCTTGACAACAAACAAGCCCCTATATGGGTCTTGTTTCGTTTCTGACTTGTGACTAGAATGACCGGCTCAACTGAGCCAGGGTATTTTCCCAAAc                                                                       >  1:319661/1‑149 (MQ=255)
                                                              aTCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCaaa                                                                        <  4:328942/148‑1 (MQ=255)
                                                                        gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCgg                                                                                                    <  4:371432/110‑1 (MQ=255)
                                                                        gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCgg                                                                                                    >  3:371432/1‑110 (MQ=255)
                                                                        gCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtg                                                             <  4:298850/149‑1 (MQ=255)
                                                                         cAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACAGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtgg                                                            >  3:305708/1‑149 (MQ=255)
                                                                          aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtggt                                                           >  2:372054/1‑149 (MQ=255)
                                                                          aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGggtggt                                                           >  2:61001/1‑149 (MQ=255)
                                                                           aGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt                                                          <  7:38386/149‑1 (MQ=255)
                                                                              cgcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACAGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTc                                                       <  4:305708/149‑1 (MQ=255)
                                                                               gcgGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCa                                                      >  3:67203/1‑149 (MQ=255)
                                                                                  ggTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATg                                                                                                         >  3:132898/1‑95 (MQ=255)
                                                                                  ggTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATg                                                                                                         <  4:132898/95‑1 (MQ=255)
                                                                                      gAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCg                                                >  2:22446/1‑148 (MQ=255)
                                                                                       aaGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGGTAACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGc                                               >  1:336783/1‑148 (MQ=255)
                                                                                        aggttaGCAATCTTCCGACCTCAGGTGACCTTGCCAACAAACAAGCCACTTTATGGTGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATTACGGTCTAACCTGATCCTGTGTTTTTTCCCAAACGGTTAGGGTGGTTGTCATAGCAGCaa                                             <  2:319661/144‑1 (MQ=255)
                                                                                        aGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCaa                                             <  2:305939/149‑1 (MQ=255)
                                                                                            gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGc                                          >  8:132079/1‑148 (MQ=255)
                                                                                              gCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                       <  5:128995/149‑1 (MQ=255)
                                                                                               cAATCTTCCGACCGCCGGTGACCTTGCCGACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                       >  4:352909/1‑148 (MQ=255)
                                                                                               cAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaat                                      <  1:22446/149‑1 (MQ=255)
                                                                                                aaTCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaat                                      >  4:300466/1‑148 (MQ=255)
                                                                                                 aTCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaataa                                    <  4:67203/149‑1 (MQ=255)
                                                                                                      ttCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTc                                 >  7:209396/1‑149 (MQ=255)
                                                                                                        ccaaCCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                               <  1:410591/146‑1 (MQ=255)
                                                                                                        ccGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                <  1:372054/148‑1 (MQ=255)
                                                                                                             cgccgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTAt                           >  5:42270/1‑148 (MQ=255)
                                                                                                              gccgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       <  5:344359/71‑1 (MQ=255)
                                                                                                              gccgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                       >  6:344359/1‑71 (MQ=255)
                                                                                                               ccgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt                             <  6:335481/144‑1 (MQ=255)
                                                                                                               ccgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt                             >  5:335481/1‑144 (MQ=255)
                                                                                                                cgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAgg                       <  2:336783/149‑1 (MQ=255)
                                                                                                                     aCCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCAt                  <  3:240911/149‑1 (MQ=255)
                                                                                                                      ccTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGccc                                                                           <  2:247746/91‑1 (MQ=255)
                                                                                                                       cTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                 >  2:391168/1‑148 (MQ=255)
                                                                                                                           ccAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                               <  5:351015/130‑1 (MQ=255)
                                                                                                                           ccAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                               >  6:351015/1‑130 (MQ=255)
                                                                                                                            caacaaACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGcc            >  7:265317/1‑148 (MQ=255)
                                                                                                                             aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCa          >  5:13342/1‑149 (MQ=255)
                                                                                                                             aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCa          >  4:57349/1‑149 (MQ=255)
                                                                                                                              acgaaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCa          >  3:147575/4‑149 (MQ=255)
                                                                                                                              acaaacaaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGc             <  1:87523/145‑1 (MQ=255)
                                                                                                                              acaaacaaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGc             >  2:87523/1‑145 (MQ=255)
                                                                                                                              acaaacaaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt         <  7:179994/149‑1 (MQ=255)
                                                                                                                                 aaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt         >  4:358190/1‑147 (MQ=255)
                                                                                                                                  aacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                       <  6:291461/116‑1 (MQ=255)
                                                                                                                                  aacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                       >  5:291461/1‑116 (MQ=255)
                                                                                                                                  aacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt         >  2:248467/1‑146 (MQ=255)
                                                                                                                                  aacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa       <  1:248467/148‑1 (MQ=255)
                                                                                                                                     aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt                             <  4:116411/123‑1 (MQ=255)
                                                                                                                                     aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGt                             >  3:116411/1‑123 (MQ=255)
                                                                                                                                      aGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATg  <  6:42270/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
AACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATG  >  NC_002947/557195‑557474

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: