Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,367,779 G→C 23.4% A249A (GCC→GCG ftsH ← integral membrane ATP‑dependent zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,367,7790GC23.4% 479.2 / 156.9 231A249A (GCC→GCGftsHintegral membrane ATP‑dependent zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (90/87);  new base C (24/30);  total (114/117)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

Only 100 of 235 total aligned reads displayed.

TGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACC  >  NC_002947/5367631‑5367928
                                                                                                                                                    |                                                                                                                                                     
 gacgatgaTGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGc                                                                                                                                                     >  1:4028049/1‑150 (MQ=255)
    gatgatGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTg                                                                                                                                                  <  2:2860684/150‑1 (MQ=255)
          gCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGcttctt                                                                                                                                            >  1:3895399/1‑150 (MQ=255)
           ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                                                                                           <  1:140041/150‑1 (MQ=255)
           ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                                                                                           <  1:3738249/150‑1 (MQ=255)
           ccATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTg                                                                                                                                           <  2:1144216/150‑1 (MQ=255)
            catcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTgg                                                                                                                                          >  1:2358133/1‑150 (MQ=255)
             atcatTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGc                                                                                                                                         <  2:4325222/150‑1 (MQ=255)
                attcatCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTg                                                                                                                                      <  1:4352492/150‑1 (MQ=255)
                   catCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTc                                                                                                                                   <  2:191058/150‑1 (MQ=255)
                    atCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCg                                                                                                                                  >  2:751045/1‑150 (MQ=255)
                     tctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGa                                                                                                                                 >  1:3663537/1‑150 (MQ=255)
                      ctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                                                                                                                <  2:4724830/150‑1 (MQ=255)
                      ctcGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGaa                                                                                                                                <  1:1590878/150‑1 (MQ=255)
                         gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                                                                                             >  2:3519665/1‑150 (MQ=255)
                         gAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAt                                                                                                                             >  1:2484175/1‑150 (MQ=255)
                           aaccatCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGt                                                                                                                           <  1:3450083/145‑1 (MQ=255)
                              catCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcac                                                                                                                        <  2:2259156/148‑3 (MQ=255)
                               atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                                                                                       >  1:5042428/2‑150 (MQ=255)
                               atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                                                                                       >  1:2969049/2‑150 (MQ=255)
                               atCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacg                                                                                                                       >  2:366916/2‑150 (MQ=255)
                                tCCATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAATATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgg                                                                                                                      >  1:3194783/1‑149 (MQ=255)
                                 ccATTTCGACCAGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTcacgga                                                                                                                     <  1:3095643/150‑3 (MQ=255)
                                 ccATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGca                                                                                                                     >  1:59818/1‑150 (MQ=255)
                                   atttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGgacg                                                                                                                   >  1:3948617/7‑149 (MQ=255)
                                     ttcgaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGAcgcg                                                                                                                 <  1:1116798/146‑1 (MQ=255)
                                      tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGACCATGTCGCGCACACgg                                                                                                                >  2:3383527/1‑150 (MQ=255)
                                      tcCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACgg                                                                                                                >  1:1534990/1‑150 (MQ=255)
                                        gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCTGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTg                                                                                                                                           <  1:4939635/120‑1 (MQ=255)
                                        gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAAc                                                                                                                               <  1:857924/132‑1 (MQ=255)
                                        gaccaGAAGCTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAAc                                                                                                                               >  2:857924/2‑133 (MQ=255)
                                        caccaGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                                                                                              <  1:668047/150‑1 (MQ=255)
                                        caccaGCAACTGGGTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACAGGTCGCGCACACGGCt                                                                                                              >  2:366128/1‑150 (MQ=255)
                                           cagaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCTAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGc                                                                                                           >  1:410381/7‑150 (MQ=255)
                                           caGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGc                                                                                                           <  2:2417946/150‑1 (MQ=255)
                                            agaagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCa                                                                                                          <  1:3837570/145‑1 (MQ=255)
                                            aGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGACGCGCACACGGCTGGCa                                                                                                          <  1:499561/150‑1 (MQ=255)
                                             cgctcTGGTTCAGCGTCTGGTCGCGTTCGTCGTTGCCGCCGCCCAGCCCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                                                                                         <  2:5204631/146‑1 (MQ=255)
                                              aagcTGGTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAcc                                                                                                        <  1:2178079/147‑1 (MQ=255)
                                                aCTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGa                                                                                                      <  2:3367626/150‑1 (MQ=255)
                                                aCTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGa                                                                                                      >  1:474135/1‑150 (MQ=255)
                                                   ggTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGacac                                                                                                   >  1:1877979/1‑150 (MQ=255)
                                                    gTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAcc                                                                                                  >  2:2202634/1‑150 (MQ=255)
                                                    gTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAcc                                                                                                  >  2:452025/1‑150 (MQ=255)
                                                    gTTCAGCGTCTGTTCGCGTTCGTCGTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACAcc                                                                                                  >  1:4868992/1‑150 (MQ=255)
                                                     ttGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAccc                                                                                                 >  1:1248016/1‑150 (MQ=255)
                                                     ttGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACAccc                                                                                                 >  2:637096/1‑150 (MQ=255)
                                                       gAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCAc                                                                                               >  2:194270/1‑150 (MQ=255)
                                                            tttGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACa                                                                                          >  2:4605649/1‑150 (MQ=255)
                                                             ttGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACAt                                                                                         <  2:3277994/150‑1 (MQ=255)
                                                                ttcgcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTc                                                                                      <  2:3068023/143‑1 (MQ=255)
                                                                cTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCt                                                                                                              >  1:927690/1‑126 (MQ=255)
                                                                cTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTc                                                                                      >  1:2769028/1‑150 (MQ=255)
                                                                 tCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTcc                                                                                     >  2:5380919/1‑150 (MQ=255)
                                                                  cACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCa                                                                                    <  2:2942133/150‑1 (MQ=255)
                                                                   gcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                                                                                         >  1:303162/5‑128 (MQ=255)
                                                                   gcgttcgtcgTTGCCGCCGCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCAc                                                                                                         <  2:303162/124‑1 (MQ=255)
                                                                   aCGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCAAc                                                                                   >  2:4128347/1‑150 (MQ=255)
                                                                     gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa                                                                                 >  1:4609798/1‑150 (MQ=255)
                                                                     gCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGa                                                                                 >  1:784520/1‑150 (MQ=255)
                                                                      cTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGaa                                                                                <  2:3949227/150‑1 (MQ=255)
                                                                          tcgttgccgccgCCAAGACCGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTca                                                                            <  2:2184548/144‑2 (MQ=255)
                                                                             tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa                                                                         >  2:4964605/1‑150 (MQ=255)
                                                                             tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa                                                                         >  2:806113/1‑150 (MQ=255)
                                                                             tGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGaa                                                                         <  1:4241960/150‑1 (MQ=255)
                                                                                ccgccgCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCg                                                                      <  1:4438259/150‑1 (MQ=255)
                                                                                 cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                     >  1:2008691/1‑150 (MQ=255)
                                                                                 cgccgcCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTGCGTCGAGGACGATGATGCACGGGGCGTGCGGCGTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCgg                                                                     >  2:3450888/1‑150 (MQ=255)
                                                                                      cccATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATg                                                                >  2:2691911/1‑150 (MQ=255)
                                                                                        cATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGt                                                              <  1:2132463/150‑1 (MQ=255)
                                                                                          tGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTga                                                            <  1:4462083/150‑1 (MQ=255)
                                                                                          gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                            <  1:4572227/148‑1 (MQ=255)
                                                                                          gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                            <  1:1759342/148‑1 (MQ=255)
                                                                                          gaccGGCACCACGGTGGCGACCGACGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTga                                                            <  1:3400275/148‑1 (MQ=255)
                                                                                           gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                           <  1:4508804/150‑1 (MQ=255)
                                                                                           gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                           <  1:2147941/150‑1 (MQ=255)
                                                                                           gCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTgaa                                                           <  1:283497/150‑1 (MQ=255)
                                                                                                      gATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACtt                                                >  2:3859653/1‑150 (MQ=255)
                                                                                                             aCCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTttcgcttc                                         <  2:4229784/150‑1 (MQ=255)
                                                                                                                aaCGGCGTCGATTTCGTCGATGAAGATGATGCATGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCGCCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGcc                                      >  1:5176776/1‑150 (MQ=255)
                                                                                                                 aCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGccc                                     >  2:3564164/1‑150 (MQ=255)
                                                                                                                 aCGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGccc                                     <  1:5380919/150‑1 (MQ=255)
                                                                                                                 aCGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGcca                                     <  1:2560415/150‑2 (MQ=255)
                                                                                                                 aCGGCGTCGATTTCGTCGATGAAGATGATGCAAGGCGCATTCTTCTTGGCCTGTTCGAACATGTCACGGACGCGGCTTGCACCAACACCGACAAACATTTCAACAAAGTCAGAACCAGAAATCGTGAAGAACGGCACATTGGCTTCGcca                                     <  2:3769776/150‑2 (MQ=255)
                                                                                                                   ggCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCcgc                                   <  1:5416588/150‑1 (MQ=255)
                                                                                                                    gCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCcgcg                                  >  1:3563486/1‑150 (MQ=255)
                                                                                                                       tCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATg                               <  2:1259437/150‑1 (MQ=255)
                                                                                                                       tCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATg                               <  1:3988993/150‑1 (MQ=255)
                                                                                                                            ttCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCtt                          <  1:4964605/150‑1 (MQ=255)
                                                                                                                               gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                       >  1:68653/1‑150 (MQ=255)
                                                                                                                               gtcgATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGc                       <  2:68653/150‑1 (MQ=255)
                                                                                                                                   aTGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGACCATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcagc                   <  2:4028049/150‑1 (MQ=255)
                                                                                                                                     gAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGACCCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCcagcag                 >  2:1375853/1‑150 (MQ=255)
                                                                                                                                        gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGACCCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt              >  2:1635959/1‑150 (MQ=255)
                                                                                                                                        gatgatGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGt              <  2:2122955/150‑1 (MQ=255)
                                                                                                                                             tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac         >  2:4876099/1‑150 (MQ=255)
                                                                                                                                             tgCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac         >  2:3914652/1‑150 (MQ=255)
                                                                                                                                             tgCACGGGGAGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtac         <  2:1726290/150‑1 (MQ=255)
                                                                                                                                              gCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtacc        >  1:1588760/1‑150 (MQ=255)
                                                                                                                                               cACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTtaccg       >  1:5231861/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
TGACGATGATGCCATCATTCATCTCGAAGCCGTCCATCTCCACCAGCAACTGGTTGAGGGTTTGCTCACGCTCGTCGTGACCGCCGCCCATGCCGGCGCCACGATGGCGACCAACGGCGTCGATTTCGTCGATGAAGATGATGCACGGGGCGTGCTTCTTGGCCTGCTCGAACATGTCGCGCACACGGCTGGCACCGACACCCACGAACATTTCCACGAAGTCCGAACCGGAAATGGTGAAGAACGGTACTTTCGCTTCGCCCGCGATGGCCTTGGCCAGCAGGGTCTTACCGGTACC  >  NC_002947/5367631‑5367928

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: