Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,377,749 T→G 8.8% S175R (AGC→CGC)  grpE ← heat shock protein GrpE

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,377,7490TG 8.8% 408.4 / 25.9 148S175R (AGC→CGC) grpEheat shock protein GrpE
Reads supporting (aligned to +/- strand):  ref base T (79/56);  new base G (6/7);  total (85/63)
Fisher's exact test for biased strand distribution p-value = 3.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 149 total aligned reads displayed.

TTTACCCATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGAACGGCT  >  NC_002947/5377604‑5377898
                                                                                                                                                 |                                                                                                                                                     
   aCCCATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGcc                                                                                                                                                >  1:717183/1‑150 (MQ=255)
      cATGTTAACTCTCCCGAAACTTGCATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTg                                                                                                                                             <  1:2308283/150‑1 (MQ=255)
      cATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGAGCGCGgcgcaacc                                                                                                                                             <  1:2118634/150‑5 (MQ=255)
       aTGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGc                                                                                                                                            <  2:1340174/150‑1 (MQ=255)
       aTGTTAACTCTCCCGAAACCTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGc                                                                                                                                            >  1:188361/1‑150 (MQ=255)
               tctcCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGaccacc                                                                                                                                    >  2:2281920/1‑150 (MQ=255)
               tctcCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGaccacc                                                                                                                                    >  1:192806/1‑150 (MQ=255)
               tctcCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGaccacc                                                                                                                                    >  1:989615/1‑150 (MQ=255)
                   ccGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGacccccaccc                                                                                                                                >  2:2067119/1‑149 (MQ=255)
                   ccGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGAccaccacca                                                                                                                                >  1:254619/1‑150 (MQ=255)
                    cGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCAt                                                                                                                               >  1:176259/1‑150 (MQ=255)
                      aaaCTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATgg                                                                                                                             <  2:157468/150‑1 (MQ=255)
                          ttGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCggg                                                                                                                         >  1:1420308/1‑150 (MQ=255)
                           tGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCgggg                                                                                                                        <  2:1273485/150‑1 (MQ=255)
                               tttGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGGAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCgca                                                                                                                    >  1:591346/1‑150 (MQ=255)
                                   gTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGgcg                                                                                                                      <  2:1775150/144‑1 (MQ=255)
                                   gTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGgcg                                                                                                                      >  1:1775150/1‑144 (MQ=255)
                                    tAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAgg                                                                                                               >  2:883471/1‑150 (MQ=255)
                                      gcagcGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGAGCGCGgcgcaac                                                                                                                                              >  1:1329699/1‑114 (MQ=255)
                                      gcagcGACTACTTTGGCCAACTGAGGTACTACTTGAAGACTTTACTCCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAggcg                                                                                                             <  1:539460/150‑1 (MQ=255)
                                       cagcGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAggcgg                                                                                                            >  2:911689/1‑150 (MQ=255)
                                        agcGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAggcggc                                                                                                           <  2:1780767/150‑1 (MQ=255)
                                         gcGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGACTTGCTGACCACCACCATGCCGGGGCCCACCAggcggct                                                                                                          >  2:465657/1‑149 (MQ=255)
                                          cGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCg                                                                                                         <  2:989615/150‑1 (MQ=255)
                                           gACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCg                                                                                                                                                     >  2:1047521/1‑105 (MQ=255)
                                           gACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCg                                                                                                                                                     <  1:1047521/105‑1 (MQ=255)
                                                ctTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTcagc                                                                                                   <  2:254619/150‑1 (MQ=255)
                                                ctTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTcagc                                                                                                   >  1:1935800/1‑150 (MQ=255)
                                                 tttGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGcc                                                                                                          >  2:1564277/1‑142 (MQ=255)
                                                 tttGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTcagca                                                                                                  <  1:2281920/150‑1 (MQ=255)
                                                     gCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGAAGGCGGCCGTTcagccgcta                                                                                              >  2:262545/1‑145 (MQ=255)
                                                      ccAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCgca                                                                                                                    <  1:2345164/127‑1 (MQ=255)
                                                      ccAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCgca                                                                                                                    >  2:2345164/1‑127 (MQ=255)
                                                      ccAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTcagc                                                                                                   <  2:547941/144‑1 (MQ=255)
                                                        aaCTGCGGTAATAATTGAAGGCTTGACACCTAGATGGGGATGCCCCGGAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAAcc                                                                                           >  1:1632596/1‑150 (MQ=255)
                                                          cTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCt                                                                                         >  2:1472757/1‑150 (MQ=255)
                                                            gCGGTAATACTTGAACGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTGAACCCAGGTGAGTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTc                                                                                       >  2:695633/1‑150 (MQ=255)
                                                             cGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCt                                                                                      >  2:1972791/1‑150 (MQ=255)
                                                              ggTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCGCCGATCGAGGGTTGCACCACGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTg                                                                                     <  1:781081/150‑1 (MQ=255)
                                                              ggTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGAGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTg                                                                                     <  1:1059339/150‑1 (MQ=255)
                                                                tAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTg                                                                                     <  1:694579/148‑1 (MQ=255)
                                                                 aaTACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGaa                                                                                  <  2:710533/150‑1 (MQ=255)
                                                                   tACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGacca                                                                                                                                      >  1:1579832/1‑96 (MQ=255)
                                                                    aCTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTa                                                                                              <  2:608568/135‑1 (MQ=255)
                                                                    aCTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTa                                                                                              >  1:608568/1‑135 (MQ=255)
                                                                    aCTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAacac                                                                               <  2:1975097/150‑1 (MQ=255)
                                                                      ttGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCgca                                                                                                                    <  2:1727385/111‑1 (MQ=255)
                                                                      ttGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCgca                                                                                                                    >  1:1727385/1‑111 (MQ=255)
                                                                      ttGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGt                                                                             >  2:2276971/1‑150 (MQ=255)
                                                                      ttGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGt                                                                             <  2:1935800/150‑1 (MQ=255)
                                                                       tGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGtt                                                                            <  2:192806/150‑1 (MQ=255)
                                                                         aaGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGa                                                                                                                                         >  1:39501/1‑87 (MQ=255)
                                                                         aaGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGa                                                                                                                                         <  2:39501/87‑1 (MQ=255)
                                                                                gACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCt                                                                   >  1:1238226/1‑150 (MQ=255)
                                                                                   acCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGtt                                                                <  1:911689/150‑1 (MQ=255)
                                                                                     cTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCACCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCt                                                                                      >  2:1627041/1‑126 (MQ=255)
                                                                                      tAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGac                                                                                                                                        >  1:1514416/1‑75 (MQ=255)
                                                                                       aGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATg                                                                                                                              <  2:1004219/84‑1 (MQ=255)
                                                                                       aGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATg                                                                                                                              >  1:1004219/1‑84 (MQ=255)
                                                                                         aTGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTc                                                          <  2:2383708/150‑1 (MQ=255)
                                                                                              gATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGt                                                                                               <  1:1658466/108‑1 (MQ=255)
                                                                                                  ccccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGc                                                 >  1:2318222/1‑150 (MQ=255)
                                                                                                   cccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCg                                                >  1:386069/1‑150 (MQ=255)
                                                                                                    ccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCgc                                               >  2:793649/1‑149 (MQ=255)
                                                                                                    ccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCg                                                >  1:2370002/1‑149 (MQ=255)
                                                                                                    ccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCg                                                <  2:2370002/149‑1 (MQ=255)
                                                                                                    ccGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGt                                               >  2:213825/1‑150 (MQ=255)
                                                                                                     cGAGGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGtt                                              >  1:1778637/1‑150 (MQ=255)
                                                                                                     cGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGtt                                              >  1:647542/1‑150 (MQ=255)
                                                                                                          gATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCt                                         >  1:2139271/1‑150 (MQ=255)
                                                                                                          gATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTACCt                                         >  1:846078/1‑150 (MQ=255)
                                                                                                            tttCAAGCCTTTTCATTGATCGAGATTTGCGCCGGGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCTGCGTTTTCCTGc                                       >  2:1355031/1‑150 (MQ=255)
                                                                                                                   ccTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGt                                                                 <  2:1846195/117‑1 (MQ=255)
                                                                                                                   ccTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGt                                                                 >  1:1846195/1‑117 (MQ=255)
                                                                                                                    cTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCt                                                           <  1:443486/122‑1 (MQ=255)
                                                                                                                    cTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCt                                                           >  2:443486/1‑122 (MQ=255)
                                                                                                                       ttCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGcatggccatggc                            <  2:1708847/150‑1 (MQ=255)
                                                                                                                          aTTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGtt                                                                <  2:2096750/111‑1 (MQ=255)
                                                                                                                          aTTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGtt                                                                >  1:2096750/1‑111 (MQ=255)
                                                                                                                          aTTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCtg                         <  1:883471/150‑1 (MQ=255)
                                                                                                                          aTTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCtg                         >  1:343187/1‑150 (MQ=255)
                                                                                                                          aTTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCtg                         >  2:1668820/1‑150 (MQ=255)
                                                                                                                           ttGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCtgg                        >  2:1490744/1‑150 (MQ=255)
                                                                                                                              aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGt                                                                             <  2:1249740/94‑1 (MQ=255)
                                                                                                                              aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGcat                                     >  1:185625/1‑134 (MQ=255)
                                                                                                                              aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGcat                                     <  2:185625/134‑1 (MQ=255)
                                                                                                                              aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCatggccatggcc                           <  2:1594815/144‑1 (MQ=255)
                                                                                                                              aTCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCatggccatggcc                           >  1:1594815/1‑144 (MQ=255)
                                                                                                                               tCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTggtggt                    >  1:437554/1‑150 (MQ=255)
                                                                                                                                     gTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCt                                                                                      <  1:511837/78‑1 (MQ=255)
                                                                                                                                     gTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCt                                                                                      >  2:511837/1‑78 (MQ=255)
                                                                                                                                     aaggcgcCGAGCGCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTc                                           <  2:433088/118‑1 (MQ=255)
                                                                                                                                     aaggcgcCGAGCGCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAg              >  1:1885679/4‑150 (MQ=255)
                                                                                                                                     aaggcgcCGAGCGCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAg              >  1:1084736/4‑150 (MQ=255)
                                                                                                                                       tGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGACATGGCCTGGTGGTGCTcccag            <  2:644419/150‑5 (MQ=255)
                                                                                                                                           ccgcGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGa        >  2:2137179/1‑150 (MQ=255)
                                                                                                                                            cgcgCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGaa       <  2:2342577/150‑1 (MQ=255)
                                                                                                                                             gcgcTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGAAc      >  2:1667768/1‑150 (MQ=255)
                                                                                                                                                cTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTggtggt                    <  2:1081024/133‑1 (MQ=255)
                                                                                                                                                 tCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACGCGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGAACGGCt  <  1:1472757/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TTTACCCATGTTAACTCTCCCGAAACTTGAATTTGGTAGCAGCGACTACTTTGGCCAACTGCGGTAATACTTGAAGGCTTGACACCTAGATGGGGATGCCCCGAAGGATTTCAAGCCTTTTCATTGATCGAGGGTTGCGCCGCGCTCGGCGCCTTGCTGACCACCACCATGGCGGGGCGCAGCAGGCGGCCGTTCAGCAGGTAACCCTTCTGGAACACGTTCAACACGCTGTTTGGCTCGACCTCGGCGTTTTCCTGCATGGCCATGGCCTGGTGGTGCTCAGGGTTGAACGGCT  >  NC_002947/5377604‑5377898

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