Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,424,904 A→G 5.2% *409W (TGA→TGG
S413P (TCA→CCA) 
moeA →
PP_2124 ←
molybdenum::molybdopterin ligase
glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,424,9040AG 5.2% 552.7 / 11.6 173*409W (TGA→TGG
S413P (TCA→CCA) 
moeA
PP_2124
molybdenum::molybdopterin ligase
glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (89/74);  new base G (5/4);  total (94/79)
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.97e-01

Only 100 of 176 total aligned reads displayed.

AGGGCCGCGTGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAACTGGCTCA  >  NC_002947/2424756‑2425053
                                                                                                                                                    |                                                                                                                                                     
aGGGCCGCGTGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAg                                                                                                                                                      >  2:1835891/1‑150 (MQ=255)
   gCCGCGTGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCgg                                                                                                                                                   <  1:2315777/150‑1 (MQ=255)
    ccGCGTGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGc                                                                                                                                                  >  1:4888004/1‑150 (MQ=255)
         tGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCccagc                                                                                                                                             <  1:2294602/150‑1 (MQ=255)
         tGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCccagc                                                                                                                                             >  1:714335/1‑150 (MQ=255)
            gcATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCccagccag                                                                                                                                          >  1:1251588/1‑150 (MQ=255)
             cATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAgg                                                                                                                                         <  2:4499261/150‑1 (MQ=255)
              aTCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGc                                                                                                                                        <  2:1153444/150‑1 (MQ=255)
               tCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACCTGGCTGGGCCGCtcagccaagta                                                                                                                                       <  2:708068/150‑5 (MQ=255)
               cccacacGAACCAGAGCTCGGGGGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCGGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACCTGGCTGGGCCGCtcagccaagta                                                                                                                                       <  2:4194787/147‑5 (MQ=255)
                cTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCAc                                                                                                                                      >  1:265264/1‑150 (MQ=255)
                cTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCAc                                                                                                                                      >  1:307843/1‑150 (MQ=255)
                 tACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACt                                                                                                                                     >  1:2045833/1‑150 (MQ=255)
                 tACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACt                                                                                                                                     >  1:1043969/1‑150 (MQ=255)
                 tACAAGAACCAGAGCTAGGGTGTGCTGCGCAGCGCAGACTGGGCGGAGGGCCAGGTGGAAGTGCGTGAAGGCAGAACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACt                                                                                                                                     >  2:4577717/1‑150 (MQ=255)
                  aCAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTg                                                                                                                                    >  1:2812653/1‑150 (MQ=255)
                   cAAGAACCAGAGCACGGGTGAGCTGCGCAGCGCAGCCTGGGAGGAGGGCCTCGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGCCGACATCGTGCCATACATCCCCTTCAGCGAGCTACTAGGCTGAGCGGCCCACCCAGGCACAgg                                                                                                                                   >  2:2171985/1‑150 (MQ=255)
                    aaGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGc                                                                                                                                  >  2:201694/1‑150 (MQ=255)
                      gAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCa                                                                                                                                >  1:5047339/1‑150 (MQ=255)
                       aaCCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGTGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAg                                                                                                                               >  1:3683514/1‑150 (MQ=255)
                         ccAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAgcg                                                                                                                             >  1:3773865/1‑150 (MQ=255)
                         ccAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAgcg                                                                                                                             <  1:2936412/150‑1 (MQ=255)
                           agagCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAgcggc                                                                                                                           >  2:2468497/1‑150 (MQ=255)
                           agagCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAgcggc                                                                                                                           >  1:625862/1‑150 (MQ=255)
                                  gggTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCt                                                                                                                    >  1:3603259/1‑150 (MQ=255)
                                    gtgtGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCg                                                                                                                  >  1:4482985/1‑150 (MQ=255)
                                       tgctgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTcc                                                                                                               <  1:1234049/150‑1 (MQ=255)
                                        gctgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGccagtgcct                                                                                                              <  1:1608918/150‑8 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:2191559/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:777826/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:1281659/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:3484843/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:3350865/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  2:2386301/1‑150 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            <  2:3271208/150‑1 (MQ=255)
                                          tgcGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGt                                                                                                            >  1:3208070/1‑150 (MQ=255)
                                            cgcagcgcagCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGGGCCGCtcagccaa                                                                                                                                          >  2:4893745/1‑117 (MQ=255)
                                             gcagcgcagcCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTgcc                                                                                                         >  2:3115740/1‑150 (MQ=255)
                                               agcgcagcCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTgccgc                                                                                                       >  2:1188004/1‑150 (MQ=255)
                                                gcgcagcCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTgccgct                                                                                                      <  2:300801/150‑2 (MQ=255)
                                                    agcCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTggg                                                                                                  <  2:2819346/150‑1 (MQ=255)
                                                     gcCTGGGCGGAGGGCCTGGTGGTAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGt                                                                                                 <  2:3803815/150‑1 (MQ=255)
                                                        tGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAg                                                                                              >  2:3089688/1‑150 (MQ=255)
                                                        tGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAg                                                                                              <  1:4536058/150‑1 (MQ=255)
                                                          ggCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGc                                                                                            <  1:4110240/150‑1 (MQ=255)
                                                                 ggCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTAtgtg                                                                                     <  2:214330/150‑1 (MQ=255)
                                                                 ggCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTAtgtg                                                                                     >  1:1776614/1‑150 (MQ=255)
                                                                    cTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGcac                                                                                  <  1:4288953/150‑1 (MQ=255)
                                                                    cTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGcac                                                                                  >  1:2956246/1‑150 (MQ=255)
                                                                    cTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGcac                                                                                  >  2:2281558/1‑150 (MQ=255)
                                                                     tggtggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCACCTGGGTCAGGCGTATGTGcaca                                                                                 >  1:2188284/1‑150 (MQ=255)
                                                                       gtggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACAtt                                                                               >  2:2214149/1‑150 (MQ=255)
                                                                        tggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATtc                                                                              <  1:1652448/150‑1 (MQ=255)
                                                                        tggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATtc                                                                              <  1:659631/150‑1 (MQ=255)
                                                                        tggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATtc                                                                              <  2:856355/150‑1 (MQ=255)
                                                                        tggAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATtc                                                                              <  2:508171/150‑1 (MQ=255)
                                                                             gTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAg                                                                         <  1:2171985/150‑1 (MQ=255)
                                                                               gCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAgtt                                                                       <  1:4850839/150‑1 (MQ=255)
                                                                                  tGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAgttgtt                                                                    <  2:449778/150‑1 (MQ=255)
                                                                                   gAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAgttgtt                                                                    >  1:414755/1‑149 (MQ=255)
                                                                                   gAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAgttgtt                                                                    <  2:414755/149‑1 (MQ=255)
                                                                                    aaGGCAGCTCGCCGAAGCAGGGCGACTTCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGTGCGGCACAGCCTGGCACTGGCCAGCGTCGCGGCCTCGTCCGTTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGttgttgt                                                                  <  1:3655007/150‑2 (MQ=255)
                                                                                    aaGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTgg                                                                  <  2:1786303/150‑1 (MQ=255)
                                                                                    aaGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTgg                                                                  >  1:4069082/1‑150 (MQ=255)
                                                                                    aaGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTgg                                                                  <  2:2685098/150‑1 (MQ=255)
                                                                                    aaGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTgg                                                                  <  1:2967598/150‑1 (MQ=255)
                                                                                    aaGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTgg                                                                  <  1:1457558/150‑1 (MQ=255)
                                                                                          gcaCGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGcc                                                            >  1:1664137/1‑150 (MQ=255)
                                                                                                cGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTc                                                      >  1:1980372/1‑150 (MQ=255)
                                                                                                cGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTc                                                      >  1:287784/1‑150 (MQ=255)
                                                                                                 gAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTcc                                                     <  2:1664137/150‑1 (MQ=255)
                                                                                                       gggCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGt                                               >  2:1078129/1‑150 (MQ=255)
                                                                                                       gggCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGt                                               <  1:107177/150‑1 (MQ=255)
                                                                                                       gggCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGt                                               >  2:698982/1‑150 (MQ=255)
                                                                                                       gggCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGt                                               <  1:2290527/150‑1 (MQ=255)
                                                                                                       gggCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGt                                               <  2:1497302/150‑1 (MQ=255)
                                                                                                          cGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGgc                                            <  2:4147902/150‑1 (MQ=255)
                                                                                                          cGACATCGTGCCATTCATCCCCTACAGCGAGCTACTGGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGgc                                            <  1:274461/150‑1 (MQ=255)
                                                                                                           gACATCGTGCCATACATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGgcg                                           >  2:2029900/1‑150 (MQ=255)
                                                                                                              aTCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTg                                        <  1:3363015/150‑1 (MQ=255)
                                                                                                                 gTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCg                                     <  1:2319826/150‑1 (MQ=255)
                                                                                                                  tGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCgg                                    >  2:694394/1‑150 (MQ=255)
                                                                                                                   gCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCggg                                   <  2:2391726/150‑1 (MQ=255)
                                                                                                                          aTCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTg                            >  2:2794195/1‑150 (MQ=255)
                                                                                                                          aTCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTg                            <  2:1559552/150‑1 (MQ=255)
                                                                                                                           tCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTgg                           >  1:1093442/1‑150 (MQ=255)
                                                                                                                               cTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTggg                                                                 >  2:1064615/1‑108 (MQ=255)
                                                                                                                                ttCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACAGTCTCCAGTGGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGccc                      <  2:3733276/150‑1 (MQ=255)
                                                                                                                                 tCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCg                     <  2:1724923/150‑1 (MQ=255)
                                                                                                                                   aGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGc                   >  1:2419321/1‑150 (MQ=255)
                                                                                                                                   aGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGc                   >  2:5017109/1‑150 (MQ=255)
                                                                                                                                   aGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGc                   >  2:5017162/1‑150 (MQ=255)
                                                                                                                                         cTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCt             <  1:1835891/150‑1 (MQ=255)
                                                                                                                                           aCTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCa           >  1:1013476/1‑150 (MQ=255)
                                                                                                                                             ttGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAAc         >  2:950053/1‑150 (MQ=255)
                                                                                                                                             ttGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAAc         >  2:614336/1‑150 (MQ=255)
                                                                                                                                             ctGGCTGGGCCGCTCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAAc         >  2:4295363/2‑150 (MQ=255)
                                                                                                                                                 cTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAACTGGc     <  1:2479665/150‑1 (MQ=255)
                                                                                                                                                  tGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAACTGGCt    <  2:265264/150‑1 (MQ=255)
                                                                                                                                                    aGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAACTGGCTCa  >  1:3648822/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGGGCCGCGTGCGCATCTACAAGAACCAGAGCTCGGGTGTGCTGCGCAGCGCAGCCTGGGCGGAGGGCCTGGTGGAAGTGCGTGAAGGCAGCACGCCGAAGCAGGGCGACATCGTGCCATTCATCCCCTTCAGCGAGCTACTTGGCTGAGCGGCCCAGCCAGGCACTGGCCAGCGGCGCGGCCTCGTCCGGTGCCGCCTGGGTCAGGCGTATGTGCACATTCTCCAGTTGTTGGGCAGCCACGCTCCAGTCGTGGCGCTGGCGGGCAAACTGGCGCCCGGCATCGCTCAACTGGCTCA  >  NC_002947/2424756‑2425053

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