Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,949,824 T→C 6.2% N189D (AAC→GAC)  fliM ← flagellar biosynthesis switching/energizing protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,949,8240TC 6.2% 284.5 / 5.9 96N189D (AAC→GAC) fliMflagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base T (42/48);  new base C (3/3);  total (45/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

GTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATC  >  NC_002947/4949681‑4949972
                                                                                                                                               |                                                                                                                                                    
gTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCAt                                                                                                                                                >  1:1553835/1‑150 (MQ=255)
gTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCAACATTGTCGGcccc                                                                                                                                                >  2:888374/1‑148 (MQ=255)
 tCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCAACATTGTCGGcccca                                                                                                                                               <  1:888374/150‑4 (MQ=255)
   cGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGc                                                                                                                                             >  2:1642903/1‑150 (MQ=255)
         gAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGtt                                                                                                                                       >  2:1713999/1‑150 (MQ=255)
           aaCCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGGCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGa                                                                                                                                     >  1:770170/1‑150 (MQ=255)
            aCCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGAc                                                                                                                                    <  2:2081830/150‑1 (MQ=255)
            aCCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGAc                                                                                                                                    >  1:139575/1‑150 (MQ=255)
              cGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCt                                                                                                                                  >  2:1952104/1‑150 (MQ=255)
               ggCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCCCCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGCCCTc                                                                                                                                 <  2:644943/150‑1 (MQ=255)
                  aTCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGa                                                                                                                              <  1:526644/150‑1 (MQ=255)
                  aTCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGa                                                                                                                              <  2:1553835/150‑1 (MQ=255)
                  aTCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGGGGAAGTTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCGGTTTTGACCTCGGa                                                                                                                              <  2:938063/150‑1 (MQ=255)
                  aTCGAGCATTTATCGCACCGACTCAATCATCGATTACGGCATGGTCACGTGCAGATCGCCGGCGCGATGTTTGTTTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGa                                                                                                                              <  1:1861097/150‑1 (MQ=255)
                     gAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGTCGACCATTTTGGCCATGGCCGGGTTGACCTCGGAGtt                                                                                                                           >  1:854268/1‑150 (MQ=255)
                      aGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGTATTTc                                                                                                                          >  1:1294310/1‑150 (MQ=255)
                       gCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCa                                                                                                                         <  1:2241056/150‑1 (MQ=255)
                        cATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCAt                                                                                                                        <  1:1082563/150‑1 (MQ=255)
                        cATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCAt                                                                                                                        <  1:1783540/150‑1 (MQ=255)
                            tctcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCTGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTAc                                                                                                                    <  1:1729279/150‑1 (MQ=255)
                             ctcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACt                                                                                                                   <  1:1952104/150‑1 (MQ=255)
                             ctcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACt                                                                                                                   <  1:693352/150‑1 (MQ=255)
                             ctcGCACCGGCTCAATCATCGAGTACGGCATGGACACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTGCATGTACt                                                                                                                   >  1:2410396/1‑150 (MQ=255)
                              tcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTc                                                                                                                  >  2:1292068/1‑150 (MQ=255)
                                       cTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTgg                                                                                                                                                    <  2:30778/107‑1 (MQ=255)
                                       cTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTgg                                                                                                                                                    >  1:30778/1‑107 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTc                                                                                                                          >  2:141759/1‑132 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTc                                                                                                                          <  1:141759/132‑1 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGAcc                                                                                                        <  2:1871184/150‑1 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGAcc                                                                                                        <  1:1292068/150‑1 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAAGGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATg                                                                                                                       >  1:144340/1‑135 (MQ=255)
                                        tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCCTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATg                                                                                                                       <  2:144340/135‑1 (MQ=255)
                                          aaTCATCGAGTACGGCATGGTCACGTGGAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTGCATGTACTCGAAGCTGACCgg                                                                                                      >  1:1775081/1‑150 (MQ=255)
                                          aaTCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCgg                                                                                                                               <  1:641850/125‑1 (MQ=255)
                                          aaTCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCgg                                                                                                                               >  2:641850/1‑125 (MQ=255)
                                           atcatcGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCAACGCAGCTTGTCAAGGTGGAAGGTGGACACCACTTCCGCGTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGc                                                                                                     <  1:1116735/150‑1 (MQ=255)
                                            tcatcGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCa                                                                                                    >  2:474173/1‑150 (MQ=255)
                                            tcatcGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTGCATGTACTCGAAGCTGACCGGCa                                                                                                    >  1:1633120/1‑150 (MQ=255)
                                               tcGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGGTTCGCTGGGGCCGACAATGATGGCCATCACCGGGTTGACCTGGCGCATCATGTACTCGAAGCTGACCGGCATGa                                                                                                 <  1:1882008/150‑1 (MQ=255)
                                               tcGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGTCCGACATTGTTTGCCATGGCCGGTGTTACCTCGGTGTTTGTCTACGCGAAGCTGACCGGCATGa                                                                                                 >  2:1706916/1‑150 (MQ=255)
                                                 gAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATc                                                                                               <  1:488407/150‑1 (MQ=255)
                                                     aCGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCt                                                                                           <  2:1235856/150‑1 (MQ=255)
                                                        gCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGTGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCCTGTACTCGAAGCTGCCCGGCATGCTCGCCTGcc                                                                                        >  2:436717/1‑150 (MQ=255)
                                                           tGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGGTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAgg                                                                                     <  1:748404/150‑1 (MQ=255)
                                                           tGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGAACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAgg                                                                                     <  1:745235/150‑1 (MQ=255)
                                                            ggTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGc                                                                                    <  1:1063005/150‑1 (MQ=255)
                                                             gTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCt                                                                                   <  1:1713999/150‑1 (MQ=255)
                                                              tCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTc                                                                                                                  <  2:900912/118‑1 (MQ=255)
                                                              tCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTc                                                                                                                  >  1:900912/1‑118 (MQ=255)
                                                                   tGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGt                                                                                                                      >  1:1975804/1‑109 (MQ=255)
                                                                   tGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGt                                                                                                                      <  2:1975804/109‑1 (MQ=255)
                                                                       gATCGCCGGCGCCACCGTCCAGTGCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCAGATAATTTTGTCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACAGGCATGATCGCCTGCCAGGCTTCCTTGAGGt                                                                         <  1:436717/150‑1 (MQ=255)
                                                                           gccgccgccACCGTCCAGTTCGATGTGGAAGGTGGATCGGACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTcgac                                                                     <  1:1706916/150‑1 (MQ=255)
                                                                               ccgccACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGaaa                                                                 <  2:2008011/150‑1 (MQ=255)
                                                                                 gccACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAaca                                                               >  1:1389067/1‑150 (MQ=255)
                                                                                 gccACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAaca                                                               <  1:2272729/150‑1 (MQ=255)
                                                                                   cACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACt                                                             <  2:347966/150‑1 (MQ=255)
                                                                                    aCCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTg                                                            <  1:1642903/150‑1 (MQ=255)
                                                                                       gTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGTCTTCCTTGAGGTCGACGAAACACTGGTc                                                         <  2:1077441/150‑1 (MQ=255)
                                                                                       gTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTc                                                         >  2:1838181/1‑150 (MQ=255)
                                                                                             ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTc                                                                        >  2:1988107/1‑129 (MQ=255)
                                                                                             ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTc                                                                        <  1:1988107/129‑1 (MQ=255)
                                                                                             ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCAc                                                   >  1:372566/1‑150 (MQ=255)
                                                                                             ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCAc                                                   >  1:2139329/1‑150 (MQ=255)
                                                                                                gATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCAt                                                <  2:372566/150‑1 (MQ=255)
                                                                                                 aTGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATg                                               <  1:1838181/150‑1 (MQ=255)
                                                                                                      gAAGGTGGACACCAGCACCGCTTCGCTGGCGCCGACACTTTTGGCCATGGCGGGCTTGACCTCGGAGTTCATGTGCTCTAAGCTGCGCGGCACGATCGCCAGCGAGGCATCCTTGAGGAAGACGAAACACTGGTCCAGCACCATGCGcac                                          >  2:94233/1‑150 (MQ=255)
                                                                                                      gAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATg                                                                                                                       <  2:721120/73‑1 (MQ=255)
                                                                                                      gAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATg                                                                                                                       >  1:721120/1‑73 (MQ=255)
                                                                                                        aGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGcacca                                        >  2:538280/1‑150 (MQ=255)
                                                                                                          gTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTAATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCAt                                                <  1:1626029/140‑1 (MQ=255)
                                                                                                          gTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTAATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCAt                                                >  2:1626029/1‑140 (MQ=255)
                                                                                                              acacCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATc                                                                                               <  2:905006/89‑1 (MQ=255)
                                                                                                              acacCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATc                                                                                               >  1:905006/1‑89 (MQ=255)
                                                                                                              acacCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCa                                  >  1:528519/1‑150 (MQ=255)
                                                                                                               caccaccacCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCa                                                 <  1:2400249/134‑1 (MQ=255)
                                                                                                               caccaccacCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCa                                                 >  2:2400249/1‑134 (MQ=255)
                                                                                                                accaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTg                                                            >  1:528633/1‑122 (MQ=255)
                                                                                                                accaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTg                                                            <  2:528633/122‑1 (MQ=255)
                                                                                                                 ccaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGtt                               <  2:139575/150‑1 (MQ=255)
                                                                                                                 ccaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGAGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTAGTCCAGCACCATGCGCACCACACGCAGtt                               >  1:483569/1‑150 (MQ=255)
                                                                                                                  caccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTtc                              >  2:1371384/1‑150 (MQ=255)
                                                                                                                       ccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTtcggtcg                         >  1:545195/1‑150 (MQ=255)
                                                                                                                           ttCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGt                     <  1:538280/150‑1 (MQ=255)
                                                                                                                              gCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGaa                  >  2:1801068/1‑150 (MQ=255)
                                                                                                                              gCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGaa                  <  2:1559577/150‑1 (MQ=255)
                                                                                                                                tGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCAAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACt                >  2:1485149/1‑150 (MQ=255)
                                                                                                                                 atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                           <  1:753246/104‑1 (MQ=255)
                                                                                                                                 atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                           >  2:753246/3‑106 (MQ=255)
                                                                                                                                 atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                           >  2:753796/3‑106 (MQ=255)
                                                                                                                                 atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                           <  1:753796/104‑1 (MQ=255)
                                                                                                                                         cAATGTTGGCCATGGCCGGGTTGACCTCGGCGTCGATATACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCtt       <  2:2410396/150‑1 (MQ=255)
                                                                                                                                            tGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGTGGGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGa    <  1:982777/150‑1 (MQ=255)
                                                                                                                                            tGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGa    >  2:882769/1‑150 (MQ=255)
                                                                                                                                            tGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGa    >  1:250827/1‑150 (MQ=255)
                                                                                                                                            tGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGACATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGa    >  2:2316657/1‑150 (MQ=255)
                                                                                                                                              ttGGCCATGGCCGGGTTTACCTCGGAGTTCAGGTACTCGAAGCTGACCGGCATGATCGCCTGCCAACAATCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATc  <  1:748371/150‑1 (MQ=255)
                                                                                                                                              ttGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATc  >  2:859061/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTTCGATC  >  NC_002947/4949681‑4949972

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