Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,389,951 C→G 6.4% P25A (CCG→GCG)  PP_1212 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,389,9510CG 6.4% 280.5 / 5.4 89P25A (CCG→GCG) PP_1212hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (49/34);  new base G (3/3);  total (52/37)
Fisher's exact test for biased strand distribution p-value = 6.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.93e-01

CTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCC  >  NC_002947/1389807‑1390094
                                                                                                                                                 |                                                                                                                                               
cttgctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCt                                                                                                                                             <  4:31886/149‑1 (MQ=255)
   gctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGa                                                                                                                                           <  6:356035/148‑1 (MQ=255)
    ctAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATc                                                                                                                                         <  8:208486/149‑1 (MQ=255)
         cATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGGCGCGGCGCTGATcttctg                                                                                                                                    <  1:27481/149‑6 (MQ=255)
          atatatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCTAGCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTg                                                                                                                                    <  1:300359/149‑1 (MQ=255)
          atatatAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGc                                                                                                                                   >  3:151577/1‑149 (MQ=255)
                 aTCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCt                                                                                                                            >  2:261610/1‑149 (MQ=255)
                  tCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTg                                                                                                                           >  1:54548/1‑149 (MQ=255)
                   cccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  7:274516/1‑149 (MQ=255)
                    ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGCCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGAGGAGCGATTGCCCGGCATGTc                                                                                                                         >  4:273790/1‑149 (MQ=255)
                    ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  1:251892/1‑148 (MQ=255)
                    ccGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  4:181621/1‑148 (MQ=255)
                    ccGCTTTTTGCTGCCCCCAGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGAGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATAAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  4:273771/1‑148 (MQ=255)
                     cGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGt                                                                                                                          >  6:187011/1‑147 (MQ=255)
                     cGCTTTTTGCTGACCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCAAGCCCCTTTATGACTATCAATGTGCGTCTT‑GCGAGCAACGAATGGAAATGCTGCAGAAGATCAGCGCCGCGCCGCTCATCGATTGCACGGACTgacc                                                                                                                        >  8:156395/1‑146 (MQ=255)
                         ttttGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCGGATCGATTGCCCGGCCTGTCAGgcg                                                                                                                    >  5:325781/1‑149 (MQ=255)
                              ctgcCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCgg                                                                                                                <  2:251892/148‑1 (MQ=255)
                                 ccccccGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGACAGgccgggcaa                                                                                                              <  1:258910/147‑8 (MQ=255)
                                        gCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                     <  1:244703/133‑1 (MQ=255)
                                        gCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGgc                                                                                                                     >  2:244703/1‑133 (MQ=255)
                                             cagcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAgctgct                                                                                                >  1:103915/1‑149 (MQ=255)
                                              agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCgg                                                                                                                               >  7:211100/1‑117 (MQ=255)
                                              agcGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCgg                                                                                                                               <  8:211100/117‑1 (MQ=255)
                                                              agagaCACGCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAACGATGCTGTCGGTGCCAGGCTTTc                                                                               >  2:90152/1‑149 (MQ=255)
                                                              agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                               <  3:181621/149‑1 (MQ=255)
                                                              agagaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTc                                                                               <  6:118985/149‑1 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCTCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGGTGTCGGTGCCAGGCTTTCGc                                                                             <  1:243575/149‑1 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGTCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                             <  1:90152/149‑1 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAgg                                                                                    <  4:23029/142‑1 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGc                                                                             >  1:198775/1‑149 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTACTGTCGGTGCCAGGCTTTCGc                                                                             <  1:261610/149‑1 (MQ=255)
                                                                agaCACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGATGCTGTCGGTGCCAGGcagctcg                                                                             >  3:23029/1‑143 (MQ=255)
                                                                    acCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCt                                                                           >  1:234016/1‑147 (MQ=255)
                                                                      ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                              >  5:208299/1‑110 (MQ=255)
                                                                      ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCa                                                                                                              <  6:208299/110‑1 (MQ=255)
                                                                      ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGc                                                                       >  3:154934/1‑149 (MQ=255)
                                                                      ccATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGc                                                                       >  3:267433/1‑149 (MQ=255)
                                                                       cATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCg                                                                      >  5:12865/1‑149 (MQ=255)
                                                                             ccTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACg                                                                >  5:109309/1‑149 (MQ=255)
                                                                             ccTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACg                                                                >  2:251217/1‑149 (MQ=255)
                                                                              cTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACgg                                                               >  1:351562/1‑149 (MQ=255)
                                                                                 tATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                              <  8:274516/147‑1 (MQ=255)
                                                                                 tATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGt                                                              <  1:335519/147‑1 (MQ=255)
                                                                                    gACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         >  5:210196/1‑149 (MQ=255)
                                                                                    gACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         >  5:339040/1‑149 (MQ=255)
                                                                                     aCTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTa                                                         <  6:12865/148‑1 (MQ=255)
                                                                                     aCTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        >  2:13610/1‑149 (MQ=255)
                                                                                      cTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        >  2:30296/1‑148 (MQ=255)
                                                                                       tATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCAGCAGAAGATCAGCGCCGCGCTGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACg                                                       >  1:263909/1‑148 (MQ=255)
                                                                                           aaTGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAAcc                                                  >  3:375852/1‑149 (MQ=255)
                                                                                           aaTGTGCCTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACAGTTGGTACGAAAcc                                                  <  8:184952/149‑1 (MQ=255)
                                                                                               tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCTGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGAc                                               <  2:263909/148‑1 (MQ=255)
                                                                                               tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACt                                              <  6:109309/149‑1 (MQ=255)
                                                                                               tgCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACt                                              <  3:273790/149‑1 (MQ=255)
                                                                                                 cGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            <  1:251217/149‑1 (MQ=255)
                                                                                                  gTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  1:7138/1‑148 (MQ=255)
                                                                                                  gTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  2:26693/1‑148 (MQ=255)
                                                                                                   tCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCGGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTc                                            >  6:22244/1‑147 (MQ=255)
                                                                                                     cT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGa                                        >  4:58483/1‑149 (MQ=255)
                                                                                                        gCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCACCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGAcc                                      >  1:25611/1‑149 (MQ=255)
                                                                                                         cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTAc                                                        <  5:359169/130‑1 (MQ=255)
                                                                                                         cGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCAGATGTCGGTGCCAGGCTTTCGCCTGAGCGCTAACGGTTGGTAc                                                        >  6:359169/1‑130 (MQ=255)
                                                                                                                cGCATGGAAGTGCTGCAGAAGATCAGCACCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaa                              <  2:25611/149‑1 (MQ=255)
                                                                                                                cGCATGGAAGTGCTGCAGAAGAGCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTTGTACGAAACCGACTTCAAGACCGGGGCGaa                              <  3:273771/149‑1 (MQ=255)
                                                                                                                 gCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCTGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaa                              >  6:293146/1‑148 (MQ=255)
                                                                                                                  cATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaa                             >  7:306612/1‑148 (MQ=255)
                                                                                                                        aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCg                                                 <  4:321440/122‑1 (MQ=255)
                                                                                                                        aGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCg                                                 >  3:321440/1‑122 (MQ=255)
                                                                                                                          tgctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTg                    <  4:151577/149‑1 (MQ=255)
                                                                                                                            ctgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGc                  <  3:58483/149‑1 (MQ=255)
                                                                                                                             tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCTCTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATc                      >  7:66946/1‑144 (MQ=255)
                                                                                                                             tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAACCCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   >  4:95624/1‑147 (MQ=255)
                                                                                                                             tgcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATc                      <  8:66946/144‑1 (MQ=255)
                                                                                                                              gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCGGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAg                >  8:75557/1‑149 (MQ=255)
                                                                                                                              gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   <  3:149738/146‑1 (MQ=255)
                                                                                                                              gcAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTgg                   >  4:149738/1‑146 (MQ=255)
                                                                                                                                agaagaTCAGCGCCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc              <  2:54548/149‑1 (MQ=255)
                                                                                                                                agaagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggc              <  6:257802/149‑1 (MQ=255)
                                                                                                                                 atcgaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcg             >  5:215860/4‑149 (MQ=255)
                                                                                                                                  aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                           >  6:354675/1‑134 (MQ=255)
                                                                                                                                  aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGaaaaa                           <  5:354675/134‑1 (MQ=255)
                                                                                                                                  aagaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACAGGGGCGAAAAAGAATCTGGCAggcgg            >  7:349576/1‑149 (MQ=255)
                                                                                                                                    gaTCAGCGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGgcggcg          >  1:361440/1‑149 (MQ=255)
                                                                                                                                    gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAggcggc           >  2:227450/1‑148 (MQ=255)
                                                                                                                                    gaTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGgcggcg          <  4:154934/149‑1 (MQ=255)
                                                                                                                                     aTCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         <  6:325781/149‑1 (MQ=255)
                                                                                                                                     aTCAGAGGCGCGGCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         >  5:136570/1‑149 (MQ=255)
                                                                                                                                      tCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGa         >  3:16168/1‑148 (MQ=255)
                                                                                                                                        aGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGAc        >  2:264309/1‑147 (MQ=255)
                                                                                                                                         gcggcgcggcGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg     >  8:257558/1‑149 (MQ=255)
                                                                                                                                         gcgccgcgccGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAg     <  6:339040/149‑1 (MQ=255)
                                                                                                                                            ccgcgccgcTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAACCTGCTGTCGGGGCCAGGCTTGCGCCTGAGCGGTAACGGGTGGTACGAAACCGACTTCCAGCCCGGGGCGAAAAAGACTCTGGCAGGCGGCGTCAAGGcc  <  6:196173/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
CTTGCTAACCATATATAATCCCGCTTTTTGCTGCCCCCCGGCAGCCAGCGTTATGGTTGACGAGAGACACCCATGCCCCTTTATGACTATCAATGTGCGTCCT‑GCGAGCACCGCATGGAAGTGCTGCAGAAGATCAGCGCCGCGCCGCTGATCGATTGCCCGGCCTGTCAGGCGCCGGCACTGAAAAAGCTGCTGTCGGTGCCAGGCTTTCGCCTGAGCGGTAACGGTTGGTACGAAACCGACTTCAAGACCGGGGCGAAAAAGAATCTGGCAGGCGGCGACAAGGCC  >  NC_002947/1389807‑1390094

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