Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,375,488 A→T 5.4% V31D (GTC→GAC)  PP_5637 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,375,4880AT 5.4% 409.4 / 9.9 142V31D (GTC→GAC) PP_5637hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (70/64);  new base T (3/5);  total (73/69)
Fisher's exact test for biased strand distribution p-value = 4.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 147 total aligned reads displayed.

CACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATAAAG  >  NC_002947/4375341‑4375637
                                                                                                                                                   |                                                                                                                                                        
cacGACTGCGCATATAGCCCGCCTCCGACCTCTTCGACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACt                                                                                                                                                        <  2:501985/150‑1 (MQ=255)
   gACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGa                                                                                                                                                     <  2:2497898/150‑1 (MQ=255)
   gACTGCGCATATAGCACGCCTTCGACCTCCTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGa                                                                                                                                                     >  1:1118053/1‑150 (MQ=255)
      tGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGAccc                                                                                                                                                  >  1:772568/1‑150 (MQ=255)
        cgcATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCag                                                                                                                                                <  2:2535200/150‑1 (MQ=255)
          cATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCagag                                                                                                                                              <  2:1118053/150‑1 (MQ=255)
           atatAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCaggcc                                                                                                                                             <  2:772568/150‑4 (MQ=255)
               agccacccTTCGACCTCTTCTACACCGATGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTAGCATCGTGAAGCCGGGGGATAATGTCGGGCCCAGGAGTCGGGATAAACACCCAAGCCTGGGTCGCGACTCGACCCAGAGGCTCg                                                                                                                                         >  2:2931634/7‑150 (MQ=255)
                gCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGc                                                                                                                                        >  2:2603844/1‑150 (MQ=255)
                   cGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCtgt                                                                                                                                     <  2:574913/150‑1 (MQ=255)
                      cTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCGCGCTGTGAt                                                                                                                                  >  2:1467460/1‑150 (MQ=255)
                       ttCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCAACGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATg                                                                                                                                 >  1:1995231/1‑150 (MQ=255)
                              tcttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCa                                                                                                                          >  2:1171844/1‑150 (MQ=255)
                               cttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCaa                                                                                                                         >  2:422003/1‑150 (MQ=255)
                               cttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCaa                                                                                                                         >  2:571582/1‑150 (MQ=255)
                                ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCaggcct                                                                                                                                            >  2:1661560/1‑126 (MQ=255)
                                ttctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGAGTCGCGACCCaggcct                                                                                                                                            <  1:1661560/130‑5 (MQ=255)
                                 tctACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACACATAGGCTCGCTGTGATGGTCGCCAAta                                                                                                                       <  1:1950874/150‑2 (MQ=255)
                                   tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGt                                                                                                                               >  2:1560763/1‑140 (MQ=255)
                                   tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGt                                                                                                                               <  1:1560763/140‑1 (MQ=255)
                                   tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                                     <  2:1216942/150‑1 (MQ=255)
                                   tACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                                     <  2:93260/150‑1 (MQ=255)
                                      cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAgg                                                                                                                                             >  1:713648/1‑123 (MQ=255)
                                      cccGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAgg                                                                                                                                             <  2:713648/123‑1 (MQ=255)
                                         gttgatgttgatGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCtgtg                                                                                                                                    <  1:749016/129‑1 (MQ=255)
                                         gttgatgttgatGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCtgtg                                                                                                                                    >  2:749016/1‑129 (MQ=255)
                                         gttgatgttgatGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTc                                                                                                               >  2:1245952/1‑150 (MQ=255)
                                              tgttgatgCCGAGCTGGGCCATGCCGTCGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                          >  2:2198730/1‑150 (MQ=255)
                                                ttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCa                                                                                                                                                 <  1:238129/109‑1 (MQ=255)
                                                ttgatgCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCa                                                                                                                                                 >  2:238129/1‑109 (MQ=255)
                                                     gCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTAGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAg‑‑ggg                                                                                                 <  1:1509602/150‑1 (MQ=255)
                                                      ccGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                        <  1:3040809/128‑1 (MQ=255)
                                                      ccGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAt                                                                                                                        >  2:3040809/1‑128 (MQ=255)
                                                       cGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCt                                                                                               <  1:2749790/150‑1 (MQ=255)
                                                       cGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCt                                                                                               >  2:1657151/1‑150 (MQ=255)
                                                       cGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCt                                                                                               >  1:179579/1‑150 (MQ=255)
                                                          gCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGa                                                                                            <  2:785127/150‑1 (MQ=255)
                                                              ggCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACttt                                                                                       >  2:1216363/1‑150 (MQ=255)
                                                                ccATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTcc                                                                                     >  1:2167119/1‑150 (MQ=255)
                                                                 cATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                                     <  1:1409786/120‑1 (MQ=255)
                                                                 cATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAAtgtg                                                                                                                     >  2:1409786/1‑120 (MQ=255)
                                                                 cATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCg                                                                                    >  1:2528429/1‑150 (MQ=255)
                                                                    gCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAg                                                                                 <  2:2057112/150‑1 (MQ=255)
                                                                    gCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAg                                                                                 <  2:1746993/150‑1 (MQ=255)
                                                                     ccGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCATAGGTAGCGACAAGACCAAGAGGAAAGCGGTGATGGTAGCGAATGTGGACGTCCTCGCGGAG‑‑GGGCGGGACTTTCCGTAgg                                                                                >  2:672422/1‑150 (MQ=255)
                                                                            cgttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTg                                                                         >  2:1287729/3‑150 (MQ=255)
                                                                             tttCGCATCGTGCAGAGGGGGGAACATGTCTTTCCCATTCCGAGAGTCAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACATCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGa                                                                        <  2:1870639/149‑1 (MQ=255)
                                                                             cTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGa                                                                        >  2:3014460/1‑150 (MQ=255)
                                                                              ttCGCATCGTGCAGCCTGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                       <  2:1257596/150‑1 (MQ=255)
                                                                              ttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                          >  1:2162808/1‑118 (MQ=255)
                                                                              ttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCcgcgc                                                                                                          <  2:2162808/118‑1 (MQ=255)
                                                                              ttCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGGTCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                       <  1:424858/150‑1 (MQ=255)
                                                                                    tCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                                 <  1:571582/150‑1 (MQ=255)
                                                                                              gggggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCAt                                                       >  1:850812/1‑150 (MQ=255)
                                                                                              gggggATAATGTCGGGCCCAGGCGTAGGCTTAAACCCCCCGGCCGGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCAt                                                       <  1:1071151/150‑1 (MQ=255)
                                                                                                gggATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGc                                                     >  2:2867559/1‑150 (MQ=255)
                                                                                                    tAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgc                                                 >  1:1853943/1‑150 (MQ=255)
                                                                                                     aaTGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCGGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgcg                                                >  2:2279160/1‑150 (MQ=255)
                                                                                                       tGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcggc                                              <  1:2603311/150‑1 (MQ=255)
                                                                                                       tGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCgcggc                                              <  1:2603844/150‑1 (MQ=255)
                                                                                                         tCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCGGTGATGGTCGCGAATATGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGt                                            >  1:1905889/1‑150 (MQ=255)
                                                                                                          cGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTc                                           <  1:1627946/150‑1 (MQ=255)
                                                                                                             gCCCAGGCGTACGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                        <  2:1474724/150‑1 (MQ=255)
                                                                                                              cccAGGCGTAGGGTTACACACCCAGGCCTGGGTCGCGACTCGACCCAGACGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GTGCTGGACTTTCCGTAGCCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAAc                                       <  2:1096968/150‑1 (MQ=255)
                                                                                                              cccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAAc                                       <  1:1287729/150‑1 (MQ=255)
                                                                                                               ccAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCGCGCAATTCTCGGGGAGTCGAACg                                      >  2:1301417/1‑150 (MQ=255)
                                                                                                                 aGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACgag                                    <  1:522566/150‑1 (MQ=255)
                                                                                                                     gTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaagg                                >  1:1396711/1‑150 (MQ=255)
                                                                                                                      tAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaagga                               >  1:3000610/1‑150 (MQ=255)
                                                                                                                       aGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                        <  1:926542/140‑1 (MQ=255)
                                                                                                                       aGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGaa                                        >  2:926542/1‑140 (MQ=255)
                                                                                                                       aGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaaggaa                              <  2:1853943/150‑1 (MQ=255)
                                                                                                                         ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTg                                                                                              <  1:696633/85‑1 (MQ=255)
                                                                                                                         ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTg                                                                                              >  2:696633/1‑85 (MQ=255)
                                                                                                                         ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                       >  1:2644172/1‑107 (MQ=255)
                                                                                                                         ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAg                                                                       <  2:2644172/107‑1 (MQ=255)
                                                                                                                         ggTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGt                            >  1:1444057/1‑150 (MQ=255)
                                                                                                                             aaaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGAc                                                                                          <  1:347763/84‑1 (MQ=255)
                                                                                                                             aaaCACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGAc                                                                                          >  2:347763/1‑84 (MQ=255)
                                                                                                                               acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                      <  1:3014460/150‑1 (MQ=255)
                                                                                                                               acacCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGc                      <  2:2733992/150‑1 (MQ=255)
                                                                                                                                 acCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCgg                    >  2:1383553/1‑150 (MQ=255)
                                                                                                                                  cccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGc                                                                               <  1:1732790/90‑1 (MQ=255)
                                                                                                                                  cccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGc                                                                               >  2:1732790/1‑90 (MQ=255)
                                                                                                                                   ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCg                                                                              >  1:1876124/1‑90 (MQ=255)
                                                                                                                                   ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCg                                                                              <  2:1876124/90‑1 (MQ=255)
                                                                                                                                   ccAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgc                                                 >  2:485427/1‑119 (MQ=255)
                                                                                                                                       gCCTGGGTCGCGACTCGACCCAGAGGCTCTATGTGATTTTTTCCAATGTGGCCCTCCGCGCGGAGCCGAGCTGGA‑‑GTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa              <  1:1301417/150‑1 (MQ=255)
                                                                                                                                       gCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACCTCATCTGCGGAGGTGTAGGGCAAGGCCACGCCCTCACCGCGGCGCGCAACGAGAAGGAAGTACATGCGGCACCaa              >  1:2065472/1‑150 (MQ=255)
                                                                                                                                       gCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCaa              >  1:2507804/1‑150 (MQ=255)
                                                                                                                                        tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgc                                                 >  2:2687412/2‑114 (MQ=255)
                                                                                                                                        tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTcgc                                                 <  1:2687412/113‑1 (MQ=255)
                                                                                                                                        tcTGGGTCGAGTCGCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAgaa                                  <  1:1129709/128‑1 (MQ=255)
                                                                                                                                         cTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATa            >  1:2292085/1‑150 (MQ=255)
                                                                                                                                          tGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATAc           >  1:338034/1‑150 (MQ=255)
                                                                                                                                          tGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATAc           >  2:2864537/1‑150 (MQ=255)
                                                                                                                                           gggTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACt          <  1:422003/150‑1 (MQ=255)
                                                                                                                                            ggTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTg         <  2:2528429/150‑1 (MQ=255)
                                                                                                                                              tCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAg                                                                 <  1:2918413/92‑1 (MQ=255)
                                                                                                                                                gcgACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGtata     >  2:1627586/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                        
CACGACTGCGCATATAGCACGCCTTCGACCTCTTCTACCCCGTTGATGTTGATGCCGAGCTGGGCCATGCCGTTGACCTTCGCATCGTGCAGCCGGGGGATAATGTCGGGCCCAGGCGTAGGGTTAAACACCCAGGCCTGGGTCGCGACTCGACCCAGAGGCTCGCTGTGATGGTCGCCAATGTGGACGTCCGCGCGGAG‑‑GGGCTGGACTTTCCGTAGGCGGTCTGAGGGCAAGGCCACGCCATGCTCGCGGCGTCGAACGAGAAGGAAGTACATGCGGCACCAATACTGTATAAAG  >  NC_002947/4375341‑4375637

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