Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 383,377 C→T 6.3% intergenic (+34/‑149) PP_0319 → / → PP_0320 hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947383,3770CT 6.3% 425.1 / 12.9 145intergenic (+34/‑149)PP_0319/PP_0320hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (62/71);  new base T (3/6);  total (66/79)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

Only 100 of 146 total aligned reads displayed.

TGTGCTCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCA  >  NC_002947/383231‑383526
                                                                                                                                                  |                                                                                                                                                     
tgtgCTCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACC                                                                                                                                                    >  1:4725950‑M1/1‑150 (MQ=255)
   gCTCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGC                                                                                                                                                 >  1:1025653‑M1/1‑150 (MQ=255)
   gCTCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGC                                                                                                                                                 >  1:3644051‑M1/1‑150 (MQ=255)
     tCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGA                                                                                                                                               >  2:2744042‑M1/1‑150 (MQ=255)
     tCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGA                                                                                                                                               >  2:4335175‑M1/1‑150 (MQ=255)
     tCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTGAACCCGCtc                                                                                                                                               <  2:1377621‑M1/150‑3 (MQ=255)
        aaGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATG                                                                                                                                            >  2:3741339‑M1/1‑150 (MQ=255)
        aaGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATG                                                                                                                                            >  2:4860650‑M1/1‑150 (MQ=255)
         aGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGC                                                                                                                                           >  2:1883659‑M1/1‑150 (MQ=255)
         aGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGC                                                                                                                                           <  2:3644051‑M1/150‑1 (MQ=255)
         aGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGC                                                                                                                                           <  2:1025653‑M1/150‑1 (MQ=255)
         aGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGC                                                                                                                                           <  1:2290851‑M1/150‑1 (MQ=255)
         aGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGC                                                                                                                                           <  1:693697‑M1/150‑1 (MQ=255)
              gTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCTCCCCCGAATGCGTCAg                                                                                                                                      >  1:2281060/1‑150 (MQ=255)
              gTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGgtgt                                                                                                                                      <  1:3026170‑M1/150‑5 (MQ=255)
              gTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAg                                                                                                                                      >  1:2587493/1‑150 (MQ=255)
               ttGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGt                                                                                                                                     <  2:2417582/150‑1 (MQ=255)
                 gAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTgg                                                                                                                                   <  2:1507730/150‑1 (MQ=255)
                 gAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCATACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTgg                                                                                                                                   <  1:12508/150‑1 (MQ=255)
                  aGCTGCAGAAGGAGCACCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGc                                                                                                                                  <  1:4272781/150‑1 (MQ=255)
                   gCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCt                                                                                                                                 <  2:3978096/150‑1 (MQ=255)
                    cTGCAGAAGGAGCATCAGCAGCCCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCtt                                                                                                                                <  2:4677016/150‑1 (MQ=255)
                    cTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCtt                                                                                                                                <  1:1883659/150‑1 (MQ=255)
                         gAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAaca                                                                                                                           <  1:4040535/150‑1 (MQ=255)
                          aaGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAacac                                                                                                                          >  2:2680271/1‑150 (MQ=255)
                          aaGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGGCTTCAacac                                                                                                                          >  2:4702539/1‑150 (MQ=255)
                             gAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcgc                                                                                                                       <  1:1396943‑M2/150‑4 (MQ=255)
                               gCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcat                                                                                                                     <  2:510301‑M1/150‑22 (MQ=255)
                                 aTCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcattc                                                                                                                   >  1:2693361‑M1/1‑127 (MQ=255)
                                  tCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAACGCGTCAGTGGCTTCAACACcgcattcg                                                                                                                  <  2:2945267‑M2/150‑9 (MQ=255)
                                   cagcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgc                                                                                                                 >  1:6749‑M2/1‑141 (MQ=255)
                                    agcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcattcgcg                                                                                                                <  2:2693361‑M1/150‑27 (MQ=255)
                                     gcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcattcgcgg                                                                                                               <  1:2633779‑M1/150‑28 (MQ=255)
                                     gcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgcgg                                                                                                               >  2:655481‑M2/1‑139 (MQ=255)
                                     gcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgcgg                                                                                                               >  2:1327756‑M2/1‑139 (MQ=255)
                                     gcagACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgcgg                                                                                                               >  1:1895425‑M2/1‑139 (MQ=255)
                                           ccTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAc                                                                                                         <  2:2190975/150‑1 (MQ=255)
                                           ccTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAc                                                                                                         <  1:2776012/150‑1 (MQ=255)
                                           ccTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTAAAAACCCCATTAGCGGGTAAAc                                                                                                         >  1:930954/1‑150 (MQ=255)
                                             ttGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAccc                                                                                                       >  1:4764128/1‑150 (MQ=255)
                                                 gcgAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTc                                                                                                   <  1:5336298/150‑1 (MQ=255)
                                                 gcgAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCAACCGCGAATTCGTAACTGGCTTCAACACCGCAGTCGCGGGTAAAACCGCTc                                                                                                   >  1:5190221/1‑150 (MQ=255)
                                                   gAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAAAGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCt                                                                                                 >  1:3661373/1‑150 (MQ=255)
                                                     aCTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcattcgc                                                                                                                 >  1:2981200‑M1/1‑107 (MQ=255)
                                                     aCTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTTACCCGCGAATGCGgtgttgaagccactgacgcattcgc                                                                                                                 <  2:2981200‑M1/132‑26 (MQ=255)
                                                       tCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGgtgttgaagccactgacgcattcgcgg                                                                                                               <  1:5031657‑M1/132‑28 (MQ=255)
                                                       tCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGgtgttgaagccactgacgcattcgcgg                                                                                                               >  2:5031657‑M1/1‑105 (MQ=255)
                                                       tCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAg                                                                                             >  2:5125687/1‑150 (MQ=255)
                                                        cAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTAGCGGGTAAACCCGCTCCTACAgg                                                                                            <  1:1520631/150‑1 (MQ=255)
                                                            cgcAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGc                                                                                        >  1:3159345/1‑150 (MQ=255)
                                                               aGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCTGTAGGAGCGGGTTTACCCG                                                                                                                                                  >  2:752472‑M1/1‑89 (MQ=255)
                                                               aGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATc                                                                                     >  2:5434722/1‑150 (MQ=255)
                                                                gATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCt                                                                                    <  2:4585744/150‑1 (MQ=255)
                                                                 aTCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGgtgttgaagccactgacgcattcgcgggtgaa                                                                                                          >  1:2188483‑M1/1‑95 (MQ=255)
                                                                  tCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGa                                                                                  <  2:1897753/150‑1 (MQ=255)
                                                                     aCGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAAAACCGCATTAGCGGGTAAAACAGAAACTACAGGGTGCATCTGATTg                                                                               >  2:1871191/1‑150 (MQ=255)
                                                                      cGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGc                                                                              <  2:1596755/150‑1 (MQ=255)
                                                                        ccACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCtt                                                                            <  2:2424298/150‑1 (MQ=255)
                                                                        ccACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCtt                                                                            <  1:3761806/150‑1 (MQ=255)
                                                                        ccACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCtt                                                                            <  2:3973546/150‑1 (MQ=255)
                                                                         cACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTc                                                                           >  2:5310577/1‑150 (MQ=255)
                                                                          aCGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCa                                                                          <  2:1706004/150‑1 (MQ=255)
                                                                          aCGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACACGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACGCGCTCCTACAGGGTTCATCTGATTGCTTCa                                                                          >  2:2753963/1‑150 (MQ=255)
                                                                             gTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGAt                                                                       >  1:4001224/1‑150 (MQ=255)
                                                                                     cATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACcg                                                                                                                        <  2:408542‑M2/93‑3 (MQ=255)
                                                                                          tGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACt                                                          <  1:2488110/150‑1 (MQ=255)
                                                                                            cGGACCATGGCTGAGCAATACGCCCCGTACACCTTCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGa                                                       <  2:122347/150‑1 (MQ=255)
                                                                                            cGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTgg                                                        <  2:1562220/150‑1 (MQ=255)
                                                                                                ccATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATAc                                                    <  1:3301554/150‑1 (MQ=255)
                                                                                                    cGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCaa                                                <  1:3222293/150‑1 (MQ=255)
                                                                                                    cGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGGCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCaa                                                >  2:4715770/1‑150 (MQ=255)
                                                                                                      cTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAgt                                              <  1:2074945/150‑1 (MQ=255)
                                                                                                      cTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCCCTCCTACAGGGTGCAACTGATTGCTTCAGATCACGGGCTTAACTGGATACGCAAgt                                              >  1:4631330/1‑150 (MQ=255)
                                                                                                       tGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAcc                                                                                                        >  1:3758204/1‑91 (MQ=255)
                                                                                                         aGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGta                                           <  2:5210763/150‑1 (MQ=255)
                                                                                                           cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGtata                                         >  2:2053008/1‑150 (MQ=255)
                                                                                                           cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGtata                                         <  2:4725950/150‑1 (MQ=255)
                                                                                                           cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGATAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGAt                                                                       >  2:903832/1‑120 (MQ=255)
                                                                                                           cAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGATAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGAt                                                                       <  1:903832/120‑1 (MQ=255)
                                                                                                            aaTAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATac                                        <  1:1028531/150‑1 (MQ=255)
                                                                                                                    ccGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACt                                                                   >  2:4435447/1‑115 (MQ=255)
                                                                                                                    ccGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACt                                                                   <  1:4435447/115‑1 (MQ=255)
                                                                                                                    ccGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCg                                <  2:889381/150‑1 (MQ=255)
                                                                                                                       tACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCaa                             >  2:3723864/1‑150 (MQ=255)
                                                                                                                          acCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAc                                                                                                         >  1:1482694/1‑71 (MQ=255)
                                                                                                                          acCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAAc                                                                                                         <  2:1482694/71‑1 (MQ=255)
                                                                                                                                cccccGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACaataat                    <  1:2134528/150‑1 (MQ=255)
                                                                                                                                 ccccGTAGGAGCGGGTTCACCCCCGAATGCGTCAGTGGCTTCAACACCGCCGTCGCGGGGAAACCCGCTCCAACAGGGTGCATCTGATTTCTTCAGATCACTGGATTAACTGGATACGCACCTGGTTACACTTTCCGCAAAACAATAAtg                   >  1:4235215/1‑150 (MQ=255)
                                                                                                                                  cccGTAGGAGCGGGTTCAACCGCGAATGCGTCAGTGGCTTCAAAAACGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCATCAGATCACTGGATTAACTGGATAAGAAAGTGTATACACTTTCCGCAAAACAATAAtgt                  >  1:4985120/1‑150 (MQ=38)
                                                                                                                                    cGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGc                <  1:1453092/150‑1 (MQ=255)
                                                                                                                                     gTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCg               <  2:1418164/150‑1 (MQ=255)
                                                                                                                                     gTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCg               <  2:1418213/150‑1 (MQ=255)
                                                                                                                                     gTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCg               <  2:1577631/150‑1 (MQ=255)
                                                                                                                                      tAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCgg              <  1:2053008/150‑1 (MQ=255)
                                                                                                                                       attcgCGGGTGAACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCggg             >  2:2073091/5‑150 (MQ=255)
                                                                                                                                       aGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCggg             <  2:1143621/150‑1 (MQ=255)
                                                                                                                                         gAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGa           <  1:2460610/150‑1 (MQ=255)
                                                                                                                                            cGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCgggaggtg        >  1:4798640/1‑145 (MQ=255)
                                                                                                                                                ttCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgcgggtgaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgat                           >  1:2188483‑M2/1‑32 (MQ=255)
                                                                                                                                                ttCACCCGCGAATGCGTCAGTGGCTTCAACACcgcattcgcgggtgaacccgctcctacggggggacggtgtacggggcttattgctcagcgatggtccgcacgatggcatccaccgtggcgtcgat                           <  2:2188483‑M2/127‑96 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
TGTGCTCCAAGCAGGTTGAGCTGCAGAAGGAGCATCAGCAGACCCTTGCGCGAACTCAAGCGCAGATCGACGCCACGGTGGATGCCATCGTGCGGACCATCGCTGAGCAATAAGCCCCGTACACCGTCCCCCCGTAGGAGCGGGTTCACCCGCGAATGCGTCAGTGGCTTCAACACCGCATTCGCGGGTAAACCCGCTCCTACAGGGTGCATCTGATTGCTTCAGATCACTGGATTAACTGGATACGCAAGTGTATACACTTTCCGCAAAACAATAATGTGCGGGGAGTAGGGGCA  >  NC_002947/383231‑383526

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