Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,105,663 T→C 6.6% S817G (AGC→GGC)  ydiJ ← putative oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,105,6630TC 6.6% 222.6 / 4.4 76S817G (AGC→GGC) ydiJputative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (37/34);  new base C (3/2);  total (40/36)
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.35e-01

GCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCC  >  NC_002947/5105517‑5105811
                                                                                                                                                   |                                                                                                                                                    
gCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTg                                                                                                                                                     >  1:331736/1‑148 (MQ=255)
  gATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTa                                                                                                                                                   <  3:73111/148‑1 (MQ=255)
  gATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTa                                                                                                                                                   >  4:73111/1‑148 (MQ=255)
     acacCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgc                                                                                                                                               >  5:106892/1‑149 (MQ=255)
       acCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcgc                                                                                                                                             <  2:331736/149‑1 (MQ=255)
         cAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGgcgc                                                                                                                                             >  6:116259/1‑147 (MQ=255)
          aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCATTCGGCCAGGTCCTTGAGCTGCGTGTCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                          <  7:190468/149‑1 (MQ=255)
          aGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAg                                                                                                                                          <  6:83133/149‑1 (MQ=255)
                aTCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAc                                                                                                                                    >  2:66589/1‑149 (MQ=255)
                 tCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAcc                                                                                                                                   >  8:154629/1‑149 (MQ=255)
                   gCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAccc                                                                                                                                  <  4:321935/148‑1 (MQ=255)
                          ccAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAccc                                                                                                                                  <  6:354884/141‑1 (MQ=255)
                          ccAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATAccc                                                                                                                                  >  5:354884/1‑141 (MQ=255)
                            aGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAg                                                                                                                        >  5:275049/1‑149 (MQ=255)
                              ccccACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGCGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAgg                                                                                                                       <  1:66589/149‑1 (MQ=255)
                                 caccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGt                                                                                                                   >  3:169103/1‑149 (MQ=255)
                                  accaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTg                                                                                                                  >  5:212350/1‑149 (MQ=255)
                                   ccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                 <  6:191851/149‑1 (MQ=255)
                                   ccaGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                 <  3:87693/149‑1 (MQ=255)
                                         gCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaacg                                                                                                                                        >  4:407704/1‑117 (MQ=255)
                                         gCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCCGTACAGCGCCaacg                                                                                                                                        <  3:407704/120‑4 (MQ=255)
                                             acCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGa                                                                                                       <  1:1693/149‑1 (MQ=255)
                                                 cAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATg                                                                                                     >  8:190395/1‑147 (MQ=255)
                                                 cAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                   <  7:306155/149‑1 (MQ=255)
                                                  aaTCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAg                                                                                                  >  5:295429/1‑149 (MQ=255)
                                                       gCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGa                                                                                             >  6:173731/1‑149 (MQ=255)
                                                               cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                 >  3:401136/1‑121 (MQ=255)
                                                               cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTgg                                                                                                                 <  4:401136/121‑1 (MQ=255)
                                                               cTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAgcg                                                                                     <  6:106892/149‑1 (MQ=255)
                                                                    gCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGtt                                                                                <  4:77291/149‑1 (MQ=255)
                                                                       gCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTc                                                                              >  3:346210/1‑148 (MQ=255)
                                                                        cGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCg                                                                             >  2:318319/1‑148 (MQ=255)
                                                                         gTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                           <  6:295429/149‑1 (MQ=255)
                                                                         gTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                           <  6:212350/149‑1 (MQ=255)
                                                                              gTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGc                                                                      >  3:369011/1‑149 (MQ=255)
                                                                                tGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    <  7:154629/149‑1 (MQ=255)
                                                                                 gCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                           >  7:162906/1‑141 (MQ=255)
                                                                                 gCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGaa                                                                           <  8:162906/141‑1 (MQ=255)
                                                                                 gCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    <  5:173731/148‑1 (MQ=255)
                                                                                  cGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCggg                                                                   <  4:346210/148‑1 (MQ=255)
                                                                                       cGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAgg                                                             >  3:212826/1‑149 (MQ=255)
                                                                                         cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGc                                                                                                >  4:344073/1‑112 (MQ=255)
                                                                                         cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCg                                                           <  7:190395/149‑1 (MQ=255)
                                                                                         cTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCg                                                           >  5:408917/1‑149 (MQ=255)
                                                                                            ccTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTcc                                                        >  7:95142/1‑149 (MQ=255)
                                                                                            ccTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTcc                                                        >  5:302593/1‑149 (MQ=255)
                                                                                              ttGGCACACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGc                                                            <  2:209649/143‑1 (MQ=255)
                                                                                              ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGc                                                            >  1:209649/1‑143 (MQ=255)
                                                                                              ttGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTg                                                      >  8:237801/1‑149 (MQ=255)
                                                                                                ggCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTg                                                      <  4:212826/147‑1 (MQ=255)
                                                                                                 gCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    <  7:77258/132‑1 (MQ=255)
                                                                                                 gCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCgg                                                                    >  8:77258/1‑132 (MQ=255)
                                                                                                            aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                   <  5:293521/90‑1 (MQ=255)
                                                                                                            aaGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGaa                                                                                                   >  6:293521/1‑90 (MQ=255)
                                                                                                              ggAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCAcc                                                  <  7:361381/137‑1 (MQ=255)
                                                                                                              ggAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCAcc                                                  >  8:361381/1‑137 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAgg                                                                                                 <  2:108010/89‑1 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAgg                                                                                                 >  1:108010/1‑89 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGt                                                                  <  3:149640/120‑1 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGt                                                                  >  4:149640/1‑120 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGaa                                                               <  4:129304/123‑1 (MQ=255)
                                                                                                               gAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGaa                                                               >  3:129304/1‑123 (MQ=255)
                                                                                                                    ccTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTcg                                <  6:275049/149‑1 (MQ=255)
                                                                                                                    ccTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGATct                                >  7:343335/1‑148 (MQ=255)
                                                                                                                           gTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTggg                         >  1:248947/1‑149 (MQ=255)
                                                                                                                             gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGc                                                     >  3:317176/1‑119 (MQ=255)
                                                                                                                             gCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGc                                                     <  4:317176/119‑1 (MQ=255)
                                                                                                                              cAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCg                      >  8:225156/1‑149 (MQ=255)
                                                                                                                                aCGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCggg                    <  3:289369/149‑1 (MQ=255)
                                                                                                                                          gTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTc                                                         >  4:128901/1‑102 (MQ=255)
                                                                                                                                          gTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTc                                                         <  3:128901/102‑1 (MQ=255)
                                                                                                                                            tGGCGCCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTgg                          >  2:4977/1‑131 (MQ=255)
                                                                                                                                            tGGCGCCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTgg                          <  1:4977/131‑1 (MQ=255)
                                                                                                                                            ctcgcgcCGTACAGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATCAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGAGAAGGCGTCCTGCACCAGCACGATGCTGCCCTCGCGCTGGGCCGGGGTCAGCTcgcg         >  4:74792/4‑149 (MQ=255)
                                                                                                                                               gcgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTcgcg         <  8:317787/146‑1 (MQ=255)
                                                                                                                                               gcgcTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTcgcg         >  7:317787/1‑146 (MQ=255)
                                                                                                                                                   tGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGcc  >  3:392638/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCGATACACCAGCGTCATCGCCGGGTCCAGCCCCACCAGCGGCACACCGCAATCGGCCAGGGCCTTGAGCTGCGTGGCGTTGCGGATCGCTGCCTTGGCAAACGCGCCAAGGAAACCCTGTACGTGCAACGGCTTGCCGTTGGCGCTGTACGGCGCCAGGAATACCCGGTGGCCCAGGCGGTGGGCAAGGTCGATGAAGGCCGACAGCAGCGGCGTTTCGAAGTAGCGGGTGAAGGCGTCCTGCACCAGCACGATGCTGCGCTCGCGCTGGGCCGGGGTCAGCTCGCGCAGGGCC  >  NC_002947/5105517‑5105811

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: