Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,949,816 G→A 6.3% V191V (GTC→GTT 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,8160GA 6.3% 275.8 / 6.9 95V191V (GTC→GTTfliMflagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base G (42/47);  new base A (3/3);  total (45/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

CCTGGTCGTCCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACT  >  NC_002947/4949668‑4949958
                                                                                                                                                    |                                                                                                                                              
ccTGGTCGTCCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGa                                                                                                                                               <  1:2349567/150‑1 (MQ=255)
 cTGGTCGTCCATGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGAc                                                                                                                                              >  2:693352/1‑150 (MQ=255)
    gtcgtcCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAAt                                                                                                                                           >  1:2008011/1‑150 (MQ=255)
    gtcgtcCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAAt                                                                                                                                           >  2:2272729/1‑150 (MQ=255)
             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                                                                                                                      >  1:139575/1‑150 (MQ=255)
                         aCCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGAc                                                                                                                      <  2:2081830/150‑1 (MQ=255)
                           cGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCt                                                                                                                    >  2:1952104/1‑150 (MQ=255)
                            ggCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCCCCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGCCCTc                                                                                                                   <  2:644943/150‑1 (MQ=255)
                               aTCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGa                                                                                                                <  2:1553835/150‑1 (MQ=255)
                               aTCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGa                                                                                                                <  1:526644/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:693352/150‑1 (MQ=255)
                                          ctcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACt                                                                                                     <  1:1952104/150‑1 (MQ=255)
                                          ctcGCACCGGCTCAATCATCGAGTACGGCATGGACACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTGCATGTACt                                                                                                     >  1:2410396/1‑150 (MQ=255)
                                           tcGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTc                                                                                                    >  2:1292068/1‑150 (MQ=255)
                                                    cTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTgg                                                                                                                                      >  1:30778/1‑107 (MQ=255)
                                                    cTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTgg                                                                                                                                      <  2:30778/107‑1 (MQ=255)
                                                     tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTc                                                                                                            <  1:141759/132‑1 (MQ=255)
                                                     tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTc                                                                                                            >  2:141759/1‑132 (MQ=255)
                                                     tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGAcc                                                                                          <  1:1292068/150‑1 (MQ=255)
                                                     tCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGAcc                                                                                          <  2:1871184/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                                                                                                    >  1:900912/1‑118 (MQ=255)
                                                                           tCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTc                                                                                                    <  2:900912/118‑1 (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:2272729/150‑1 (MQ=255)
                                                                                              gccACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAaca                                                 >  1:1389067/1‑150 (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                                                          <  1:1988107/129‑1 (MQ=255)
                                                                                                          ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTc                                                          >  2:1988107/1‑129 (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                                                                                 >  1:905006/1‑89 (MQ=255)
                                                                                                                           acacCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATc                                                                                 <  2:905006/89‑1 (MQ=255)
                                                                                                                           acacCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCa                    >  1:528519/1‑150 (MQ=255)
                                                                                                                            caccaccacCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCa                                   >  2:2400249/1‑134 (MQ=255)
                                                                                                                            caccaccacCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCa                                   <  1:2400249/134‑1 (MQ=255)
                                                                                                                             accaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTg                                              <  2:528633/122‑1 (MQ=255)
                                                                                                                             accaccaccGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTg                                              >  1:528633/1‑122 (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:1559577/150‑1 (MQ=255)
                                                                                                                                           gCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGaa    >  2:1801068/1‑150 (MQ=255)
                                                                                                                                             tGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCAAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACt  >  2:1485149/1‑150 (MQ=255)
                                                                                                                                              atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                             >  2:753246/3‑106 (MQ=255)
                                                                                                                                              atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                             >  2:753796/3‑106 (MQ=255)
                                                                                                                                              atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                             <  1:753246/104‑1 (MQ=255)
                                                                                                                                              atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                             <  1:753796/104‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
CCTGGTCGTCCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACT  >  NC_002947/4949668‑4949958

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: