Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 985,937 G→C 8.7% L233F (TTG→TTC pilF → type IV pili biogenesis protein PilF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947985,9370GC 8.7% 185.9 / 4.6 65L233F (TTG→TTCpilFtype IV pili biogenesis protein PilF
Reads supporting (aligned to +/- strand):  ref base G (28/30);  new base C (3/3);  total (31/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985790‑986083
                                                                                                                                                            |                                                                                                                                                  
tGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑Tgg                                                                                                                                                   >  3:237401/1‑149 (MQ=255)
tGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑Tgg                                                                                                                                                   <  8:297183/149‑1 (MQ=255)
  aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑Tgg                                                                                                                                                   <  4:274564/147‑1 (MQ=255)
  aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑Tgg                                                                                                                                                   <  8:189455/147‑1 (MQ=255)
  aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑Tgg                                                                                                                                                   <  1:66032/147‑1 (MQ=255)
  aGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGc                                                                                                                                                 <  8:216625/149‑1 (MQ=255)
      gTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAg                                                                                                                                              <  2:213560/148‑1 (MQ=255)
       tCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAgt                                                                                                                                             <  3:67665/148‑2 (MQ=255)
        ccTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAgtgt                                                                                                                                           <  7:106375/149‑4 (MQ=255)
         cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAAt                                                                                                                                          >  2:126379/1‑149 (MQ=255)
         cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAgtgt                                                                                                                                           >  8:100820/1‑145 (MQ=255)
         cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAgtgt                                                                                                                                           >  8:106375/1‑145 (MQ=255)
         cTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAgtgt                                                                                                                                           <  7:100820/148‑4 (MQ=255)
          tACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAAt                                                                                                                                          >  6:329184/1‑148 (MQ=255)
            cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGCTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAAtt                                                                                                                                         <  1:270705/148‑1 (MQ=255)
            cGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTAc                                                                                                                                       <  2:234854/149‑1 (MQ=255)
             gAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTAc                                                                                                                                       >  1:129200/1‑148 (MQ=255)
               aaaCATTCATTATGTGCCGGCCCTGGACTACTACGTTCGTTTCATCCATTTTAGCGAACACGACGCCCTTAGCCTGCTGCTGGGCAGCCGCCTTGCTATGGTGTTCGACGATCATGTCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTAc                                                                                                                                       <  7:242947/146‑1 (MQ=255)
                 aCAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGa                                                                                                                                  >  7:281399/1‑149 (MQ=255)
                                    ccGGGACTACTAGGATCGTTTCAGACAGTTGAGCGACCACGACGCCAGTAGACTGCTGCTGGGCAGCCGCCTTGCCAAGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGtcatcac                                                                                                                >  8:214509/1‑142 (MQ=255)
                                          ctactaCGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCCTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATc                                                                                                         <  5:329184/149‑1 (MQ=255)
                                             ctacGATCGTTTCATCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTTGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCa                                                                                                        <  7:325193/147‑1 (MQ=255)
                                             ctacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAg                                                                                                       <  4:237401/148‑1 (MQ=255)
                                              tacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGc                                                                                                      >  1:384397/1‑148 (MQ=255)
                                              tacGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGc                                                                                                      >  1:370004/1‑148 (MQ=255)
                                                 gATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAAt                                                                                                   >  2:313634/1‑148 (MQ=255)
                                                  aTCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAGCGACt                                                                                                                                <  6:220621/118‑1 (MQ=255)
                                                      tttCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGt                                                                                             >  3:50203/1‑149 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAg                                                                                        >  5:353725/1‑149 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAg                                                                                        <  2:296966/149‑1 (MQ=255)
                                                           gCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAg                                                                                        >  8:105288/1‑149 (MQ=255)
                                                               gTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAAc                                                                                    >  4:87489/1‑149 (MQ=255)
                                                                   aGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGAt                                                                                 <  5:155452/148‑1 (MQ=255)
                                                                       aCCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAAcg                                                                            >  6:271559/1‑149 (MQ=255)
                                                                            gacgCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCcgcg                                                                       >  8:92006/1‑149 (MQ=255)
                                                                             acccccGAAGCCTGCTCCTGGGGAGCCGCCTAGCCACGCTGTTCGACGAGCAGGGCCCACTGGCC‑‑‑GA‑‑‑GA‑‑TGGGCCAGGAAGTACAACGACTTTAACCCGGTACCCCGGAATATCAGCAATACCTGTCGGGGCAACGCTGCACGCCgagc                                                                      >  6:351772/4‑149 (MQ=255)
                                                                                 ccGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATcc                                                                  <  2:384397/149‑1 (MQ=255)
                                                                                      gCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGaa                                                              <  2:351118/148‑1 (MQ=255)
                                                                                       ccTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGAGGAACGCCGCGCATCCCGAagt                                                            >  4:47843/1‑149 (MQ=255)
                                                                                        ctgctgctgGGCAGCCGCCTTGCCAGGGTTTTCCACGCGCTGGTCACACTGGCCTCACATCCATCCTACACCCGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATc                                                                   <  1:126379/149‑1 (MQ=255)
                                                                                        ctgctgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCaa                                                                                                    >  2:3454/1‑108 (MQ=255)
                                                                                        ctgctgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCaa                                                                                                    <  1:3454/108‑1 (MQ=255)
                                                                                           ctgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtag                                                          <  1:313634/147‑1 (MQ=255)
                                                                                            tgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagcag                                                       >  3:333867/1‑149 (MQ=255)
                                                                                            tgctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTCGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGGATCCCGAagtag                                                          >  4:223119/1‑146 (MQ=255)
                                                                                              ctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagcag                                                       <  3:87489/147‑1 (MQ=255)
                                                                                              ctgGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagcag                                                       >  3:155864/1‑147 (MQ=255)
                                                                                                 ggCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagtagcagtagc                                                   <  8:281399/148‑1 (MQ=255)
                                                                                                   cAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGcc                                                <  7:92006/149‑1 (MQ=255)
                                                                                                      ccgccTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagt                                                            >  2:209820/1‑134 (MQ=255)
                                                                                                      ccgccTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAagt                                                            <  1:209820/134‑1 (MQ=255)
                                                                                                        gccTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           <  7:105288/149‑1 (MQ=255)
                                                                                                          cTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGc                                            >  6:378486/1‑146 (MQ=255)
                                                                                                               cAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAAccc                                    >  1:98620/1‑149 (MQ=255)
                                                                                                                  ggTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCgg                                 >  2:311769/1‑149 (MQ=255)
                                                                                                                       ttCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAg                             >  2:354904/1‑149 (MQ=255)
                                                                                                                         cgacgaGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGc                            <  2:370004/148‑1 (MQ=255)
                                                                                                                            cgaGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCGAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCtttt                        <  8:97522/149‑1 (MQ=255)
                                                                                                                            cgaGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCGGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATGCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCtttt                        <  2:129200/149‑1 (MQ=255)
                                                                                                                              aGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTgc                      <  4:155864/149‑1 (MQ=255)
                                                                                                                                 aGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           >  7:294909/1‑125 (MQ=255)
                                                                                                                                 aGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           >  7:289702/1‑125 (MQ=255)
                                                                                                                                 aGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           <  8:289702/125‑1 (MQ=255)
                                                                                                                                 aGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGcc                                           <  8:294909/125‑1 (MQ=255)
                                                                                                                                   ggCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTg                 <  6:353725/149‑1 (MQ=255)
                                                                                                                                     attgcTGGCC‑‑‑CA‑‑‑AT‑‑TCGGCCAGCAATTACAACGAATTTATCCCGGTAAGCCGGAATATCAGCAATACCTGTCGGATCAACGATGAACGCCGAGCATACCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGAGCCTGGc               >  2:147365/5‑149 (MQ=255)
                                                                                                                                      acacTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGcc              >  8:239272/1‑149 (MQ=255)
                                                                                                                                            gCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTga         <  5:378486/148‑1 (MQ=255)
                                                                                                                                                   a‑‑‑aT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAAc    <  4:50203/149‑1 (MQ=255)
                                                                                                                                                        t‑‑tGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACgg  <  4:333867/149‑1 (MQ=255)
                                                                                                                                                            |                                                                                                                                                  
TGAGTTGTCCTACGAAAACAGGCATTATGTGCCGGCCCGGGACTACTACGATCGTTTCAGCCAGTTGAGCGACCACGACGCCCGTAGCCTGCTGCTGGGCAGCCGCCTTGCCAGGGTGTTCGACGAGCAGGGCACACTGGCC‑‑‑GA‑‑‑AT‑‑TGGGCCAGCAATTACAACGACTTTATCCCGGTACGCCGGAATATCAGCAATACCTGTCGGAGCAACGATGAACGCCGCGCATCCCGAAGTAGCAGTAGCGCCTGGCCAGAACCCCGGTGAGCTTTTGCGCCTGGCCCGTGAGAAACGG  >  NC_002947/985790‑986083

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: