Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,089,094:1 +T 8.6% coding (102/510 nt) PP_2697 ← flavin reductase domain‑containing oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,089,0941.T 8.6% 341.6 / 27.3 93coding (102/510 nt)PP_2697flavin reductase domain‑containing oxidoreductase
Reads supporting (aligned to +/- strand):  ref base . (42/43);  new base T (4/4);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

GACCTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTCCT  >  NC_002947/3088947‑3089243
                                                                                                                                                    |                                                                                                                                                     
gACCTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTATCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTc                                                                                                                                                     >  2:2104443/1‑150 (MQ=255)
   cTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCAcc                                                                                                                                                  >  1:1696239/1‑150 (MQ=255)
   cTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCttcacc                                                                                                                                                 <  2:1981303/150‑6 (MQ=255)
   cTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGCGAGCCGATCGGCttcacc                                                                                                                                                 <  1:806030/150‑6 (MQ=255)
    tGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCTTCACCt                                                                                                                                                 >  1:2305819/1‑150 (MQ=255)
        cACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgt                                                                                                                                             <  1:1316950/150‑1 (MQ=255)
        cACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCTTCACCtgcca                                                                                                                                             <  2:2305819/150‑5 (MQ=255)
        cACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCTTCACCtgcca                                                                                                                                             <  2:666382/150‑5 (MQ=255)
         aCCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgtg                                                                                                                                            >  1:803336/1‑150 (MQ=255)
         aCCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgtg                                                                                                                                            >  1:175047/1‑150 (MQ=255)
            tGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgtggtg                                                                                                                                         >  1:629768/1‑150 (MQ=255)
              aCAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCg                                                                                                                                       >  2:884543/1‑150 (MQ=255)
                     gTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAAt                                                                                                                                >  2:470311/1‑150 (MQ=255)
                        gACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCAc                                                                                                                             >  1:36855/1‑150 (MQ=255)
                         aCCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCAcc                                                                                                                            <  2:784899/150‑1 (MQ=255)
                          ccGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGGCATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCg                                                                                                                           <  1:522567/150‑1 (MQ=255)
                           cGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGt                                                                                                                          <  1:1988043/150‑1 (MQ=255)
                                     ccAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGAt                                                                                                                >  1:2443092/1‑150 (MQ=255)
                                       aGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCATGATGCCGGATGc                                                                                                              <  2:803336/150‑1 (MQ=255)
                                         ccTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGt                                                                                                            <  1:947833/150‑1 (MQ=255)
                                         ccTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGAGGCCTGATGctt                                                                                                            <  2:1473663/150‑3 (MQ=255)
                                         ccTGGCGGATTTGCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGt                                                                                                            <  1:1425946/150‑1 (MQ=255)
                                            ggCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGt                                                                                                         <  2:1311761/150‑1 (MQ=255)
                                             gCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTg                                                                                                        >  1:1348344/1‑150 (MQ=255)
                                                        ggTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCt                                                                                             <  1:1522402/150‑1 (MQ=255)
                                                        ggTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCt                                                                                             >  1:1800370/1‑150 (MQ=255)
                                                         gTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTg                                                                                            <  2:1461994/150‑1 (MQ=255)
                                                                gCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGaaaa                                                                                     <  1:2104443/150‑1 (MQ=255)
                                                                gCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGaaaa                                                                                     >  2:1948109/1‑150 (MQ=255)
                                                                 ccGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGcc                                                                                                                    <  2:914268/118‑1 (MQ=255)
                                                                 ccGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGcc                                                                                                                    >  1:914268/1‑118 (MQ=255)
                                                                   gAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTc                                                                                  <  1:1482003/150‑1 (MQ=255)
                                                                      cTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGAc                                                                               <  2:175047/150‑1 (MQ=255)
                                                                        ggACATCACACTGAATGACACCAGCGGCGGGCTCACCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACgg                                                                             <  1:1948109/150‑1 (MQ=255)
                                                                         gACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGc                                                                            >  1:1704270/1‑150 (MQ=255)
                                                                                cacTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCg                                                                     >  2:245425/1‑150 (MQ=255)
                                                                                 acTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGc                                                                    >  1:2242249/1‑150 (MQ=255)
                                                                                  cTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCGCCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   >  1:627300/1‑150 (MQ=255)
                                                                                  cTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAATTCGACGGCTCGATCGcc                                                                   <  2:1696239/150‑1 (MQ=255)
                                                                                  cTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   <  1:203748/150‑1 (MQ=255)
                                                                                   tGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCg                                                                  >  1:821557/1‑150 (MQ=255)
                                                                                     aaTGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCTTCACCtgc                                                                                                                                               <  1:513119/71‑3 (MQ=37)
                                                                                     aaTGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAGCCGATCGGCTTCACCtgc                                                                                                                                               >  2:513119/1‑69 (MQ=37)
                                                                                      aTGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgtg                                                                                                                                            >  2:2385388/1‑73 (MQ=255)
                                                                                      aTGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCgtg                                                                                                                                            <  1:2385388/73‑1 (MQ=255)
                                                                                           tgcAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGaa                                                          >  1:82549/3‑150 (MQ=255)
                                                                                            ccAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCCCCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGaaa                                                         >  2:2031845/1‑150 (MQ=255)
                                                                                             cAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   <  1:1352659/139‑1 (MQ=255)
                                                                                             cAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   >  2:1352659/1‑139 (MQ=255)
                                                                                                  gcggGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACa                                                   <  1:2453278/150‑1 (MQ=255)
                                                                                                   cggGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACAt                                                  <  1:180280/150‑1 (MQ=255)
                                                                                                    gggCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATg                                                 <  2:1704270/150‑1 (MQ=255)
                                                                                                    gggCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATg                                                 <  2:1931368/150‑1 (MQ=255)
                                                                                                       cTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTa                                              <  1:1069379/150‑1 (MQ=255)
                                                                                                       cTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTa                                              <  2:1800370/150‑1 (MQ=255)
                                                                                                        tCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTg                                                                                            >  1:2391808/1‑103 (MQ=255)
                                                                                                        tCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTg                                                                                            <  2:2391808/103‑1 (MQ=255)
                                                                                                           tCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTc                                          <  1:245425/150‑1 (MQ=255)
                                                                                                               tACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCTAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCtgt                                      <  1:847240/150‑1 (MQ=255)
                                                                                                                 cGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTCCTGGAAAGCAACATGGTAGGTCCCGTgg                                    >  2:1129682/1‑150 (MQ=255)
                                                                                                                  gCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCTGTACTGGAAAGCAACATGGTAGGTCCTGTGGc                                   <  2:1348344/150‑1 (MQ=255)
                                                                                                                  gCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACt                                                              <  1:265544/123‑1 (MQ=255)
                                                                                                                  gCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACt                                                              >  2:265544/1‑123 (MQ=255)
                                                                                                                     gTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGACGTACTGGAAAACAACAGCGAAGGGCCTATggcggc                                >  2:1097202/1‑150 (MQ=255)
                                                                                                                       aGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGcat                              <  2:2443092/150‑1 (MQ=255)
                                                                                                                         aaCGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTTCCTGGCGAAAAGTCGACGGCTCCATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGcattc                            >  2:265572/1‑150 (MQ=255)
                                                                                                                            accacgaaGGTGAGCCGATCGGCTTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGATCCTGGGGCGGcagtcatt                         >  2:1635797/8‑150 (MQ=255)
                                                                                                                              cacgaaGGTGAGCCGATCGGCTTCACCTGCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   >  2:1788109/6‑106 (MQ=38)
                                                                                                                              cacgaaGGTGAGCCGATCGGCTTCACCTGCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGcc                                                                   <  1:1788109/101‑1 (MQ=38)
                                                                                                                               tGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCa                      <  1:2104207/150‑1 (MQ=255)
                                                                                                                               atctaGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGCTGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTct                      <  1:1134495/146‑2 (MQ=255)
                                                                                                                                 gCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAgc                    <  2:1906446/150‑1 (MQ=255)
                                                                                                                                   aGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCtgt                                      >  1:2364425/1‑130 (MQ=255)
                                                                                                                                   aGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCtgt                                      <  2:2364425/130‑1 (MQ=255)
                                                                                                                                   aGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATCTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAgcgt                  >  2:1276162/1‑150 (MQ=255)
                                                                                                                                    ggTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGAGGCGTAGTGCCCGAGCGGCTGGCGAAAAGTCGACGGGTCCATCGCCGTACTGAAAAGCAACATGATAGGAGCTGTGGGTGGAGTCATGTCAgagat                 >  2:1134495/1‑146 (MQ=255)
                                                                                                                                     gTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGCGCGAATATTGAGAAGTCGACGGCGCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAgcgtgc                <  1:884543/150‑1 (MQ=255)
                                                                                                                                     gTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAgcgtgc                <  2:450936/150‑1 (MQ=255)
                                                                                                                                     gTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTACTGTGGCTGCATTCATTTCAgcgtgc                >  1:1950351/1‑150 (MQ=255)
                                                                                                                                       gAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAgcgtgcgt              <  1:1097202/150‑1 (MQ=255)
                                                                                                                                          gCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAg           >  1:1975536/1‑150 (MQ=255)
                                                                                                                                            cGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGaaaa                                                                                     >  2:1960613/1‑74 (MQ=255)
                                                                                                                                            cGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGaaaa                                                                                     <  1:1960613/74‑1 (MQ=255)
                                                                                                                                             gATCGGCTCACCTTCGTGGTGGGATTTAATCGCCTGGATCGCGGAGGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACAGGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGa        <  2:629768/150‑1 (MQ=255)
                                                                                                                                              aTCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAg       >  1:2101589/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGc                                                                                          <  2:2165446/63‑1 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGc                                                                                          >  1:2165446/1‑63 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGc                                                                    >  2:1924859/1‑85 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGc                                                                    <  1:1924859/85‑1 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  2:93363/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  1:1501624/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  2:712775/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  2:463812/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  2:463582/1‑150 (MQ=255)
                                                                                                                                                  gCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTcc   >  2:27039/1‑150 (MQ=255)
                                                                                                                                                   cTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTCCt  >  2:613606/1‑150 (MQ=255)
                                                                                                                                                   cTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCAGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTCCt  >  2:1237231/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACCTGCTCACCTGACAGGATGTTGACCGTGAAGCGGCCAGCCTGGCGGATTTTCGGGTAGCTGGCCGAACTGGACATCACACTGAATGACACCAGCGGCGGGCTCATCGATACGCTGTAGAACGACTGGCAGGTGAAGCCGATCGGCTCACCTTCGTGGTGCGATGTAATCACCGTGATGCCGGATGCGTAGTGCCCGAGTGCCTGGCGAAAAGTCGACGGCTCGATCGCCGTACTGGAAAGCAACATGGTAGGTCCTGTGGCTGCATTCATTTCAGCGTGCGTTAGCGAGTTCCT  >  NC_002947/3088947‑3089243

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