Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,128,188 C→G 5.7% A171A (GCC→GCG htrG → signal transduction protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,128,1880CG 5.7% 261.4 / 4.7 83A171A (GCC→GCGhtrGsignal transduction protein
Reads supporting (aligned to +/- strand):  ref base C (37/41);  new base G (3/2);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCG  >  NC_002947/4128040‑4128331
                                                                                                                                                    |                                                                                                                                               
ccATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGccc                                                                                                                                                <  2:1519168/150‑1 (MQ=255)
   tCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCTGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCGCGcc                                                                                                                                             >  1:2004724/1‑150 (MQ=255)
       cgacAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTggg                                                                                                                                         <  1:31881/150‑1 (MQ=255)
        gacAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACACAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTgggg                                                                                                                                        <  2:739847/150‑1 (MQ=255)
            gCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGAGGATg                                                                                                                                    <  2:276512/150‑1 (MQ=255)
             cTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                                                                   <  1:469779/150‑1 (MQ=255)
               ggAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATgaga                                                                                                                                 >  1:261696/1‑150 (MQ=255)
                gaagaaCCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCTCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGaa                                                                                                                                <  2:2300357/150‑1 (MQ=255)
                     aCCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                                                           <  2:261696/150‑1 (MQ=255)
                      ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGAt                                                                                                                                     <  2:901766/139‑1 (MQ=255)
                      ccGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGAt                                                                                                                                     >  1:901766/1‑139 (MQ=255)
                               gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                                                                   >  1:1531560/1‑132 (MQ=255)
                               gggCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATga                                                                                                                                   <  2:1531560/132‑1 (MQ=255)
                                   aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                                                             >  1:298550/1‑150 (MQ=255)
                                   aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                                                             >  1:365507/1‑150 (MQ=255)
                                   aTGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATgcgc                                                                                                             >  2:1000984/1‑150 (MQ=255)
                                    tGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                                                     >  2:360326/1‑141 (MQ=255)
                                    tGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                                                     <  1:360326/141‑1 (MQ=255)
                                               cTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTAt                                                                                                 >  1:1824863/1‑150 (MQ=255)
                                                   aTTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcg                                                                                             >  1:475923/1‑150 (MQ=255)
                                                      cgcgCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgca                                                                                          <  1:1632132/150‑1 (MQ=255)
                                                       gcgcAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCAATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgcag                                                                                         >  2:664050/1‑150 (MQ=255)
                                                        cgcAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGgcagc                                                                                        >  2:1995370/1‑150 (MQ=255)
                                                             gcagcTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                                                                   >  2:1946491/1‑150 (MQ=255)
                                                                 cTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGt                                                                               <  1:1816265/150‑1 (MQ=255)
                                                                 cTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGt                                                                               >  1:1384017/1‑150 (MQ=255)
                                                                   gattgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGc                                                                             >  2:591753/1‑150 (MQ=255)
                                                                    attgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCa                                                                            >  2:1582885/1‑150 (MQ=255)
                                                                     ttgatGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAg                                                                           <  2:1824863/150‑1 (MQ=255)
                                                                              tGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTgg                                                                  <  2:365507/150‑1 (MQ=255)
                                                                                     aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                                                         >  2:326447/1‑136 (MQ=255)
                                                                                     aCGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                                                         <  1:326447/136‑1 (MQ=255)
                                                                                       gCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGa                                                         >  2:820297/1‑150 (MQ=255)
                                                                                       gCAACAAGGCACTGAACGAGCAGCTCGATCAAAGCCAGTCTGATCTGCGCGATACCCAGGCCCGCCTGGGGGAGGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGTCGGTGCAGGCTTGCTGGCGGGGCTGa                                                         <  1:1102528/150‑1 (MQ=255)
                                                                                        caacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                                                     <  2:176288/89‑1 (MQ=255)
                                                                                        caacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGtga                                                                                                                     >  1:176288/1‑89 (MQ=255)
                                                                                         aacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                                                           >  1:1372668/1‑82 (MQ=255)
                                                                                         aacaaGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGc                                                                                                                           <  2:1372668/82‑1 (MQ=255)
                                                                                           caaGGCACTGAACGAGCAGCTCGATCCGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCt                                                     <  1:591753/150‑1 (MQ=255)
                                                                                             aGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGGATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGc                                                   >  1:1077685/1‑150 (MQ=255)
                                                                                                 aCTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtg                                                                                                    <  2:1671598/97‑1 (MQ=255)
                                                                                                 aCTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtg                                                                                                    >  1:1671598/1‑97 (MQ=255)
                                                                                                   tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcggc                                                                                           <  2:1752913/104‑1 (MQ=255)
                                                                                                   tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATggcggc                                                                                           >  1:1752913/1‑104 (MQ=255)
                                                                                                   tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCgg                                                <  1:308686/147‑1 (MQ=255)
                                                                                                   tGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCgg                                                >  2:308686/1‑147 (MQ=255)
                                                                                                    gatcGAGCTGCTCGTTCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCgg                                                               >  1:1289521/4‑131 (MQ=255)
                                                                                                    gatcGAGCTGCTCGTTCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCgg                                                               <  2:1289521/128‑1 (MQ=255)
                                                                                                    gAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACt                                            <  1:1582885/150‑1 (MQ=255)
                                                                                                     aaCGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCAGTTGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTg                                           >  1:1919978/1‑150 (MQ=255)
                                                                                                        gAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                                                                   >  2:1343083/1‑107 (MQ=255)
                                                                                                        gAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGc                                                                                   <  1:1343083/107‑1 (MQ=255)
                                                                                                          gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCg                                      >  1:2311501/1‑150 (MQ=255)
                                                                                                          gcagcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCg                                      <  1:1946491/150‑1 (MQ=255)
                                                                                                             gcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTgg                                   <  1:1285694/150‑1 (MQ=255)
                                                                                                             gcTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACCGACCCGTgg                                   <  2:1077685/150‑1 (MQ=255)
                                                                                                                 gATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGgcgc                               >  1:1649968/1‑150 (MQ=255)
                                                                                                                   tCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtgt                                                                                                   >  1:1733530/1‑80 (MQ=255)
                                                                                                                   tCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGgtgt                                                                                                   <  2:1733530/80‑1 (MQ=255)
                                                                                                                      gagCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCaagaa                          >  1:1383310/1‑150 (MQ=255)
                                                                                                                        gCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAagaaga                        >  2:164905/1‑150 (MQ=255)
                                                                                                                          cAGTCGGACCTGGGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCTCTACATCCTGTGTGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAAc                      <  1:1569871/150‑1 (MQ=255)
                                                                                                                          cAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAAc                      <  2:1383310/150‑1 (MQ=255)
                                                                                                                            gTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                                                         <  1:1974498/97‑1 (MQ=255)
                                                                                                                            gTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGc                                                                         >  2:1974498/1‑97 (MQ=255)
                                                                                                                                gACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGc                <  1:664050/150‑1 (MQ=255)
                                                                                                                                gACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGc                <  1:1995370/150‑1 (MQ=255)
                                                                                                                                  ccTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTg              <  2:1935014/150‑1 (MQ=255)
                                                                                                                                   cTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTgg             <  1:820297/150‑1 (MQ=255)
                                                                                                                                    tGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGgcg                                >  1:1412435/1‑130 (MQ=255)
                                                                                                                                    tGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGgcg                                <  2:1412435/130‑1 (MQ=255)
                                                                                                                                     gcgcGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGAt                                                        <  2:155183/105‑1 (MQ=255)
                                                                                                                                     gcgcGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGAt                                                        >  1:155183/1‑105 (MQ=255)
                                                                                                                                        cgATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTg        <  2:298550/150‑1 (MQ=255)
                                                                                                                                         gATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGa       <  2:2004724/150‑1 (MQ=255)
                                                                                                                                         gATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGa       <  2:1384017/150‑1 (MQ=255)
                                                                                                                                         atcccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCtt                                                                       >  2:1780434/5‑86 (MQ=255)
                                                                                                                                         atcccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCtt                                                                       <  1:1780434/82‑1 (MQ=255)
                                                                                                                                            ccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAAc                      >  2:1506507/2‑132 (MQ=255)
                                                                                                                                            ccccAGGCGGGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAAc                      <  1:1506507/131‑1 (MQ=255)
                                                                                                                                             cccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTc   >  1:1379490/1‑150 (MQ=255)
                                                                                                                                              ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGc                >  2:436102/1‑136 (MQ=255)
                                                                                                                                              ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGc                <  1:436102/136‑1 (MQ=255)
                                                                                                                                              ccAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCg  >  1:604903/1‑150 (MQ=255)
                                                                                                                                                    ccGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCt         <  1:2299556/137‑1 (MQ=255)
                                                                                                                                                    ccGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCt         >  2:2299556/1‑137 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCATCGACGACAGCTGGAAGAACCGCGTACAGGGCATGCAGGAAACGCTTGATTCGCGCAAGCAGCTGATTGATGAGTTGGAGGCACGCAACAAGGCACTGAACGAGCAGCTCGATCAGAGCCAGTCGGACCTGCGCGATACCCAGGCCCGCCTGGGGGATGAGAACAAGCAGGTGATGATGCGCTACATGGTGTATGGCGGCAGCATTGCCGGTGCAGGCTTGCTGGCGGGGCTGATCCTGCCGGCACTGACCCGTGGGCGCAAGAAGAACGATCGCTGGTTCTGAGCTCG  >  NC_002947/4128040‑4128331

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: