Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,082,339 A→G 5.9% L232S (TTG→TCG)  purK ← 5‑(carboxyamino)imidazole ribonucleotide synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,082,3390AG 5.9% 193.9 / 2.1 67L232S (TTG→TCG) purK5‑(carboxyamino)imidazole ribonucleotide synthase
Reads supporting (aligned to +/- strand):  ref base A (34/29);  new base G (2/2);  total (36/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

AGATGGTTCTCGAACTGGCTGCACTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCCACGCACGGCCACCA  >  NC_002947/6082192‑6082483
                                                                                                                                                   |                                                                                                                                                
aGATGGTTCTCGAACTGGCTGCACTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCaa                                                                                                                                                 >  2:217586/1‑149 (MQ=255)
         tcGAACTGGCTGCACTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGcc                                                                                                                                        >  7:3368/1‑149 (MQ=255)
          cGAACTGGCTGCACTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGcc                                                                                                                                        >  7:307541/1‑148 (MQ=255)
             aCTGGCTGCACTCGTCTCATTCGATGGTCCAGTGCCCGGAGTTTTGCACGCGCGGGGCTATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACg                                                                                                                                    <  8:307541/149‑1 (MQ=255)
                      aCTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCg                                                                                                                           <  3:92991/149‑1 (MQ=255)
                        tctgcgcTTTCGATGTTCTAGTGCCGGGAGTTGTGCACGCGCGGGGCTATTTCGTTGGCCTTCAAGCGACCTTAAATTTGGTAGAAATCTATGTCAATCACGCAAACGTATTCGAGCTTCTGCAACTCGCGGCCGACGTAATCTTCGGc                                                                                                                         <  4:115726/146‑1 (MQ=255)
                        tCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGc                                                                                                                         >  3:115726/1‑149 (MQ=255)
                          ggCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGc                                                                                                                         >  1:165586/1‑147 (MQ=255)
                             gcCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCagag                                                                                                                    >  5:184316/1‑149 (MQ=255)
                              ccTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGc                                                                                                                   <  8:3368/149‑1 (MQ=255)
                                  cGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCt                                                                                                                 <  1:217586/147‑1 (MQ=255)
                                             tGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGc                                                                                                    <  1:217266/149‑1 (MQ=255)
                                               cccGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCct                                                                                                  >  4:269123/1‑149 (MQ=255)
                                                    aGTTGAGCACGCGCGGGGCGAGTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGACGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCctggctg                                                                                             <  6:178393/149‑1 (MQ=255)
                                                           cACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTAcc                                                                                       >  5:97288/1‑148 (MQ=255)
                                                           cACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTAc                                                                                        >  6:86072/1‑147 (MQ=255)
                                                           cACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTAc                                                                                        <  5:86072/147‑1 (MQ=255)
                                                                      cGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGa                                                                           <  5:260331/149‑1 (MQ=255)
                                                                       gATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGaa                                                                          >  5:32728/1‑149 (MQ=255)
                                                                         tttCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATg                                                                        >  1:329288/1‑149 (MQ=255)
                                                                               tGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGctc                                                                  >  6:278686/1‑149 (MQ=255)
                                                                                   cTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTAc                                                                                        <  3:148421/123‑1 (MQ=255)
                                                                                   cTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTAc                                                                                        >  4:148421/1‑123 (MQ=255)
                                                                                      cAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            <  6:184316/148‑1 (MQ=255)
                                                                                      cAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            <  6:32728/148‑1 (MQ=255)
                                                                                       aaGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAgcgc                                                                              <  4:243532/129‑1 (MQ=255)
                                                                                       aaGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAgcgc                                                                              >  3:243532/1‑129 (MQ=255)
                                                                                       aaGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTggg                                                           >  3:340738/1‑148 (MQ=255)
                                                                                         gCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTAt                                                        >  5:289839/1‑149 (MQ=255)
                                                                                           caccaTCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGt                                                              <  7:328076/141‑1 (MQ=255)
                                                                                           caccaTCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGt                                                              >  8:328076/1‑141 (MQ=255)
                                                                                                tCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCAc                                                 <  2:329288/149‑1 (MQ=255)
                                                                                                 cAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCAc                                                 >  2:144883/1‑148 (MQ=255)
                                                                                                  aaCTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCGACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGa                                                                                     <  8:238554/111‑1 (MQ=255)
                                                                                                  aaCTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCGACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGa                                                                                     >  7:238554/1‑111 (MQ=255)
                                                                                                  aaCTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCGACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGa                                                                                     <  2:198972/111‑1 (MQ=255)
                                                                                                  aaCTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCGACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGa                                                                                     >  1:198972/1‑111 (MQ=255)
                                                                                                  aaCTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCAc                                                 >  1:35717/1‑147 (MQ=255)
                                                                                                        gaagaaCTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGt                                         <  6:289839/149‑1 (MQ=255)
                                                                                                         aagaaCTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCggg                                          >  4:245643/1‑147 (MQ=255)
                                                                                                            aaCTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGaa                                     <  1:144883/149‑1 (MQ=255)
                                                                                                              cTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGc                                   >  2:31198/1‑149 (MQ=255)
                                                                                                               tCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCg                                  >  3:170488/1‑149 (MQ=255)
                                                                                                                 gAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGt                                         >  3:3171/1‑140 (MQ=255)
                                                                                                                 gAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGt                                         <  4:3171/140‑1 (MQ=255)
                                                                                                                 gAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCggg                                >  2:37161/1‑149 (MQ=255)
                                                                                                                  aaGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGt                               >  1:301556/1‑149 (MQ=255)
                                                                                                                     gCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTc                            <  2:165586/149‑1 (MQ=255)
                                                                                                                        aTCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            <  5:86682/114‑1 (MQ=255)
                                                                                                                        aTCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            >  6:86682/1‑114 (MQ=255)
                                                                                                                           aCGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATc                      <  2:35717/149‑1 (MQ=255)
                                                                                                                             gCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGCATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATcgc                    >  7:122256/1‑149 (MQ=255)
                                                                                                                               cAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCgg                  <  5:278686/149‑1 (MQ=255)
                                                                                                                                aaCGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGc                 <  4:170488/149‑1 (MQ=255)
                                                                                                                                 aCGTAATCCAGCTTCTGCAACACTCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGAGAAGCGCAGAATGCCGCTCTCGTTGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGcc                <  3:269123/149‑1 (MQ=255)
                                                                                                                                   gTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCcac              <  3:245643/149‑1 (MQ=255)
                                                                                                                                   gTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCcac              >  2:345866/1‑149 (MQ=255)
                                                                                                                                       tCCAGATTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCcacgcac          >  7:225800/1‑149 (MQ=255)
                                                                                                                                        ccAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATcgc                    <  6:295334/138‑1 (MQ=255)
                                                                                                                                        ccAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATcgc                    >  5:295334/1‑138 (MQ=255)
                                                                                                                                        ccAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCcacgcacg         >  8:79665/1‑149 (MQ=255)
                                                                                                                                         cAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCCACGCACgg        >  6:192909/1‑149 (MQ=255)
                                                                                                                                           gCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTc                            >  6:315811/1‑127 (MQ=255)
                                                                                                                                           gCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTc                            <  5:315811/127‑1 (MQ=255)
                                                                                                                                               cTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATtg                                                      >  2:301261/1‑97 (MQ=255)
                                                                                                                                                  cAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            >  7:188228/1‑88 (MQ=255)
                                                                                                                                                  cAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTgg                                                            <  8:188228/88‑1 (MQ=255)
                                                                                                                                                  cAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCCACGCACGGccacca  <  4:340738/146‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
AGATGGTTCTCGAACTGGCTGCACTCGGCGCCTTCGATGGTCCAGTGCCCGGAGTTGTGCACGCGCGGGGCGATTTCGTTGGCCTTCAAGCCACCATCAACTTCGAAGAACTCGAAGGCCATCACGCCAACGTAATCCAGCTTCTGCAACACGCGGCCGACGTAATCTTCGGCCAGAGCCTGCAGCGGGTGTGCCTGGCTGGCTACCGACAAGCGCAGAATGCCGCTCTCGTGGGTATTGTGCACCAGCGGGTAGAAGCGGGTTTCACCATCGCGGCCACGCACGGCCACCA  >  NC_002947/6082192‑6082483

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: ‑