Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,708,217 G→C 6.7% S154R (AGC→AGG)  PP_1503 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,708,2170GC 6.7% 235.5 / 6.3 74S154R (AGC→AGG) PP_1503hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (35/34);  new base C (2/3);  total (37/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

GATTCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTCGTCGTACTTCTT  >  NC_002947/1708069‑1708366
                                                                                                                                                    |                                                                                                                                                     
gATTCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAgcgc                                                                                                                                                      <  2:812500/150‑1 (MQ=255)
 aTTCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCt                                                                                                                                                     <  1:493741/150‑1 (MQ=255)
   tCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTgc                                                                                                                                                   <  2:154386/150‑1 (MQ=255)
   tCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGCGCTCAGTGAGCCCTgc                                                                                                                                                   <  2:1538022/150‑1 (MQ=255)
      aGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGTCGCACGGCTCACTAAGCGCTGCGCa                                                                                                                                                >  2:665183/1‑150 (MQ=255)
      aGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCa                                                                                                                                                <  1:1692055/150‑1 (MQ=255)
       ggTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAg                                                                                                                                               <  1:1175190/150‑1 (MQ=255)
               ttCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCt                                                                                                                                       >  1:347074/1‑150 (MQ=255)
               ttCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCt                                                                                                                                       >  2:1070775/1‑150 (MQ=255)
               ttCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCt                                                                                                                                       >  1:184179/1‑150 (MQ=255)
                      cATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCacgtac                                                                                                                                <  2:1674127/150‑1 (MQ=255)
                       aTGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAgg                                                                                                                                              <  1:1401053/135‑1 (MQ=255)
                       aTGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAgg                                                                                                                                              >  2:1401053/1‑135 (MQ=255)
                           gAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCt                                                                                                                           >  1:287415/1‑150 (MQ=255)
                               aggcagCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAAc                                                                                                                       <  1:1070775/150‑1 (MQ=255)
                                ggcagCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACg                                                                                                                      >  2:1690610/1‑150 (MQ=255)
                                 gcagcaAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGc                                                                                                                     >  1:2276014/1‑150 (MQ=255)
                                         aCCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCgg                                                                                                             <  2:14597/150‑1 (MQ=255)
                                               gCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGTGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATcc                                                                                                       <  2:549087/150‑1 (MQ=255)
                                                 tCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTg                                                                                                     >  1:48899/1‑150 (MQ=255)
                                                    gggCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCAc                                                                                                  >  2:2048185/1‑150 (MQ=255)
                                                      gCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACgg                                                                                                <  2:2276014/150‑1 (MQ=255)
                                                       ccAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGc                                                                                               <  1:733159/150‑1 (MQ=255)
                                                        cagcagGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGCGCTCAGTGAGCCCTgc                                                                                                                                                   >  2:74958/1‑97 (MQ=255)
                                                        cagcagGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGCGCTCAGTGAGCCCTgc                                                                                                                                                   <  1:74958/97‑1 (MQ=255)
                                                        cagcagGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTATGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAAAGACCTGGCGAGCATCCTGCACGTcc                                                                                              >  2:1123254/1‑150 (MQ=255)
                                                          gcagGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAgg                                                                                                                                              <  1:2409177/100‑1 (MQ=255)
                                                          gcagGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAgg                                                                                                                                              >  2:2409177/1‑100 (MQ=255)
                                                              ggCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCgg                                                                                                             <  2:715280/129‑1 (MQ=255)
                                                              ggCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCgg                                                                                                             >  1:715280/1‑129 (MQ=255)
                                                                 ggggATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCt                                                                                     >  2:406750/1‑150 (MQ=255)
                                                                  ggAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCtt                                                                                    >  2:1393962/1‑150 (MQ=255)
                                                                     aTCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTGCGGGTCTTTgg                                                                                 >  2:1341537/1‑150 (MQ=255)
                                                                         ttACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGc                                                                             >  2:833892/1‑150 (MQ=255)
                                                                         ttACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGc                                                                             <  1:406750/150‑1 (MQ=255)
                                                                         ttACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGc                                                                             >  1:2162823/1‑150 (MQ=255)
                                                                           aCGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCt                                                                           >  2:2336406/1‑150 (MQ=255)
                                                                                 ttGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCt                                                                     <  1:2336406/150‑1 (MQ=255)
                                                                                  tGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCtt                                                                    <  1:2006548/150‑1 (MQ=255)
                                                                                   gcgcATAGTGCTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTgg                                                                                        >  2:464165/1‑129 (MQ=255)
                                                                                   gcgcATAGTGCTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTgg                                                                                        <  1:464165/129‑1 (MQ=255)
                                                                                           tgtTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGc                                                                       >  1:318376/1‑138 (MQ=255)
                                                                                           tgtTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGc                                                                       <  2:318376/138‑1 (MQ=255)
                                                                                             ttCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGa                                                         >  1:2210491/1‑150 (MQ=255)
                                                                                                 tCCTTTGATCGTGTGCGAAGGCTAAGGTTGGAGCGCGAGGGGCCTTTGAGCGATACGCAGGTCTTCCTCACGGACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCg                                                     >  1:658749/1‑150 (MQ=255)
                                                                                                     ttGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCa                                                 <  2:184179/150‑1 (MQ=255)
                                                                                                              ggCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCg                                        <  2:985277/149‑1 (MQ=255)
                                                                                                                cGAAGGCTAAGGTTGGCAAGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGc                                      >  2:1234508/1‑150 (MQ=255)
                                                                                                                 gAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGcc                                     >  1:1576968/1‑150 (MQ=255)
                                                                                                                  aaGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCgg                                                                                                             <  2:1568047/77‑1 (MQ=255)
                                                                                                                  aaGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCgg                                                                                                             >  1:1568047/1‑77 (MQ=255)
                                                                                                                  aaGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGc                                      <  1:2145410/148‑1 (MQ=255)
                                                                                                                  aaGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGc                                      >  2:2145410/1‑148 (MQ=255)
                                                                                                                     gCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGtt                                 <  2:1576968/150‑1 (MQ=255)
                                                                                                                     gCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGtt                                 <  2:48899/150‑1 (MQ=255)
                                                                                                                     gCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGtt                                 >  1:817560/1‑150 (MQ=255)
                                                                                                                     gCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGtt                                 <  1:665183/150‑1 (MQ=255)
                                                                                                                       tAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGc                                                      >  2:1095500/1‑127 (MQ=255)
                                                                                                                       tAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGc                                                      <  1:1095500/127‑1 (MQ=255)
                                                                                                                        aaGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTggg                              <  1:1393962/150‑1 (MQ=255)
                                                                                                                               gCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTg                                >  2:204593/1‑141 (MQ=255)
                                                                                                                               gCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTg                                <  1:204593/141‑1 (MQ=255)
                                                                                                                                ctgcgCAGCGCTCAGTGAGCCCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGc                      >  2:1933465/5‑150 (MQ=255)
                                                                                                                                cGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGc                      <  2:2162823/150‑1 (MQ=255)
                                                                                                                                   ggCAGGGCTCACTGAGCGCTGCGCAGGCNNTCCTCACGTACGTCNTNANCGCCCTGGCGGGCATACTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGCGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTg                   >  1:1446412/1‑150 (MQ=255)
                                                                                                                                   ggCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTg                   <  1:2048185/150‑1 (MQ=255)
                                                                                                                                   ggCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCGGACCGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTg                   <  1:1341537/150‑1 (MQ=255)
                                                                                                                                      aGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGcc                >  2:1551277/1‑150 (MQ=255)
                                                                                                                                        ggCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCtc              <  2:658749/150‑1 (MQ=255)
                                                                                                                                          cTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCtcgt            >  2:751396/1‑150 (MQ=255)
                                                                                                                                           tCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCtcgtc           <  1:1690610/150‑1 (MQ=255)
                                                                                                                                             aCTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTg                                >  1:604595/1‑127 (MQ=255)
                                                                                                                                             aCTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTg                                <  2:604595/127‑1 (MQ=255)
                                                                                                                                             aCTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTcgtcgt         >  2:942545/1‑150 (MQ=255)
                                                                                                                                               tGAGCCCTGCGCAGGCCTTCCTCACGTACGTCCTGAACTTCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTCGTCGTAc       <  1:1933465/150‑1 (MQ=255)
                                                                                                                                                  gcgcTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTCGTCGTActtc    <  1:1123254/150‑1 (MQ=255)
                                                                                                                                                    gcTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTCGTCGTActtctt  >  2:367749/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GATTCCAGGTTGGCATTCGGATCATGCGAGCAGGCAGCAAGACCGACGCTCAGGGCCAGCAGGGCGGGAATCATTACGTGTTTGCGCATAGTGTTCATCCTTTGATCGTGTGCGAAGGCTAAGGTTGGCGGGGCAGGGCTCACTGAGCGCTGCGCAGGCCTTCCTCACGTACGTCCTGAACGCCCTGGCGGGCATCCTGCACGGCCTTCTGGGCTTTGGCTGCCTGGGCCTTGCGTTCGGCGAGGCGTGCATCCCACTCGGCCTGTTGGGCCAGTCGCTTGGCCTCGTCGTACTTCTT  >  NC_002947/1708069‑1708366

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