Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,949,820 A→T 7.3% I190N (ATT→AAT)  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,8200AT 7.3% 280.5 / 8.6 94I190N (ATT→AAT) fliMflagellar biosynthesis switching/energizing protein
Reads supporting (aligned to +/- strand):  ref base A (39/47);  new base T (4/3);  total (44/50)
Fisher's exact test for biased strand distribution p-value = 7.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGTCCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTT  >  NC_002947/4949672‑4949967
                                                                                                                                                    |                                                                                                                                                   
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                                                                                                                         <  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                                                                                                                 <  2:1975804/109‑1 (MQ=255)
                                                                            tGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGt                                                                                                                 >  1:1975804/1‑109 (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                                                                   <  1:1988107/129‑1 (MQ=255)
                                                                                                      ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTc                                                                   >  2:1988107/1‑129 (MQ=255)
                                                                                                      ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCAc                                              >  1:2139329/1‑150 (MQ=255)
                                                                                                      ttCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCAc                                              >  1:372566/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                                                       <  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                                                      <  1:753796/104‑1 (MQ=255)
                                                                                                                                          atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                      <  1:753246/104‑1 (MQ=255)
                                                                                                                                          atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                      >  2:753246/3‑106 (MQ=255)
                                                                                                                                          atggCCAACATTGTCGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTgg                                                      >  2:753796/3‑106 (MQ=255)
                                                                                                                                                  cAATGTTGGCCATGGCCGGGTTGACCTCGGCGTCGATATACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCtt  <  2:2410396/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTCGTCCAGGTCCGACTGGAAACCGGCATCGAGCATTTCTCGCACCGGCTCAATCATCGAGTACGGCATGGTCACGTGCAGATCGCCGCCGCCACCGTCCAGTTCGATGTGGAAGGTGGACACCACCACCGCTTCGCTGGGGCCGACAATGTTGGCCATGGCCGGGTTGACCTCGGAGTTCATGTACTCGAAGCTGACCGGCATGATCGCCTGCCAGGCTTCCTTGAGGTCGACGAAACACTGGTCCAGCACCATGCGCACCACACGCAGTTCGGTCGGGGTGAACTCACGCCCTT  >  NC_002947/4949672‑4949967

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: