Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,492,276 G→A 6.0% intergenic (+30/+65) cysNC → / ← imm sulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,492,2760GA 6.0% 614.2 / 18.3 198intergenic (+30/+65)cysNC/immsulfate adenylyltransferase subunit 1/adenylyl‑sulfate kinase/colicin‑E7 immunity protein
Reads supporting (aligned to +/- strand):  ref base G (93/93);  new base A (5/7);  total (98/100)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 207 total aligned reads displayed.

AAAGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCATTTTCC  >  NC_002947/1492128‑1492425
                                                                                                                                                    |                                                                                                                                                      
aaaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcgg                                                                                                                                                       >  2:826607/1‑150 (MQ=255)
 aaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCGGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                                      >  1:2599027/1‑150 (MQ=255)
 aaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCGGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                                      >  1:2106609/1‑150 (MQ=255)
 aaGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTcggc                                                                                                                                                      >  2:2342534/1‑150 (MQ=255)
    cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                   <  2:2628447/150‑1 (MQ=255)
    cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                   >  2:3755999/1‑150 (MQ=255)
    cTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCggg                                                                                                                                                   <  1:2404811/150‑1 (MQ=255)
         cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttt                                                                                                                                              >  1:3596873/1‑150 (MQ=255)
         cGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGtttt                                                                                                                                              <  1:3923217/150‑1 (MQ=255)
           ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGAGGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTGCGGCGGGGttgtgt                                                                                                                                            >  2:4425556/1‑146 (MQ=255)
           ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttttt                                                                                                                                            <  2:2486133/150‑1 (MQ=255)
           ttCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGtttttt                                                                                                                                            <  2:5335213/150‑1 (MQ=255)
            tCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                           <  1:3556481/150‑1 (MQ=255)
            tCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                           <  2:4290947/150‑1 (MQ=255)
              gATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTa                                                                                                                                         >  2:676947/1‑150 (MQ=255)
                   gCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCt                                                                                                                                    <  1:348412/150‑1 (MQ=255)
                          gATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGttttccttg                                                                                                                                         >  1:4087445/1‑133 (MQ=255)
                          gATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCCGAAGCGGGGttttccttg                                                                                                                                         <  2:4087445/138‑6 (MQ=255)
                           aTGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCgg                                                                                                                           <  1:634390/150‑1 (MQ=255)
                             gCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGt                                                                                                                         >  1:317667/1‑150 (MQ=255)
                             gCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGt                                                                                                                         >  2:3486903/1‑150 (MQ=255)
                                gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGaa                                                                                                                      <  2:64194/150‑1 (MQ=255)
                                gATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGaa                                                                                                                      <  1:234362/150‑1 (MQ=255)
                                  tCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGaaaa                                                                                                                    <  1:4521006/150‑1 (MQ=255)
                                    ttGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAgg                                                                                                                  <  2:770946/150‑1 (MQ=255)
                                    ttGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGAAAAACCCCGCCGAAGCGGGGttt                                                                                                                                               >  2:379143/1‑122 (MQ=255)
                                     tgtgATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGc                                                                                                                 >  1:1899836/1‑150 (MQ=255)
                                     tgtgATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGc                                                                                                                 <  2:3873892/150‑1 (MQ=255)
                                       tgATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGccc                                                                                                               <  1:3959373/150‑1 (MQ=255)
                                        gATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                              >  2:2556560/1‑150 (MQ=255)
                                        gATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCt                                                                                                              >  1:179417/1‑150 (MQ=255)
                                          tCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCa                                                                                                            <  1:2789969/150‑1 (MQ=255)
                                           cGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCag                                                                                                           >  1:2486640/1‑150 (MQ=255)
                                               acCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGtt                                                                                                       >  2:1272757/1‑150 (MQ=255)
                                                 ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCt                                                                                                                                          <  1:5259305/114‑1 (MQ=255)
                                                 ccAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCa                                                                                                     >  2:2874390/1‑150 (MQ=255)
                                                  cAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAg                                                                                                    >  1:1253898/1‑150 (MQ=255)
                                                    ggCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAggg                                                                                                  >  1:3034661/1‑150 (MQ=255)
                                                     gCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGc                                                                                                 >  1:2018364/1‑150 (MQ=255)
                                                      ccaccaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGcc                                                                                                <  2:4760359/150‑1 (MQ=255)
                                                      ccaccaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGcc                                                                                                <  2:1051244/150‑1 (MQ=255)
                                                        accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGGTGAAAAGGCCCTCAGAGTTCATGGCCt                                                                                               <  2:4129030/150‑1 (MQ=255)
                                                        accaGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCtt                                                                                              >  2:1340366/1‑150 (MQ=255)
                                                              gTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgt                                                                                        <  1:4470673/150‑1 (MQ=255)
                                                              gTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgt                                                                                        <  2:2105491/150‑1 (MQ=255)
                                                               tGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtg                                                                                       <  2:2452404/150‑1 (MQ=255)
                                                                ggACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGTCCTTGCTCgtgt                                                                                      <  2:2746052/150‑1 (MQ=255)
                                                                ggACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtgt                                                                                      <  1:1398597/150‑1 (MQ=255)
                                                                ggACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCgtgt                                                                                      >  2:4825311/1‑150 (MQ=255)
                                                                 gACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGtt                                                                                     >  1:925324/1‑150 (MQ=255)
                                                                   cGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAg                                                                                   <  2:321547/150‑1 (MQ=255)
                                                                    gAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGc                                                                                  <  2:4182815/150‑1 (MQ=255)
                                                                     aaGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                                 >  1:2674648/1‑150 (MQ=255)
                                                                     aaGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCt                                                                                 <  1:4629453/150‑1 (MQ=255)
                                                                      aGGCGTGAAACAGGTGCTGGACGTGATGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTg                                                                                <  1:890965/150‑1 (MQ=255)
                                                                         cGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGc                                                                             <  1:1617203/150‑1 (MQ=255)
                                                                          gTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCt                                                                            >  2:2985788/1‑150 (MQ=255)
                                                                           tGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCtt                                                                           <  1:1223041/150‑1 (MQ=255)
                                                                             aaaCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAACAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGa                                                                         <  2:4218811/150‑1 (MQ=255)
                                                                               aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGaaa                                                                       >  1:2598212/1‑150 (MQ=255)
                                                                               aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGaaa                                                                       <  2:2388135/150‑1 (MQ=255)
                                                                               aCAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGaaa                                                                       <  1:4096278/150‑1 (MQ=255)
                                                                                 aGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAAcc                                                                     <  1:514165/150‑1 (MQ=255)
                                                                                    tGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTc                                                                  <  2:524637/150‑1 (MQ=255)
                                                                                      cTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCtt                                                                <  1:4678756/150‑1 (MQ=255)
                                                                                       tGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCttt                                                               >  1:1542026/1‑150 (MQ=255)
                                                                                        ggACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTa                                                              <  1:2640952/150‑1 (MQ=255)
                                                                                          aCGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTAcc                                                            <  1:2985788/150‑1 (MQ=255)
                                                                                            gTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGt                                                          <  1:724406/150‑1 (MQ=255)
                                                                                             tgctgcGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGtt                                                         >  1:3495065/1‑150 (MQ=255)
                                                                                                     cAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGACTTTACCGTTGGCTAgcc                                                 >  2:3970746/1‑150 (MQ=255)
                                                                                                      aGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAgccg                                                >  1:1452098/1‑150 (MQ=255)
                                                                                                         gTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCa                                             >  1:4755799/1‑150 (MQ=255)
                                                                                                           ggCGCGATCTAAGCGCTAGCCTGCAAGGAAATCCCCGCTTCGGCGGGGTTTTGTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCAtt                                           >  2:3159761/1‑150 (MQ=255)
                                                                                                             cgcgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGc                                         <  2:5299048/150‑1 (MQ=255)
                                                                                                               cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                       <  2:4593342/150‑1 (MQ=255)
                                                                                                               cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                       <  1:5046496/150‑1 (MQ=255)
                                                                                                               cgATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                       <  2:1198975/150‑1 (MQ=255)
                                                                                                                gATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                       <  1:4898622/149‑1 (MQ=255)
                                                                                                                gATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCtt                                       >  2:4898622/1‑149 (MQ=255)
                                                                                                                   cTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                           >  2:2070470/1‑47 (MQ=255)
                                                                                                                   cTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTc                                                                                                                                           <  1:2070501/47‑1 (MQ=255)
                                                                                                                    tAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTa                                  >  1:3792173/1‑150 (MQ=255)
                                                                                                                    tAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTa                                  >  2:3151091/1‑150 (MQ=255)
                                                                                                                     aaGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTAt                                 >  2:5156412/1‑150 (MQ=255)
                                                                                                                       gcgcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTAttt                               <  1:104901/150‑1 (MQ=255)
                                                                                                                        cgcTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTc                              >  2:5414459/1‑150 (MQ=255)
                                                                                                                          cgaGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                            <  2:3580507/148‑1 (MQ=255)
                                                                                                                          cTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTa                            <  1:4425556/150‑1 (MQ=255)
                                                                                                                           tAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGc                                                                                                 <  2:1838568/80‑1 (MQ=255)
                                                                                                                                 gCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttt                     >  2:499349/1‑150 (MQ=255)
                                                                                                                                 gCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACttt                     >  1:1408021/1‑150 (MQ=255)
                                                                                                                                     ggAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCa                 <  1:826607/150‑1 (MQ=255)
                                                                                                                                     ggAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCa                 <  1:398309/150‑1 (MQ=255)
                                                                                                                                      gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCt                                                                                           <  2:4795930/75‑1 (MQ=255)
                                                                                                                                      gAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCt                                                                                           >  1:4795930/1‑75 (MQ=255)
                                                                                                                                         aaCCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGtt             >  1:4995035/1‑150 (MQ=255)
                                                                                                                                              cGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCa        <  1:2136534/150‑1 (MQ=255)
                                                                                                                                                 ttCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCAttt     >  2:458490/1‑150 (MQ=255)
                                                                                                                                                   cggcggGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCATTTTc   <  2:1489903/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
AAAGCTTCCCGTTCGATGTGCCGCTGGATGCTGATCTTGTGATCGACACCCAGGCCACCAGCGTGGACGAAGGCGTGAAACAGGTGCTGGACGTGCTGCGTCAGCGTGGCGCGATCTAAGCGCTAGCCTGCAAGGAAAACCCCGCTTCGGCGGGGTTTTTTCTATTGCTCAGTTCGGGTGAAAAGGCCCTCAGAGTTCAGGGCCTTGCTCGTGTTAGCTGGGCTTGAAACCGTCTTTACCGTTGGCTAGCCGCCATTGCTTGACTATTTCTAGTACTTTTTCAGGTTTACCATTTTCC  >  NC_002947/1492128‑1492425

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: