Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,378,190 G→C 6.1% E66Q (GAG→CAG)  PP_2978 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,378,1900GC 6.1% 311.1 / 8.1 98E66Q (GAG→CAG) PP_2978hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (47/45);  new base C (4/2);  total (51/47)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CAAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCACG  >  NC_002947/3378042‑3378339
                                                                                                                                                      |                                                                                                                                                     
cAAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGa                                                                                                                                                      >  2:1468771/1‑150 (MQ=255)
cAAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGa                                                                                                                                                      >  2:452399/1‑150 (MQ=255)
 aaaGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAg                                                                                                                                                     >  1:2314660/1‑150 (MQ=255)
  aaGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGc                                                                                                                                                    >  1:770970/1‑150 (MQ=255)
  aaGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGACAAGCACCCAACATCAGGTTGTTCCAGATTTCCCGCCTGAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGc                                                                                                                                                    <  1:1596712/150‑1 (MQ=255)
   aGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGAGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCa                                                                                                                                                   <  2:2913991/150‑1 (MQ=255)
          gAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGcccc                                                                                                                                            >  2:1578013/1‑150 (MQ=255)
          gAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGcccc                                                                                                                                            >  2:2453512/1‑150 (MQ=255)
          gAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGcccc                                                                                                                                            >  1:2194560/1‑150 (MQ=255)
            gCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGt                                                                                                                                          <  2:1859648/150‑1 (MQ=255)
             cGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                         <  2:392506/150‑1 (MQ=255)
             cGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                         <  2:1368061/150‑1 (MQ=255)
              gTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCGGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                        >  1:2263211/1‑150 (MQ=255)
              gTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                        <  2:1943054/150‑1 (MQ=255)
                   gAAAGGCATGATGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtg                                                                                                                                   <  1:1666768/150‑1 (MQ=255)
                   gAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCTGCGCTGAGCACCGCCCCGTGGTAGtg                                                                                                                                   >  2:2598574/1‑150 (MQ=255)
                     aaGGAATGGCTAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtgat                                                                                                                                 <  1:519510/150‑1 (MQ=255)
                     aaGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGGAAGATAGTGCCACAGGCAGCGCAGTCAGCATCCCTTTCGAGCGAATTGTCGAGCGTAACGTCAAAGTCCGTTGTGCTGGCGGCGCTGGGCACCGCCCCGTGGTAGtgat                                                                                                                                 <  2:498012/150‑1 (MQ=255)
                        gAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGaa                                                                                                                              <  2:770970/150‑1 (MQ=255)
                           tGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTCTGAACTGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAt                                                                                                                           >  1:2464838/1‑150 (MQ=255)
                            ggCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtg                                                                                                                          <  1:298370/150‑1 (MQ=255)
                                      gggTGCGTTTCAAAGTTGAAGTGAAGGACAGTGACACGGGCAGAGCAGCCAGCATCCCTTTTGCGCGAAGGGTCGCGCGTCACGTCAAAGTACGTGCAGCTGTCGGCGCTGAGCACCGCCCCGCGGTAGTGATGAAGATGTGCATTGGCa                                                                                                                >  2:1105326/1‑150 (MQ=255)
                                       ggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                         >  2:2482267/1‑124 (MQ=255)
                                       ggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTg                                                                                                                                         <  1:2482267/124‑1 (MQ=255)
                                       ggTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCaa                                                                                                               <  2:2263211/150‑1 (MQ=255)
                                       ggTGCGTTTCAAAGCTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGTATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCaa                                                                                                               <  1:2606780/150‑1 (MQ=255)
                                        gTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAAc                                                                                                              >  1:2040250/1‑150 (MQ=255)
                                        gTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAAc                                                                                                              <  2:2131694/150‑1 (MQ=255)
                                         tGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgt                                                                                                                         >  1:748473/1‑138 (MQ=255)
                                         tGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGAtgt                                                                                                                         <  2:748473/138‑1 (MQ=255)
                                          gCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                                                                                                                                  <  1:2121180/128‑1 (MQ=255)
                                          gCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtga                                                                                                                                  >  2:2121180/1‑128 (MQ=255)
                                             tttCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCt                                                                                                         <  2:2314660/150‑1 (MQ=255)
                                                 aaaGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTa                                                                                                     <  1:180391/150‑1 (MQ=255)
                                                     ttGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAAc                                                                                                 <  1:911700/150‑1 (MQ=255)
                                                        aaGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGaa                                                                                              <  2:1639858/150‑1 (MQ=255)
                                                           tGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACaa                                                                                           >  1:2870013/1‑150 (MQ=255)
                                                               aGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATtct                                                                                       >  1:164195/1‑150 (MQ=255)
                                                                 aCAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     <  1:1468771/150‑1 (MQ=255)
                                                                 aCAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     >  2:1908246/1‑150 (MQ=255)
                                                                  cAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTccc                                                                                    <  1:656624/150‑1 (MQ=255)
                                                                   aGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCt                                                                                   >  2:2750904/1‑150 (MQ=255)
                                                                    gTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGTGCTCAGCgccgcc                                                                                                                                              >  2:1146890/1‑90 (MQ=255)
                                                                    gTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGTGCTCAGCgccgcc                                                                                                                                              <  1:1146890/90‑1 (MQ=255)
                                                                      gACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCTGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTgaa                                                                                <  1:2598574/150‑1 (MQ=255)
                                                                         acaGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCTGCGCTGAGCACCGACCCGCGGTAGTGCAGAAGATGTGCATTGGGAACCAGCTTTTACAACGAAAAATTTTCCCTgcggca                                                                             >  2:492308/1‑145 (MQ=255)
                                                                           aGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                        <  1:765835/89‑1 (MQ=255)
                                                                           aGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTgg                                                                                                                                        >  2:765835/1‑89 (MQ=255)
                                                                            ggCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCg                                                                          <  1:2453512/150‑1 (MQ=255)
                                                                            ggCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCg                                                                          >  2:2140906/1‑150 (MQ=255)
                                                                             gCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCgg                                                                         <  1:1578013/150‑1 (MQ=255)
                                                                               agagCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTCCTGGCGGCGCTGAGCACCGCCCCGTGGGAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGTGa                                                                       >  2:2779240/1‑150 (MQ=255)
                                                                                       cagcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTTAATAACCGGGAAAAAATGc                                                               <  2:2040250/150‑1 (MQ=255)
                                                                                        agcATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCt                                                              >  1:268028/1‑150 (MQ=255)
                                                                                             cccTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     <  1:2552407/122‑1 (MQ=255)
                                                                                             cccTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     >  2:2552407/1‑122 (MQ=255)
                                                                                             cccTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATc                                                         >  2:1600777/1‑150 (MQ=255)
                                                                                               cTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCa                                                                                                            >  1:3042569/1‑97 (MQ=255)
                                                                                               cTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCa                                                                                                            <  2:3042569/97‑1 (MQ=255)
                                                                                               cTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCg                                                       >  1:903865/1‑150 (MQ=255)
                                                                                               cTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCg                                                       >  2:316294/1‑150 (MQ=255)
                                                                                                tttCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtg                                                                                                                                   >  2:2242797/1‑73 (MQ=255)
                                                                                                tttCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGtg                                                                                                                                   <  1:2242797/73‑1 (MQ=255)
                                                                                                tttCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCgg                                                      >  1:2627932/1‑150 (MQ=255)
                                                                                                tttCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCgg                                                      >  2:1949518/1‑150 (MQ=255)
                                                                                                  tCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGtt                                                    <  1:1908246/150‑1 (MQ=255)
                                                                                                           gggTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCgg                                                      >  1:3032342/1‑139 (MQ=255)
                                                                                                           gggTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCgg                                                      <  2:3032342/139‑1 (MQ=255)
                                                                                                              tCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGt                                        >  2:2922680/1‑150 (MQ=255)
                                                                                                              tCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGt                                        >  1:1877112/1‑150 (MQ=255)
                                                                                                                          gTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCGCTGGa                             >  1:642973/1‑150 (MQ=255)
                                                                                                                          gTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGa                             >  2:1590762/1‑150 (MQ=255)
                                                                                                                          gTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCAGGGTTTGATCGGCCGTCGAGCACTGGa                             >  2:93049/1‑150 (MQ=255)
                                                                                                                           tCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCACTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAg                            >  2:555210/1‑150 (MQ=255)
                                                                                                                                gTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCt                       >  2:1804880/1‑150 (MQ=255)
                                                                                                                                gTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCt                       >  2:914274/1‑150 (MQ=255)
                                                                                                                                 tCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGa                             <  1:876550/143‑1 (MQ=255)
                                                                                                                                 tCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGa                             >  2:876550/1‑143 (MQ=255)
                                                                                                                                  ccGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAg                     <  2:164195/150‑1 (MQ=255)
                                                                                                                                  ccGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAg                     <  2:2627932/150‑1 (MQ=255)
                                                                                                                                    gtggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     <  2:1913454/84‑1 (MQ=255)
                                                                                                                                    gtggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTcc                                                                                     >  1:1913454/1‑84 (MQ=255)
                                                                                                                                     tggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTa                      <  2:945995/146‑1 (MQ=255)
                                                                                                                                     tggtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTa                      >  1:945995/1‑146 (MQ=255)
                                                                                                                                     accacggggcggtGCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAg                  >  1:274089/8‑150 (MQ=255)
                                                                                                                                      ccacggggcggtGCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGc                 >  2:1542690/7‑150 (MQ=255)
                                                                                                                                       gtgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCg                <  1:1600777/150‑1 (MQ=255)
                                                                                                                                        tgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGa                             <  2:1118119/136‑1 (MQ=255)
                                                                                                                                        tgCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGa                             >  1:1118119/1‑136 (MQ=255)
                                                                                                                                           tGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTgg            <  2:903865/150‑1 (MQ=255)
                                                                                                                                            ggcggtGCTCAGCGCCGCCCCGTGGTAGTGATTAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa           <  1:1542690/150‑1 (MQ=255)
                                                                                                                                            ggcggcGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGa           <  2:2870013/150‑1 (MQ=255)
                                                                                                                                              gcggtgCTCAGCGCCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGAt          >  1:2753802/1‑150 (MQ=255)
                                                                                                                                                ggcgCTGAGCACCGCCCCGTGGTAGGGATGCAGCTGGGCACTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGt        <  1:1105326/150‑1 (MQ=255)
                                                                                                                                                   gcTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATc     >  1:1158534/1‑150 (MQ=255)
                                                                                                                                                    cTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCa    <  2:2194560/150‑1 (MQ=255)
                                                                                                                                                     tGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCAc   >  2:2896031/1‑150 (MQ=255)
                                                                                                                                                     tGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCAc   <  1:452399/150‑1 (MQ=255)
                                                                                                                                                      gAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCACg  >  2:1546398/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CAAAGATTTGGAGCGTTTCGAAAGGAATGGCGAGAAATGGGTGCGTTTCAAAGTTGAAGTGAAAGACAGTGACACAGGCAGAGCAGCCAGCATCCCTTTCGAGCGAAGGGTCGAGCGTAACGTCAAAGTCCGTGGTGCTGGCGGCGCTGAGCACCGCCCCGTGGTAGTGATGAAGATGTGCATTGGCAACCAGCTTTTACAACGAACAATTCTCCCTGAAGAACCGGGAAAAAATGCTCCATCCGGTTTTGATCGGCCGTCGAGCACTGGAGCACCTAGGAGCGGTGGATGTATCACG  >  NC_002947/3378042‑3378339

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