Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,945,056 C→T 6.1% intergenic (+371/+427) PP_5673 → / ← flhB D‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,945,0560CT 6.1% 230.5 / 2.6 80intergenic (+371/+427)PP_5673/flhBD‑alanine‑‑D‑alanine ligase/flagellin export apparatus substrate specificity protein
Reads supporting (aligned to +/- strand):  ref base C (38/36);  new base T (2/3);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTG  >  NC_002947/4944909‑4945201
                                                                                                                                                   |                                                                                                                                                 
gCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCg                                                                                                                                                 <  2:1506830/150‑1 (MQ=39)
gCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCg                                                                                                                                                 <  1:2047419/150‑1 (MQ=39)
  gTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGc                                                                                                                                               >  1:1533146/1‑150 (MQ=39)
  gTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGc                                                                                                                                               <  1:908908/150‑1 (MQ=38)
       cacaGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGGGCTCGGccggg                                                                                                                                           >  1:1437152/1‑146 (MQ=255)
       cacaGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGGGCTCGGCCggg                                                                                                                                           <  2:1437152/149‑1 (MQ=255)
        acaGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGgtt                                                                                                                                         >  1:753209/1‑150 (MQ=38)
        acaGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGgtt                                                                                                                                         >  2:1510107/1‑150 (MQ=39)
           gACGAGCAGCACCAGGATCGACAGGGCGATCTAGAAGGTGAAAAACACGTGGAGCCGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAGCAGTGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGgttgtt                                                                                                                                      <  2:802287/150‑1 (MQ=11)
            aCCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGttgttg                                                                                                                                     >  1:1510583/1‑150 (MQ=39)
             ccAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGc                                                                                                                                    >  1:1630528/1‑150 (MQ=39)
              cagcagCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCg                                                                                                                                   <  2:753209/150‑1 (MQ=39)
              cagcagCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCg                                                                                                                                   <  2:717339/150‑1 (MQ=39)
               agcagcACCTCGATCGACAGCGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAAGATCATCATTTTCAACATCGCCAGCAGTAACAGCGGCACACCCAGGTTGCCCCCGCCTAGCCCGGCCAGGTTGTTGCgg                                                                                                                                  <  1:893428/150‑1 (MQ=11)
               agcagcACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCgg                                                                                                                                  <  1:1510107/150‑1 (MQ=38)
                       cacGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCt                                                                                                                          >  2:779042/1‑150 (MQ=39)
                       cacGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGCCGGTATCGGCAACATCATCATGGCCGACATCCCCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCt                                                                                                                          >  2:686495/1‑150 (MQ=255)
                          gATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAt                                                                                                                       <  2:289407/150‑1 (MQ=38)
                           aTCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAtt                                                                                                                      >  1:1212432/1‑150 (MQ=37)
                              gACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAAc                                                                                                                   >  2:1999921/1‑150 (MQ=39)
                                 aGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACGAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGa                                                                                                                <  1:1364125/150‑1 (MQ=37)
                                 aGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGa                                                                                                                <  2:1533146/150‑1 (MQ=38)
                                       aTGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGa                                                                                                          <  2:1630528/150‑1 (MQ=39)
                                       aTGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGa                                                                                                          <  1:1385016/150‑1 (MQ=39)
                                            gAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCAc                                                                                                     >  1:401466/1‑150 (MQ=38)
                                                 tGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGa                                                                                                <  2:1182299/150‑1 (MQ=34)
                                                     gaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTcc                                                                                            >  2:875725/1‑150 (MQ=35)
                                                     gaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTcc                                                                                            >  1:1102569/1‑150 (MQ=35)
                                                     gaaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCAGCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTcc                                                                                            >  1:801094/1‑150 (MQ=25)
                                                      aaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGa                                                                                                                        <  1:996374/121‑1 (MQ=37)
                                                      aaCACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGa                                                                                                                        >  2:996374/1‑121 (MQ=37)
                                                        cacGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCAAGCGCCCCGAGGTTGTTGTGGGTGTTGCTGATTAACTGAGTGCGATGCACCGCGAGTCctct                                                                                         <  2:84022/150‑1 (MQ=11)
                                                              gAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAg                                                                                                                                            <  1:1655800/93‑1 (MQ=39)
                                                              gAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAg                                                                                                                                            >  2:1655800/1‑93 (MQ=39)
                                                                  aGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCAc                                                                                                     <  1:114681/128‑1 (MQ=37)
                                                                  aGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCAc                                                                                                     >  2:114681/1‑128 (MQ=37)
                                                                     aaCGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACtt                                                                            <  2:1102569/150‑1 (MQ=34)
                                                                           ggTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGc                                                                      >  2:1039313/1‑150 (MQ=35)
                                                                                 ggCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTg                                                                                                             >  2:2387048/1‑105 (MQ=38)
                                                                                 ggCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTg                                                                                                             <  1:2387048/105‑1 (MQ=38)
                                                                                  gaaaCATAATCATTGCCAACAGCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAgggg                                                               >  2:643118/3‑150 (MQ=19)
                                                                                          atcatTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGc                                                       >  2:992634/1‑150 (MQ=37)
                                                                                           tcatTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGGCTGCg                                                      >  1:1548029/1‑150 (MQ=37)
                                                                                             atTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAt                                                                                                                       <  1:368202/83‑1 (MQ=39)
                                                                                             atTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAt                                                                                                                       >  2:368202/1‑83 (MQ=40)
                                                                                             atTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGc                                                    >  1:2271978/1‑150 (MQ=39)
                                                                                              ttGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCt                                                   <  1:875725/150‑1 (MQ=39)
                                                                                                    aCATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCAtt                                             <  2:1510583/150‑1 (MQ=255)
                                                                                                     cATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCct                                                                                           <  2:1899740/103‑1 (MQ=38)
                                                                                                     cATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGGCct                                                                                           >  1:1899740/1‑103 (MQ=38)
                                                                                                      aTCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTGTTGATGCAAAACTTTGACGCACCGGGGGCGCCTGCGGCTGGCATTGc                                           >  2:183367/1‑150 (MQ=255)
                                                                                                      aTCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGc                                           <  2:2271978/150‑1 (MQ=255)
                                                                                                        cACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTa                                         <  2:401466/150‑1 (MQ=255)
                                                                                                         aCCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAg                                        <  1:1534874/150‑1 (MQ=255)
                                                                                                         aCCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAg                                        >  1:2141009/1‑150 (MQ=255)
                                                                                                          ccAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGa                                       >  1:1862565/1‑150 (MQ=255)
                                                                                                                gTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCttt                                 <  2:801094/150‑1 (MQ=255)
                                                                                                                    cAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCaa                             <  1:1039313/150‑1 (MQ=255)
                                                                                                                     aGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAAc                            >  1:2064535/1‑150 (MQ=255)
                                                                                                                     aGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAAc                            >  1:2279748/1‑150 (MQ=255)
                                                                                                                      gCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGAGTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACt                           >  2:2095897/1‑150 (MQ=255)
                                                                                                                        ggCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACttt                         <  1:779042/150‑1 (MQ=255)
                                                                                                                         gCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTtg                        <  2:2141009/150‑1 (MQ=255)
                                                                                                                          cacaCCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGCGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTtgt                       <  1:802741/150‑1 (MQ=255)
                                                                                                                              cccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATc                   >  2:2080111/1‑150 (MQ=255)
                                                                                                                               ccAGGTTGCCCCGGCCGAGCTTGGCCAGGTTGTTGCGGGCGTTGCTTATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGGAAAACTATGACGCCCAGGGGGCGCCGGCGGCTGGCATTGCTAGAAGCTTCCCAACTTTGTGATct                  <  2:1548029/150‑1 (MQ=14)
                                                                                                                               ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATct                  >  2:1635392/1‑150 (MQ=255)
                                                                                                                               ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATct                  >  1:1628414/1‑150 (MQ=255)
                                                                                                                               ccAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATct                  >  1:1464839/1‑150 (MQ=255)
                                                                                                                                        cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGaa                                      <  1:1245709/121‑1 (MQ=255)
                                                                                                                                        cccGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGaa                                      >  2:1245709/1‑121 (MQ=255)
                                                                                                                                         ctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTtgt                       >  1:1186861/3‑135 (MQ=255)
                                                                                                                                         ctggCCGGGCTCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTtgt                       <  2:1186861/133‑1 (MQ=255)
                                                                                                                                           ggCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCgg      <  2:1628414/150‑1 (MQ=255)
                                                                                                                                            gCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCggg     <  2:1862565/150‑1 (MQ=255)
                                                                                                                                              cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCt                                                   <  2:325956/102‑1 (MQ=35)
                                                                                                                                              cGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCt                                                   >  1:325956/1‑102 (MQ=35)
                                                                                                                                               gAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGgtg  >  1:361816/1‑150 (MQ=255)
                                                                                                                                                 gCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGc                                           >  2:1735144/1‑107 (MQ=255)
                                                                                                                                                 gCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGc                                           <  1:1735144/107‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCGTATACACAGACCAGCAGCACCACGATCGACAGGGCGATGTTGAAGGTGAAGAACACGTCGAGCAGGAACGGCGGTATCGGCAACATCATCATTGCCAACATCACCAGCAGTAACAGCGGCACACCCAGGTTGCCCCGGCCGAGCCCGGCCAGGTTGTTGCGGGCGTTGCTGATTAACTGAGTGCGATCCACCGCGAGTCCTCTTGATGCAAAACTTTGACGCCCAGGGGGCGCCTGCGGCTGGCATTGCTAGAAGCTTTCCAACTTTGTGATCTGCTGTGAATCGGGGTG  >  NC_002947/4944909‑4945201

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