Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,198,294 G→A 6.6% intergenic (‑665/+58) PP_1940 ← / ← PP_1941 methyl‑accepting chemotaxis transducer/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,198,2940GA 6.6% 465.6 / 15.0 152intergenic (‑665/+58)PP_1940/PP_1941methyl‑accepting chemotaxis transducer/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (66/76);  new base A (5/5);  total (71/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

GTGACCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCACATCATTAG  >  NC_002947/2198147‑2198441
                                                                                                                                                   |                                                                                                                                                   
  gACCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCgctgc                                                                                                                                                 <  1:3790956/150‑1 (MQ=255)
  gACCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCgctgc                                                                                                                                                 <  1:3790795/150‑1 (MQ=255)
   aCCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCgctgct                                                                                                                                                >  2:3942096/1‑150 (MQ=255)
   aCCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCgctgct                                                                                                                                                >  2:3942115/1‑150 (MQ=255)
     cccTGGTCGGTGGGGACCTGAATCATGCAAGGATGGAGGCCCATCTGCTATAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCtgctgc                                                                                                                                              <  2:2739556/150‑1 (MQ=255)
     cccAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCtgctgc                                                                                                                                              >  2:3808077/1‑150 (MQ=255)
     cccAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCtgctgc                                                                                                                                              >  2:1187570/1‑150 (MQ=255)
     cccAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCtgctgc                                                                                                                                              <  2:1424204/150‑1 (MQ=255)
       cAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCga                                                                                                                                            <  1:2872251/150‑1 (MQ=255)
       cAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCga                                                                                                                                            <  2:3165146/150‑1 (MQ=255)
        aGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCAGCGCTGCggg                                                                                                                                           <  1:1763956/150‑3 (MQ=255)
         ggTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCTCGCAGCAGCGCTGCggga                                                                                                                                          <  1:4586741/150‑1 (MQ=255)
         ggTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCgaga                                                                                                                                          <  1:5393861/150‑1 (MQ=255)
             ggTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTaa                                                                                                                                      >  1:3566123/1‑150 (MQ=255)
             ggTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTaa                                                                                                                                      <  1:307999/150‑1 (MQ=255)
             ggTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTaa                                                                                                                                      <  1:307784/150‑1 (MQ=255)
              gTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAAc                                                                                                                                     <  1:1715851/150‑1 (MQ=255)
                  ggAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCTCGCAGCAGCGCTGCGGGAtttcattgc                                                                                                                                 <  1:111906/150‑8 (MQ=255)
                   gatcTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGTCTTCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTTcc                                                                                                                                <  1:3591052/147‑1 (MQ=255)
                         aaTCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATAAACCAACCGAgatc                                                                                                                          >  2:3591052/1‑147 (MQ=255)
                             atcCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAg                                                                                                                      <  1:108001/150‑1 (MQ=255)
                              tcCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGAtttcatt                                                                                                                                   >  1:4283125/1‑131 (MQ=255)
                               ccAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGc                                                                                                                    >  1:942711/1‑150 (MQ=255)
                                 aaGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGccc                                                                                                                  >  1:4333431/1‑150 (MQ=255)
                                 aaGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGccc                                                                                                                  >  1:1550076/1‑150 (MQ=255)
                                     aTGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGtt                                                                                                              <  1:2610058/150‑1 (MQ=255)
                                       tgacgcacATCTTCTAAAGCTCATTGATCTATACAGCAGTGATCTGAACGACTTTACTATCAGACATCCTGTTCGACTTGATTAGTCAATTCAATGAAATACCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTgg                                                                                                            <  2:3408747/143‑1 (MQ=255)
                                             ccATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAg                                                                                                      >  2:583002/1‑150 (MQ=255)
                                              cATCTGCTAAAGCGCCTGGAACGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAgg                                                                                                     <  1:1121752/150‑1 (MQ=255)
                                              cATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAgg                                                                                                     <  1:1241262/150‑1 (MQ=255)
                                                tCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGc                                                                                                   <  1:5324914/150‑1 (MQ=255)
                                                 cTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCa                                                                                                  >  2:4366841/1‑150 (MQ=255)
                                                  tGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCaa                                                                                                 <  2:2098122/150‑1 (MQ=255)
                                                  tGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCaa                                                                                                 >  2:2956835/1‑150 (MQ=255)
                                                   gCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAAt                                                                                                <  2:4246545/150‑1 (MQ=255)
                                                   gCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAAt                                                                                                <  1:4753374/150‑1 (MQ=255)
                                                     tAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATga                                                                                              >  1:2792527/1‑150 (MQ=255)
                                                     tAAAGCGCATGGATCGCGTCAGCAGTGATCGGAACGCCTTTACTCTCAGACCTCCTGTGCGCCTTGAGTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATga                                                                                              <  2:849148/150‑1 (MQ=255)
                                                      aaaGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATgag                                                                                             >  1:4237921/1‑150 (MQ=255)
                                                      aaaGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATgag                                                                                             >  1:212319/1‑150 (MQ=255)
                                                           gcctGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTaa                                                                                        <  1:5087014/147‑1 (MQ=255)
                                                           gcATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTaa                                                                                        >  1:5006072/1‑150 (MQ=255)
                                                           gcATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTaa                                                                                        <  2:4539932/150‑1 (MQ=255)
                                                             aTGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGa                                                                                      >  2:1077057/1‑150 (MQ=255)
                                                                    gCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActct                                                                               <  2:3731036/150‑1 (MQ=255)
                                                                     cTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctc                                                                              <  1:3490949/150‑1 (MQ=255)
                                                                     cTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctc                                                                              <  2:1651319/150‑1 (MQ=255)
                                                                      ttCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctct                                                                             <  2:500468/150‑1 (MQ=255)
                                                                       tCAGCAGTGATCGGAACGCCTTTACTATCAGACTTCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCtt                                                                            <  1:5193891/150‑1 (MQ=255)
                                                                       tCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCtt                                                                            >  1:5355092/1‑150 (MQ=255)
                                                                          gcagTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGAc                                                                         <  1:5293645/150‑1 (MQ=255)
                                                                           cagTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGAcc                                                                        >  1:3985028/1‑150 (MQ=255)
                                                                           cagTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGAcc                                                                        >  1:3985055/1‑150 (MQ=255)
                                                                            agTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCa                                                                       >  1:1274152/1‑150 (MQ=255)
                                                                            agTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCa                                                                       >  2:5267370/1‑150 (MQ=255)
                                                                               gATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTg                                                                    <  1:2044751/150‑1 (MQ=255)
                                                                                aTCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGc                                                                   >  1:2514654/1‑150 (MQ=255)
                                                                                   ggAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCa                                                                >  2:1423730/1‑150 (MQ=255)
                                                                                   ggAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCa                                                                <  2:3774671/150‑1 (MQ=255)
                                                                                        gCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCaa                                                           <  2:3593324/150‑1 (MQ=255)
                                                                                           tttACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGa                                                        >  1:644295/1‑150 (MQ=255)
                                                                                              aCTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGAttt                                                     <  1:663785/150‑1 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctct                                                                             <  1:4079646/122‑1 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctct                                                                             <  1:4079618/122‑1 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctct                                                                             >  2:4079618/1‑122 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAActctct                                                                             >  2:4079646/1‑122 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCaaa                                                 >  1:5228343/1‑150 (MQ=255)
                                                                                                  tCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCaaa                                                 <  1:3515508/150‑1 (MQ=255)
                                                                                                   cAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCaaaa                                                <  1:5103279/150‑1 (MQ=255)
                                                                                                    aGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAAc                                               >  2:2011701/1‑150 (MQ=255)
                                                                                                    aGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAAc                                               >  1:5020458/1‑150 (MQ=255)
                                                                                                      aCATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCt                                             >  2:1926654/1‑150 (MQ=255)
                                                                                                       cATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTg                                            <  2:2296301/150‑1 (MQ=255)
                                                                                                          ccTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTc                                                                                                                         >  1:69383/1‑70 (MQ=255)
                                                                                                          ccTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTc                                                                                                                         <  2:69383/70‑1 (MQ=255)
                                                                                                           cTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAAc                                                                                                                                     <  2:5296149/57‑1 (MQ=255)
                                                                                                           cTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAAc                                                                                                                                     >  1:5296149/1‑57 (MQ=255)
                                                                                                            tgtgCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTcgc                                       >  1:4869571/1‑150 (MQ=255)
                                                                                                              tgCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGa                                     <  1:370064/150‑1 (MQ=255)
                                                                                                                cgcCTTGATTACGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCCAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGt                                   <  2:1457832/150‑1 (MQ=255)
                                                                                                                  ccTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTaa                                 >  2:4677582/1‑150 (MQ=255)
                                                                                                                     tGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAg                              >  2:3196584/1‑150 (MQ=255)
                                                                                                                        ctaatcaaggcaaTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAgggg                           >  2:3449422/6‑150 (MQ=255)
                                                                                                                         tAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGc                          <  1:3283416/150‑1 (MQ=255)
                                                                                                                         tAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGc                          <  1:4855553/150‑1 (MQ=255)
                                                                                                                            gcaaggcaaTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGa                       <  1:2541510/150‑1 (MQ=255)
                                                                                                                                ggcaaTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCg                   >  1:703014/1‑150 (MQ=255)
                                                                                                                                 gcaaTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGc                  >  1:1214303/1‑150 (MQ=255)
                                                                                                                                   aaTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCt                <  2:3238937/150‑1 (MQ=255)
                                                                                                                                    aTGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTa               <  2:2193542/150‑1 (MQ=255)
                                                                                                                                     gttaaTCTCGCAGCAGCGCTGCGGGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAg              >  2:4451620/4‑150 (MQ=255)
                                                                                                                                     gttaaTCTCGCAGCAGCGCTGCGGGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAg              >  2:223499/4‑150 (MQ=255)
                                                                                                                                      ttaaTCTCGCAGCAGCGCTGCGGGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGc             >  2:4476542/3‑150 (MQ=255)
                                                                                                                                         aTCTCGCAGCAGCGCTGCGGGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCcac          >  1:2173274/1‑150 (MQ=255)
                                                                                                                                          tCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCcaca         <  2:3323835/150‑1 (MQ=255)
                                                                                                                                          tCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCcaca         >  1:5218650/1‑150 (MQ=255)
                                                                                                                                              gCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAATCCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCAcatcat     >  2:2605661/1‑150 (MQ=255)
                                                                                                                                              gCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCAcatcat     <  1:4275241/150‑1 (MQ=255)
                                                                                                                                                aGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCACATCATTa   <  2:4333431/150‑1 (MQ=255)
                                                                                                                                                 gcgcTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCACATCATTAg  <  1:2011701/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
GTGACCCCAGGTCGGTGGGGAACTGAATCATCCAAGGATGGAGGCCCATCTGCTAAAGCGCATGGATCGCTTCAGCAGTGATCGGAACGCCTTTACTATCAGACATCCTGTGCGCCTTGATTAGGCAAGGCAATGAAATCCCGCAGCGCTGCTGCGAGATTAACCTACCGAGCGTCCAGGCCCCGTTGGCACTAGGGCAATGAGACTAAGAAAACTCTCTTGACCACTGCTCAGCCAACGATTTCAAAACCTGCTCGCGAGTAATAGGGGCCGATCCGCCTAGCCACATCATTAG  >  NC_002947/2198147‑2198441

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