Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,015,057 C→G 7.8% R396P (CGC→CCC)  PP_0875 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,015,0570CG 7.8% 282.3 / 12.6 91R396P (CGC→CCC) PP_0875hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (45/38);  new base G (5/2);  total (51/40)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

GAAGCTCTGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGTTT  >  NC_002947/1014912‑1015205
                                                                                                                                                 |                                                                                                                                                    
gAAGCTCTGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGccc                                                                                                                                                  >  2:1382957/1‑150 (MQ=255)
gAAGCTCTGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGccc                                                                                                                                                  <  2:299437/150‑1 (MQ=255)
   gCTCTGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGTTGCGcccacc                                                                                                                                               >  2:2173835/1‑150 (MQ=255)
    ctctGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGcccaccc                                                                                                                                              <  1:1325733/150‑1 (MQ=255)
           tgatCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTgg                                                                                                                                       <  2:1903780/150‑1 (MQ=255)
           tgatCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCACCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTgg                                                                                                                                       <  1:475500/150‑1 (MQ=255)
             atCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTa                                                                                                                                     >  1:1079254/1‑150 (MQ=255)
              tCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGAGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAg                                                                                                                                    <  2:1363973/150‑1 (MQ=255)
                      tAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACg                                                                                                                            >  2:2028261/1‑150 (MQ=255)
                        gCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGcc                                                                                                                          <  1:2105611/150‑1 (MQ=255)
                         cccAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCCGGTTGGTAGCCCCgccgccc                                                                                                                         >  2:748766/1‑149 (MQ=255)
                          ccagccagCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCtt                                                                                                                        <  2:1247886/150‑1 (MQ=255)
                          ccagccagCGGAAGGGCACGGGAGGCCAGCATGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGATGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGCGAGTGCCCACCCCGGTTGGCAGCCCCGACCCCtt                                                                                                                        >  2:814025/1‑150 (MQ=255)
                                       gggCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGCGACGCTGGAGCTTCAGCGCCGCCAGGGTGCGCCCCCCCAGGTTTGTAGCCCCGACGCCTTCGCCGCCGTAgcc                                                                                                           >  1:420219/1‑150 (MQ=255)
                                         gCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAgccgc                                                                                                         <  1:1858910/150‑1 (MQ=255)
                                           tCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTg                                                                                                                                        >  2:1300998/1‑117 (MQ=255)
                                           tCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTg                                                                                                                                        <  1:1300998/117‑1 (MQ=255)
                                           tCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGccgccg                                                                                                       >  1:1864820/1‑150 (MQ=255)
                                             gggAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGc                                                                                                     >  2:1825041/1‑150 (MQ=255)
                                                  gCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGTGCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAggg                                                                                                <  1:783028/150‑1 (MQ=255)
                                                  gCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAgcc                                                                                                           >  2:2299780/1‑139 (MQ=255)
                                                  gCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAgcc                                                                                                           <  1:2299780/139‑1 (MQ=255)
                                                              ggCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCt                                                                                    >  1:677787/1‑150 (MQ=255)
                                                                    aTAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGt                                                                              <  1:1825041/150‑1 (MQ=255)
                                                                     tAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTc                                                                             >  1:1907517/1‑150 (MQ=255)
                                                                        cgggcggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCa                                                                          >  2:927456/1‑150 (MQ=255)
                                                                         gggcggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCaa                                                                         <  2:677787/150‑1 (MQ=255)
                                                                          ggcggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTg                                                                                                                                        <  2:28436/86‑1 (MQ=255)
                                                                          ggcggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTg                                                                                                                                        >  1:28436/1‑86 (MQ=255)
                                                                            cggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGc                                                                      <  1:927456/150‑1 (MQ=255)
                                                                             ggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTa                                                                                                              <  2:978169/109‑1 (MQ=255)
                                                                             ggTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTa                                                                                                              >  1:978169/1‑109 (MQ=255)
                                                                              gTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCAt                                                                    <  1:789241/150‑1 (MQ=255)
                                                                               ttGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCACCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGccgccg                                                                                                       <  1:1129663/114‑1 (MQ=255)
                                                                               ttGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCACCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGccgccg                                                                                                       >  2:1129663/1‑114 (MQ=255)
                                                                                    gAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTgg                                                              >  2:2382817/1‑150 (MQ=255)
                                                                                    gAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTgg                                                              <  1:2382817/150‑1 (MQ=255)
                                                                                     aGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGGGGTCCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTggg                                                             <  1:531865/150‑1 (MQ=255)
                                                                                      ccACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGc                                                            <  1:1193047/149‑1 (MQ=255)
                                                                                         cccAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGa                                                         >  2:1080203/1‑150 (MQ=255)
                                                                                          ccAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAGAGCATGTGTGGGCGGaa                                                        >  2:998469/1‑150 (MQ=255)
                                                                                          ccAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGaa                                                        >  1:1554555/1‑150 (MQ=255)
                                                                                          ccAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGaa                                                        >  1:2191118/1‑150 (MQ=255)
                                                                                          ccAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGGGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGaa                                                        >  1:1377277/1‑150 (MQ=255)
                                                                                               gTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTc                                                   >  2:169748/1‑150 (MQ=255)
                                                                                               gTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTc                                                   <  2:1258859/150‑1 (MQ=255)
                                                                                                ttGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTcg                                                  >  1:435987/1‑150 (MQ=255)
                                                                                                    gcCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGcc                                              >  1:877593/1‑150 (MQ=255)
                                                                                                     ccGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGAt                                                                                            >  1:73769/1‑103 (MQ=255)
                                                                                                     ccGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGAt                                                                                            <  2:73769/103‑1 (MQ=255)
                                                                                                           gTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGa                                       <  2:497787/150‑1 (MQ=255)
                                                                                                           gTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGa                                       <  2:1869605/150‑1 (MQ=255)
                                                                                                            ttCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAg                                      >  2:320398/1‑150 (MQ=255)
                                                                                                            ttCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAg                                      >  1:556061/1‑150 (MQ=255)
                                                                                                               gTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCg                                                           <  1:509823/126‑1 (MQ=255)
                                                                                                               gTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCg                                                           >  2:509823/1‑126 (MQ=255)
                                                                                                                ttGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACTCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTg                                  >  1:1388653/1‑150 (MQ=255)
                                                                                                                ttGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCgcgc                                                                                     >  1:741108/1‑99 (MQ=255)
                                                                                                                ttGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCgcgc                                                                                     <  2:741108/99‑1 (MQ=255)
                                                                                                                     cTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTgccgcc                             >  2:960807/1‑150 (MQ=255)
                                                                                                                      tGGTTGAGGATCAGCGCCGCCTGGGTGGGCGCACCCTgg                                                                                                                                           <  1:1874742/39‑1 (MQ=9)
                                                                                                                       ggTTGAGGATCAGCGCCGCCAGGTTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGcccc                           <  1:2173835/150‑1 (MQ=255)
                                                                                                                         ttGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAg                         <  1:998469/150‑1 (MQ=255)
                                                                                                                          tGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAgg                        <  1:320398/150‑1 (MQ=255)
                                                                                                                           gAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGt                       <  2:1864820/150‑1 (MQ=255)
                                                                                                                            aGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGCTCGGGAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTa                      >  1:635352/1‑150 (MQ=255)
                                                                                                                              gATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTc                                                                             >  1:2022875/1‑93 (MQ=255)
                                                                                                                              gATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTc                                                                             <  2:2022875/93‑1 (MQ=255)
                                                                                                                              gATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATg                    >  2:1867854/1‑150 (MQ=255)
                                                                                                                                tCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGtt                                   >  2:788013/1‑133 (MQ=255)
                                                                                                                                tCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGtt                                   <  1:788013/133‑1 (MQ=255)
                                                                                                                                  aGCGCCGCCAGGGGGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAACCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCTATGCCGAGGTTGCCGCCCGAGGTATGAGTg                >  1:928878/1‑150 (MQ=255)
                                                                                                                                   gcgcCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGc                                          <  2:1291478/123‑1 (MQ=255)
                                                                                                                                   gcgcCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGc                                          >  1:1291478/1‑123 (MQ=255)
                                                                                                                                   gcgcCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGa               >  2:1016261/1‑150 (MQ=255)
                                                                                                                                   gcgcCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGa               >  2:1671209/1‑150 (MQ=255)
                                                                                                                                    cgccgcCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGAt              <  1:1016261/150‑1 (MQ=255)
                                                                                                                                      caacctgggTGGGCGCACCCTGGTAGGTAGCCGCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGc            >  1:2232706/7‑150 (MQ=255)
                                                                                                                                      caacctgggTGGGCGCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCg                                        <  2:912773/116‑1 (MQ=255)
                                                                                                                                      caacctgggTGGGCGCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCg                                        >  1:912773/7‑122 (MQ=255)
                                                                                                                                      caacctgggTGGGCGCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGc            >  2:1953487/7‑150 (MQ=255)
                                                                                                                                       cgccAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCg           >  2:1552106/1‑150 (MQ=255)
                                                                                                                                       aacctgggTGGGCGCACCCTGGTTGGTAGCCCAGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCg           >  2:1296008/6‑150 (MQ=255)
                                                                                                                                       aacctgggTGGGCGCACCCTGGGTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCg           >  2:613134/6‑150 (MQ=255)
                                                                                                                                          cAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTAc        >  2:1464768/1‑150 (MQ=255)
                                                                                                                                          cAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTAc        >  1:1281226/1‑150 (MQ=255)
                                                                                                                                           aGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTAcc       <  1:1324329/150‑1 (MQ=255)
                                                                                                                                              gTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGt    <  1:2028261/150‑1 (MQ=255)
                                                                                                                                              gTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGt    >  1:386301/1‑150 (MQ=255)
                                                                                                                                              gTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCAATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCGGt    >  2:2209142/1‑150 (MQ=255)
                                                                                                                                               tGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGtt   <  1:1552106/150‑1 (MQ=255)
                                                                                                                                                gcgcCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGttt  <  1:748766/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GAAGCTCTGGATGATCGCGTTGTAGCCCAGCCAGCGGCAGGGCTCGGGAGGCCAGCTTGCCAGGCTCGATAGCGGGCGGTTGTCGAGCACCCAGGGTTGCGCCGTCAGTTCGTTGTGCTGGTTGAGGATCAGCGCCGCCAGGGTGCGCCCACCCAGGTTGGTAGCCCCGACGCCTTCGCCGCCGTAGCCGCCGGCAAGGGCGATGCCGCGCTGGCGGTCGCAAAGCATGTGTGGGCGGAAGCGTCGCGCCATGCCGAGGTTGCCGCCCCAGGTATGAGTGATGCGTACCTGTTT  >  NC_002947/1014912‑1015205

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: ‑