Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,261,125 A→C 7.2% L175R (CTG→CGG)  PP_4637 ← 5‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,261,1250AC 7.2% 311.1 / 12.9 110L175R (CTG→CGG) PP_46375‑methyltetrahydropteroyltriglutamate‑ homocysteine S‑methyltransferase family protein
Reads supporting (aligned to +/- strand):  ref base A (55/47);  new base C (3/5);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

Only 100 of 110 total aligned reads displayed.

CCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCAT  >  NC_002947/5260977‑5261272
                                                                                                                                                    |                                                                                                                                                     
ccccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                      <  1:2359620/150‑1 (MQ=255)
  ccAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGt                                                                                                                                                    >  2:1010040/1‑150 (MQ=255)
     gggTCTTGCCCAAAACACACGGCGAGGCGTTGCGGCCGTCGACAAGTCCCAGCGACAGTACCTTGTGGGCCGCCAGGCGATCGAGAATGGTCCGGTACTGCTCGGGCGCCTGAACCAGATCCACGTGCAGGCCATCAACCGGCAGGTTgg                                                                                                                                                 <  2:161375/150‑1 (MQ=255)
       gTCGTGCCCAAAACACAAGGCGAGGCGTTGAGACCGCCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                               <  2:501997/150‑4 (MQ=255)
         cTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatggc                                                                                                                                             <  2:2193052/150‑6 (MQ=255)
          ttGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcc                                                                                                                                            >  2:157525/1‑150 (MQ=255)
              ccAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGc                                                                                                                                        >  1:2100410/1‑150 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatgg                                                                                                                                              >  1:315412/1‑139 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCGCCAACCTGCCGGTTgatgg                                                                                                                                              <  2:315412/143‑5 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGcctgcac                                                                                                                                       <  1:453397/150‑6 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGcc                                                                                                                                       >  1:584751/1‑150 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGcc                                                                                                                                       >  1:578768/1‑150 (MQ=255)
               cAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGcc                                                                                                                                       >  2:1450468/1‑150 (MQ=255)
                aaaTCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCg                                                                                                                                      >  1:1388718/1‑150 (MQ=255)
                      caAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTg                                                                                                                                >  2:2194673/1‑150 (MQ=255)
                       aaGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGc                                                                                                                               <  1:1540897/150‑1 (MQ=255)
                          gcgaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatggc                                                                                                                                             >  1:2290360/1‑128 (MQ=255)
                            gaggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTcc                                                                                                                          <  1:2194673/150‑1 (MQ=255)
                              ggcgTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAg                                                                                                                        <  2:240179/150‑1 (MQ=255)
                                       aCCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGaa                                                                                                               <  1:1801103/150‑1 (MQ=255)
                                       aCCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCGTCCAGGCCACCGaa                                                                                                               <  2:740406/150‑1 (MQ=255)
                                         cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                             >  1:1708581/1‑150 (MQ=255)
                                         cGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGGTCGACGTGCAGGCAATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                             >  2:1303501/1‑150 (MQ=255)
                                          gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAgg                                                                                                                       >  2:1248934/1‑139 (MQ=255)
                                          gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAgg                                                                                                                       <  1:1248934/139‑1 (MQ=255)
                                          gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTa                                                                                                            >  2:182905/1‑150 (MQ=255)
                                          gTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTa                                                                                                            >  1:2241778/1‑150 (MQ=255)
                                             gACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGt                                                                                                         >  2:2054468/1‑150 (MQ=255)
                                               ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                      >  2:1856229/1‑103 (MQ=255)
                                               ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAg                                                                                                                                                      <  1:1856229/103‑1 (MQ=255)
                                               ccAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGCCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTgg                                                                                                       >  2:879202/1‑150 (MQ=255)
                                                 aGGCCCAGCGACAGTACCTTGTAGGACGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCg                                                                                                     >  1:2269970/1‑150 (MQ=255)
                                                     ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                                 >  2:1132923/1‑150 (MQ=255)
                                                     ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                                 >  1:932740/1‑150 (MQ=255)
                                                     ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                                 >  2:861380/1‑150 (MQ=255)
                                                     ccAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCa                                                                                                 >  1:163566/1‑150 (MQ=255)
                                                       aGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCAAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                               <  2:2301831/150‑1 (MQ=255)
                                                       aGCGACAGTACCTGGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGc                                                                                               <  2:730947/150‑1 (MQ=255)
                                                         cGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGAGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGctt                                                                                             <  1:25369/150‑1 (MQ=255)
                                                               tACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGc                                                                                                                      <  2:738081/119‑1 (MQ=255)
                                                               tACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGc                                                                                                                      >  1:738081/1‑119 (MQ=255)
                                                                  cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                               <  1:895200/91‑4 (MQ=255)
                                                                  cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCTGCCGGTTgatg                                                                                                                                               >  2:895200/1‑88 (MQ=255)
                                                                  cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGcttc                                                                                            <  1:999333/142‑1 (MQ=255)
                                                                  cTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACgg                                                                                  >  1:983301/1‑150 (MQ=255)
                                                                   ttGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGt                                                                                 <  2:942380/150‑1 (MQ=255)
                                                                      tAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTgcg                                                                              >  1:1109926/1‑150 (MQ=255)
                                                                       aGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTgcgc                                                                             <  2:2170066/150‑1 (MQ=255)
                                                                         gCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGt                                                                           >  2:243506/1‑150 (MQ=255)
                                                                             gCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTgta                                                                       <  1:1098935/150‑3 (MQ=255)
                                                                             gCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGaa                                                                       >  1:2098297/1‑150 (MQ=255)
                                                                               aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAgg                                                                     <  2:2098297/150‑1 (MQ=255)
                                                                               aGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAgg                                                                     >  1:2229403/1‑150 (MQ=255)
                                                                                 gCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGCGCGCGTTGAAGGAt                                                                   <  2:393841/150‑1 (MQ=255)
                                                                                   gATCGAGAATGGTCTGGTACTGCGCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGAGTCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAggttct                                                                 >  1:500505/1‑146 (MQ=255)
                                                                                    aTCGAGAGTGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGAtgtt                                                                <  1:1823300/150‑1 (MQ=255)
                                                                                     tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                                 <  2:1537368/86‑1 (MQ=255)
                                                                                     tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGAtgttg                                                               <  1:243506/150‑1 (MQ=255)
                                                                                     tCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGAtgttg                                                               <  2:445839/150‑1 (MQ=255)
                                                                                      cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGata                                                                                                                                   <  1:1913504/83‑3 (MQ=255)
                                                                                      cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                             <  2:710288/105‑1 (MQ=255)
                                                                                      cGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGt                                                                                                             >  1:710288/1‑105 (MQ=255)
                                                                                              gTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTAGCCGTCCAGGCCACCGAAGTCAGTGATGCTATTTCGGCGGTACTGTGCGCGTTGAAGGATGTTGTAGACGCGt                                                      >  1:2377163/1‑150 (MQ=255)
                                                                                               tCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGtt                                                     >  1:1775739/1‑150 (MQ=255)
                                                                                                 tGGTACTGCTCGGGCGCCTGAACCATATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCa                                                   <  2:1708581/150‑1 (MQ=255)
                                                                                                   gTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATa                                                 <  1:1450468/150‑1 (MQ=255)
                                                                                                    tACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAg                                                <  1:1010040/150‑1 (MQ=255)
                                                                                                     aCTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                         <  1:1174146/78‑1 (MQ=255)
                                                                                                     aCTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCa                                                                                                                         >  2:1174146/1‑78 (MQ=255)
                                                                                                      cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTggc                                                                                                                                                >  1:701706/1‑54 (MQ=255)
                                                                                                      cTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGCGGCTATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGc                                              <  2:2432681/150‑1 (MQ=255)
                                                                                                       tGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTGGTATCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCg                                             <  1:157525/150‑1 (MQ=255)
                                                                                                        gCTCCGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGt                                            <  2:949471/150‑1 (MQ=255)
                                                                                                             ggCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTc                                       <  2:2229403/150‑1 (MQ=255)
                                                                                                                  cTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTc                                  >  2:2382254/1‑150 (MQ=255)
                                                                                                                     aaCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGtt                                                                                                                                 >  2:332820/1‑54 (MQ=255)
                                                                                                                      aCCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGACGGATGTTGTAGACGCGTTCATATGCGTTTTTTCATTCCTgg                              <  2:983301/150‑1 (MQ=255)
                                                                                                                       ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTggg                             >  1:1030251/1‑150 (MQ=255)
                                                                                                                       ccAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTGTGTCCAGTCCTggg                             >  1:1696622/1‑150 (MQ=255)
                                                                                                                         aGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAg                                                                                                <  2:2059982/83‑1 (MQ=255)
                                                                                                                         aGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAg                                                                                                >  1:2059982/1‑83 (MQ=255)
                                                                                                                          gATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCa                          >  1:1582415/1‑150 (MQ=255)
                                                                                                                           aTCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAg                         >  2:1551503/1‑150 (MQ=255)
                                                                                                                                cGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCt                    <  1:1551503/150‑1 (MQ=255)
                                                                                                                                      ctccATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTAACGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGa              <  1:1453869/148‑1 (MQ=255)
                                                                                                                                       gCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAg             >  1:1816101/1‑150 (MQ=255)
                                                                                                                                         cATCAACCGGCAGGTTGTCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa           >  1:2136078/1‑150 (MQ=255)
                                                                                                                                         cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGCCGCGTTAATAGGGGATCTGCCTGCGCTGGGGCAGCGCCTGGGCGAgag           >  2:1453869/1‑148 (MQ=255)
                                                                                                                                         cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa           >  1:380078/1‑150 (MQ=255)
                                                                                                                                         cATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGa           >  1:2032709/1‑150 (MQ=255)
                                                                                                                                          aTCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGAt          >  1:904813/1‑150 (MQ=255)
                                                                                                                                           tCAACCGGCAGGTTGTCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATg         <  2:2136078/150‑1 (MQ=255)
                                                                                                                                            cAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATgg        <  1:792139/150‑1 (MQ=255)
                                                                                                                                            cAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATgg        <  2:2211461/150‑1 (MQ=255)
                                                                                                                                             aaCCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATggg       <  1:70503/150‑1 (MQ=255)
                                                                                                                                                cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGtt                                                                          >  1:1894664/1‑80 (MQ=255)
                                                                                                                                                cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGtt                                                                          <  2:1894664/80‑1 (MQ=255)
                                                                                                                                                cGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTc    >  1:1703984/1‑150 (MQ=255)
                                                                                                                                                 ggCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCa   >  1:1122922/1‑150 (MQ=255)
                                                                                                                                                  gCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCAt  >  2:2034114/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCCCAGGGTCTTGCCCAAATCACAAGGCGAGGCGTTGCGACCGTCGACCAGGCCCAGCGACAGTACCTTGTAGGCCGGCAGGCGATCGAGAATGGTCTGGTACTGCTCGGGCGCCTGAACCAGATCGACGTGCAGGCCATCAACCGGCAGGTTGGCGGCCAGGCCGAGGTTGCCTTCCAGGCCACCGAAGTACGTGGCGATCAGCTTCTTC‑‑AACGGTGCGCGTTGAAGGATGTTGTAGACGCGTTCATAGGCGTTTTTCCAGTCCTGGGGCAGCGCCTGGGCGAGGATGGGTTCAT  >  NC_002947/5260977‑5261272

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: