Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,973,319 C→A 6.4% E92* (GAG→TAG) ‡ flgL ← flagellar hook‑associated protein FlgL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,973,3190CA 6.4% 351.6 / 3.9 113E92* (GAG→TAG) ‡flgLflagellar hook‑associated protein FlgL
Reads supporting (aligned to +/- strand):  ref base C (55/49);  new base A (5/3);  total (61/52)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 113 total aligned reads displayed.

TTGGAGCCGCCAAACATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAGCGCCGACCGGG  >  NC_002947/4973174‑4973467
                                                                                                                                                 |                                                                                                                                                    
ttGGAGCCGCCAAACATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTg                                                                                                                                                  <  2:1621104/150‑1 (MQ=255)
    aGCCGCCAAACATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCGAGCTGGcga                                                                                                                                              <  2:1092147/150‑3 (MQ=255)
         ccAAACATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAACTCGCTGGCACGCTGc                                                                                                                                         >  2:1152771/1‑150 (MQ=255)
           agacgtGTACTCGCCGTTGGCATCGCGCGATTTCAGCAGACCGAACACGTTGGCTTCGATTTCCATCATCGCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAg                                                                                                                                       <  1:1731557/145‑1 (MQ=255)
           aaaCATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAg                                                                                                                                       <  2:697510/150‑1 (MQ=255)
                    gTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGc                                                                                                                                                    <  1:1053072/128‑1 (MQ=255)
                    gTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGc                                                                                                                                                    >  2:1053072/1‑128 (MQ=255)
                       gCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTg                                                                                                                                                  <  2:148520/127‑1 (MQ=255)
                       gCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTg                                                                                                                                                  >  1:148520/1‑127 (MQ=255)
                        ccGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATg                                                                                                                          <  1:832073/150‑1 (MQ=255)
                          gTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGt                                                                                                                        >  1:1823148/1‑150 (MQ=255)
                          gTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGt                                                                                                                        >  1:100406/1‑150 (MQ=255)
                                aTCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAg                                                                                                                  >  2:1349602/1‑150 (MQ=255)
                                     gcgAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGtt                                                                                                                              <  2:1110409/133‑1 (MQ=255)
                                     gcgAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGtt                                                                                                                              >  1:1110409/1‑133 (MQ=255)
                                     gcgAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGc                                                                                                             >  1:713161/1‑150 (MQ=255)
                                       gtgTTCAGCAGACCGAACAGGTTGGCTTCGATTTCCTTGATCTCGCTTGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGTCTGGATAGCCAGCTCGCAGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGttc                                                                                                           <  2:535109/148‑1 (MQ=255)
                                       gAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGctc                                                                                                           >  2:1228244/1‑150 (MQ=255)
                                       gAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGctc                                                                                                           <  1:1349602/150‑1 (MQ=255)
                                         gTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGctctc                                                                                                         >  2:2245400/1‑150 (MQ=255)
                                              gcagACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCt                                                                                                    >  1:1526242/1‑150 (MQ=255)
                                              gaagAACGAACACACTGGCTACGACTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCCCTCTTCCt                                                                                                    <  2:752306/148‑1 (MQ=255)
                                                  aCCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTTGCGATCGAAGTACGGTCTGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCGCGCTGGCACGCTGCCGGGCACCGTGGATGGTCGATAGCACCCTCTCTTCCtcgag                                                                                                >  2:1104093/1‑146 (MQ=255)
                                                  aCCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTgcag                                                                                                >  2:367719/1‑150 (MQ=255)
                                                    cGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTc                                                                                                                       <  2:2241574/125‑1 (MQ=255)
                                                    cGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTc                                                                                                                       >  1:2241574/1‑125 (MQ=255)
                                                    cGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTgcagca                                                                                              >  1:1160218/1‑150 (MQ=255)
                                                    cGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTgcagca                                                                                              >  2:1978516/1‑150 (MQ=255)
                                                    cGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTgcagca                                                                                              >  2:48779/1‑150 (MQ=255)
                                                     gAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCGCGCGGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGcagcag                                                                                             >  1:1132543/1‑150 (MQ=255)
                                                     gAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTc                                                                                                                       <  1:1811446/124‑1 (MQ=255)
                                                     gAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTc                                                                                                                       >  2:1811446/1‑124 (MQ=255)
                                                      aaCACGTTGGCTTGTTTTTCCTTGATCTCGCTCGCGATCGCGGTACGGTCGGCATCGCTCAGGCCCCCTTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGTATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAgg                                                                                            <  1:2296174/150‑1 (MQ=255)
                                                      aaCACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCACCGCCCAGCCACCCCTCACCCGCCTGGAACCACTGTGATCTGTCACTCTGGATGTCTTCGTTGATTGTCGATAGCATGCTCTCTTCCTGCAGCAgg                                                                                            <  2:276233/150‑1 (MQ=255)
                                                           gTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGGGAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGtt                                                                                       >  1:1354415/1‑150 (MQ=255)
                                                            ttGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTg                                                                                                                                                  <  1:1889115/90‑1 (MQ=255)
                                                            ttGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTg                                                                                                                                                  >  2:1889115/1‑90 (MQ=255)
                                                            ttGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTc                                                                                      >  1:1850304/1‑150 (MQ=255)
                                                               gCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTg                                                                                   >  1:1624997/1‑150 (MQ=255)
                                                               gCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTg                                                                                   <  1:524957/150‑1 (MQ=255)
                                                               gCTACGAGTTCCTTGATCTCGCTGGCGTTCGCGGTGCGGCCGTCATCGCTCAGGCCCACGTTACCCGCCGGGATCGCGATCTCGCTGGCACGCTGCAGGGCATCCTTGTAGATCGGTAGCGCGCCCTCTTCCTTGAGCAGCGCGTGCTTg                                                                                   >  2:1062259/1‑150 (MQ=255)
                                                                cTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGc                                                                                          <  1:1063660/142‑1 (MQ=255)
                                                                cTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAACTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGa                                                                                  <  1:1152771/150‑1 (MQ=255)
                                                                 ttCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAgg                                                                                            >  2:340162/1‑139 (MQ=255)
                                                                 ttCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGAc                                                                                 >  2:347564/1‑150 (MQ=255)
                                                                 ttCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGAc                                                                                 <  1:347564/150‑1 (MQ=255)
                                                                     aTTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTg                                                                             <  2:1823148/150‑1 (MQ=255)
                                                                     aTTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTg                                                                             >  2:2199208/1‑150 (MQ=255)
                                                                        tCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTc                                                                          <  2:551134/150‑1 (MQ=255)
                                                                           ttGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATg                                                                       >  1:1603958/1‑150 (MQ=255)
                                                                            tGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGt                                                                      <  2:1526242/150‑1 (MQ=255)
                                                                             gATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGtt                                                                     >  1:2360797/1‑150 (MQ=255)
                                                                                   gctggcgatccaggcgggtaacgggggcctgagcgatgccgaccgtacctcgatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCt                                                               >  1:1646742‑M2/62‑150 (MQ=255)
                                                                                      GGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATTTTTCCGCTGTa                                                            >  2:1354415‑M1/1‑150 (MQ=255)
                                                                                         gatccaggcgggtaacgggggcctgagcgatgccgaccgtacctcgatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGa                                                                                                                            <  1:1196210‑M2/28‑1 (MQ=255)
                                                                                         gatccaggcgggtaacgggggcctgagcgatgccgaccgtacctcgatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGa                                                                                                                            >  2:1196210‑M2/56‑83 (MQ=255)
                                                                                          atccaggcgggtaacgggggcctgagcgatgccgaccgtacctcgatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGc                                                        <  2:1646742‑M2/96‑1 (MQ=255)
                                                                                           TCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCt                                                       <  2:669618‑M1/150‑1 (MQ=255)
                                                                                            CGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTc                                                      >  1:669618‑M1/1‑150 (MQ=255)
                                                                                                   cGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAg                                               <  2:1603958/150‑1 (MQ=255)
                                                                                                   cGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACGGTTCCAACAAg                                               <  2:713161/150‑1 (MQ=255)
                                                                                                     gTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGgc                                             >  1:1135024/1‑150 (MQ=255)
                                                                                                     gTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGgc                                             <  2:160780/150‑1 (MQ=255)
                                                                                                         gCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGATCCTGTTGTACTGCTCCAACAAGGCGctc                                         <  1:2245400/150‑1 (MQ=255)
                                                                                                          cATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCcaacaa                                                >  1:1838509/1‑142 (MQ=255)
                                                                                                          cATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCcaacaa                                                <  2:1838509/142‑1 (MQ=255)
                                                                                                          cATCGCTCAGGCCCCCGTTACACGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGctct                                        >  2:577687/1‑150 (MQ=255)
                                                                                                           aTCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCGTGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTg                                       >  1:913720/1‑150 (MQ=255)
                                                                                                            tCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGt                                      <  1:1228244/150‑1 (MQ=255)
                                                                                                              gCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAgg                                                                                                                                      <  1:994776/52‑1 (MQ=255)
                                                                                                                  aGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTtgctgc                                <  2:1850304/150‑1 (MQ=255)
                                                                                                                    gCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTgcag                              >  2:91730/1‑150 (MQ=255)
                                                                                                                         cGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCtt                                                                                    >  1:400054/1‑91 (MQ=255)
                                                                                                                         cGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCtt                                                                                    <  2:400054/91‑1 (MQ=255)
                                                                                                                               ccGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGa                                                                                                                            <  2:2112344/45‑1 (MQ=255)
                                                                                                                               ccGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGa                                                                                                                            >  1:2112344/1‑45 (MQ=255)
                                                                                                                               ccGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGa                                                                                  >  1:465253/1‑87 (MQ=255)
                                                                                                                               ccGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCcaacaa                                                >  2:681595/1‑121 (MQ=255)
                                                                                                                               ccGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGgc                   >  2:979559/1‑150 (MQ=255)
                                                                                                                                 gCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTAc                                                           >  2:1190359/1‑108 (MQ=255)
                                                                                                                                  ccTGGATCGCCAGCTGGCTGGCACGCTGCAGGGCATAGTTGAAGGGCGAGAGCACGCAATCTACCTGCAACCGGGCGTACGTTATCGTGGGCATGTTTTCGCTGTACTGCTCCAACAAGGCGATCTGGTGCTGCCGCTCTCGCAGgcgcg                >  2:934137/1‑150 (MQ=255)
                                                                                                                                  ccTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGct                                          >  2:1389308/1‑124 (MQ=255)
                                                                                                                                  ccTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGct                                          <  1:1389308/124‑1 (MQ=255)
                                                                                                                                     ggATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTa                                                            <  2:728837/103‑1 (MQ=255)
                                                                                                                                     ggATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTa                                                            >  1:728837/1‑103 (MQ=255)
                                                                                                                                      gatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCg                                                                 >  2:582667‑M2/11‑97 (MQ=255)
                                                                                                                                      gatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCg                                                                 <  1:582667‑M2/87‑1 (MQ=255)
                                                                                                                                      gatcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCCCTGTACTGCTCCAACCCGGTGCTCTGTTGATGCAGCACCCGCAGGCGCgcacc            <  2:28147‑M2/140‑3 (MQ=255)
                                                                                                                                        tcgccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAgcgc          <  1:979559‑M2/142‑1 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGt                                                                           >  2:1709441‑M2/9‑85 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTtgctg                                 >  2:93945‑M2/9‑127 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAgcgc          >  2:353306‑M2/9‑150 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAgcgc          >  1:1370840‑M2/9‑150 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATAGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTtgctgc                                >  1:2007679‑M2/9‑128 (MQ=255)
                                                                                                                                        gtgccagcGAGCTGGCACGCTGCAGGGCATCGTTGATAGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTtgctgc                                <  2:2007679‑M2/120‑1 (MQ=255)
                                                                                                                                          gccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTc                                                                   >  2:411818‑M2/7‑91 (MQ=255)
                                                                                                                                          gccagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTc                                                                   <  1:411818‑M2/85‑1 (MQ=255)
                                                                                                                                            cagcTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAGCGCCGAc      >  2:1804118‑M2/5‑150 (MQ=255)
                                                                                                                                               cTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAGCGCCGACCgg   >  1:483158‑M2/2‑150 (MQ=255)
                                                                                                                                                tCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAGCGCCGACCggg  >  2:2075991/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TTGGAGCCGCCAAACATGTAGTCGCCGTTGGCATCGCGCGAGTTCAGCAGACCGAACACGTTGGCTTCGATTTCCTTGATCTCGCTGGCGATCGAGGTACGGTCGGCATCGCTCAGGCCCCCGTTACCCGCCTGGATCGCCAGCTCGCTGGCACGCTGCAGGGCATCGTTGATGGTCGATAGCACGCTCTCTTCCTGCAGCAGGGCGTTCTTGACCGTGGTCATGTTTCCGCTGTACTGCTCCAACAAGGCGCTCTGTTGCTGCAGCACCAGCAGGCGCGCAGCGCCGACCGGG  >  NC_002947/4973174‑4973467

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