Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 863,903 T→A 5.8% intergenic (+17/+19) hemH → / ← uraA ferrochelatase/uracil permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947863,9030TA 5.8% 368.7 / 8.3 113intergenic (+17/+19)hemH/uraAferrochelatase/uracil permease
Reads supporting (aligned to +/- strand):  ref base T (50/56);  new base A (4/3);  total (54/59)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 122 total aligned reads displayed.

GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATCA  >  NC_002947/863754‑864052
                                                                                                                                                     |                                                                                                                                                     
gCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCttt                                                                                                                                                       <  1:2760824/150‑1 (MQ=255)
gCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCttt                                                                                                                                                       >  1:4040707/1‑150 (MQ=255)
gCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCttt                                                                                                                                                       >  1:2543047/1‑150 (MQ=255)
  tGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGc                                                                                                                                                     <  1:1388518/150‑1 (MQ=255)
  tGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGc                                                                                                                                                     <  2:5100655/150‑1 (MQ=255)
    gAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCgg                                                                                                                                                   <  1:3735593/150‑1 (MQ=255)
     aaGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGc                                                                                                                                                  <  2:356132/150‑1 (MQ=255)
      agagATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTTAATGATCACCCGGAGTGGGTGATGGTGTTGGCTTATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGcc                                                                                                                                                 <  2:3942836/150‑1 (MQ=255)
        agaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGcccc                                                                                                                                               <  2:2507710/150‑1 (MQ=255)
        agaTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGcccc                                                                                                                                               <  1:413717/150‑1 (MQ=255)
          aTAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCAAAGCGGCCCCtt                                                                                                                                             <  2:4380162/150‑1 (MQ=255)
             ggCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCg                                                                                                                                          >  2:3228673/1‑150 (MQ=255)
              gCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGCGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGa                                                                                                                                         <  2:203851/150‑1 (MQ=255)
                 tGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCg                                                                                                                                      <  1:815078/150‑1 (MQ=255)
                  gCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCgg                                                                                                                                     >  2:577378/1‑150 (MQ=255)
                   cGTGGCAGCGTGCTGTTTGTCGAGGCGGGCGGCCTGGAGCTGGTGCTGGTGCCGTGCCTTAATTATCACCCGGAGTGGGTGTGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGt                                                                                                                                    <  2:2277486/150‑1 (MQ=255)
                    gTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGACGCTTTGCGGCCCCTTCCGACCGGTc                                                                                                                                   >  2:1200711/1‑150 (MQ=255)
                             aGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCa                                                                                                                          <  2:4688462/150‑1 (MQ=255)
                             aGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCa                                                                                                                          <  2:3105303/150‑1 (MQ=255)
                              gCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGACGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAg                                                                                                                         <  1:3103488/150‑1 (MQ=255)
                              gCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAg                                                                                                                         <  2:553669/150‑1 (MQ=255)
                              gCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCCCGTTGAGGCCCCTTCCGACCGGTAAAACATCCac                                                                                                                         >  1:3948305/1‑149 (MQ=255)
                                 gTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGctt                                                                                                                      >  2:2744791/1‑150 (MQ=255)
                                  tttGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGcttc                                                                                                                     <  2:4647407/150‑1 (MQ=255)
                                    tGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGcttctt                                                                                                                   <  2:875071/150‑1 (MQ=255)
                                     gTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtg                                                                                                                  <  1:817728/150‑1 (MQ=255)
                                      tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATTATCACCCGGAGTGGGTGTGGGTGTTGGCTGATATTTGTGAAAAGGCCTGAATCGACGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgt                                                                                                                 <  2:115385/150‑1 (MQ=255)
                                      tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgt                                                                                                                 <  2:1714761/150‑1 (MQ=255)
                                      tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgt                                                                                                                 <  2:2115081/150‑1 (MQ=255)
                                      tCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTtgt                                                                                                                 <  2:115397/150‑1 (MQ=255)
                                        gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGTTTTTGGCTGATATGTGTCAAAAGGCCTGAAACGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGc                                                                                                               <  1:1573462/150‑1 (MQ=255)
                                        gAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCgg                                                                                                                                     <  2:3320404/128‑1 (MQ=255)
                                                cggcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCt                                                                                                       <  2:2113326/150‑1 (MQ=255)
                                                  gcCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGt                                                                                                     >  1:2347645/1‑150 (MQ=255)
                                                   ccAGGAGCTGGTGCTGGTGCCGTGCCTGATTGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTc                                                                                                    <  1:3816829/150‑1 (MQ=255)
                                                    cAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCa                                                                                                   >  2:501293/1‑150 (MQ=255)
                                                     aGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCGTTGCGGCCCCATACGACCGCTCATACATCCAGCTTCTTGTGCTTACACCGGTCAt                                                                                                  >  1:4427529/1‑150 (MQ=255)
                                                         gctggtgctggtGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAcc                                                                                              <  2:1522571/150‑1 (MQ=255)
                                                         gctggtgctggtGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTAcc                                                                                              >  1:1522571/1‑150 (MQ=255)
                                                            ggtgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATACAGCTTCTTGTGCTTCCAGCTGTCATTACCCgg                                                                                           >  1:2383216/1‑150 (MQ=255)
                                                             gtgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGc                                                                                          >  2:2018930/1‑150 (MQ=255)
                                                              tgctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCa                                                                                         <  2:3723950/150‑1 (MQ=255)
                                                               gctggtgcCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCGTCCGACCGGTCAAACATACAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAg                                                                                        >  1:1753133/1‑150 (MQ=255)
                                                                   gtgccgtgccTGAATGATCACCCGGTGTGGGTGTGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATc                                                                                    <  1:3505447/150‑1 (MQ=255)
                                                                    ggccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGACGGGGCCGATTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCa                                                                                   <  1:98718/149‑1 (MQ=255)
                                                                      ccgtgccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAgg                                                                                 >  1:2743959/1‑150 (MQ=255)
                                                                          gccTGAATGATCACCCGGTGTGGGTGTGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCa                                                                             <  1:3210817/150‑1 (MQ=255)
                                                                           ccTGACTGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAg                                                                            <  2:204960/150‑1 (MQ=255)
                                                                           ccTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTTCGGCCCCTTCCGACCGGTCATACATCCAGCGTCTTTTGCTTCCAGCTGTCATTAACCGGCAGAATCAGGTTCAg                                                                            >  1:2575962/1‑150 (MQ=255)
                                                                            cTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAgg                                                                           <  2:36085/150‑1 (MQ=255)
                                                                             tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGACGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAggg                                                                          <  1:2420316/150‑1 (MQ=255)
                                                                             tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAggg                                                                          >  1:126857/1‑150 (MQ=255)
                                                                             tGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATAAATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATAAGGTTCAggg                                                                          >  1:2679716/1‑150 (MQ=255)
                                                                                 tGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAAt                                                                      >  1:5104844/1‑150 (MQ=255)
                                                                                   aTCACCCGGCGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGACGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATat                                                                    <  1:3641299/150‑2 (MQ=255)
                                                                                    tCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGc                                                                   >  1:4567055/1‑150 (MQ=255)
                                                                                       cccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCcac                                                                >  2:1045764/1‑150 (MQ=255)
                                                                                       cccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCcac                                                                >  2:2360921/1‑150 (MQ=255)
                                                                                        ccGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCcaca                                                               >  2:2210429/1‑150 (MQ=255)
                                                                                         cGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCATTGAGGCCCCTTACGACCGGTCATACATCCAGCTTCTTGTGATTCCAGCTGTCAATACCCGGCAGAATCAGGTTCAGGGCAATAGCCACaa                                                              >  2:630439/1‑150 (MQ=255)
                                                                                            aGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATgg                                                           <  2:5021579/150‑1 (MQ=255)
                                                                                              tGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGgcg                                                         <  2:4081302/150‑1 (MQ=255)
                                                                                                  tGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGcaca                                                     <  2:273626/150‑1 (MQ=255)
                                                                                                       gtgtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCg                                                >  2:4287085/1‑150 (MQ=255)
                                                                                                         gtTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGAt                                              <  1:3218213/150‑1 (MQ=255)
                                                                                                         gtTGGCTGATATGTGTGAAAAGGCCCGAAGCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGAt                                              <  1:1555169/150‑1 (MQ=255)
                                                                                                            ggCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCAAACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTAAGGGCAATAGCCACAATGGCGCACAGGGCGATGcc                                           >  2:552608/1‑150 (MQ=255)
                                                                                                             gCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACAACCAGCTTCTTGTGCTTCCAGCTGTCATTAACCCGCAGAAACAGGTTAAGGGCAATAGCCACAAAGGCGCACAGGGCGATGccc                                          >  1:736545/1‑150 (MQ=255)
                                                                                                               tGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCtt                                        >  2:2702718/1‑150 (MQ=255)
                                                                                                                 atatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCa                                      <  2:4395766/150‑1 (MQ=255)
                                                                                                                 atatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCa                                      <  2:4395576/150‑1 (MQ=255)
                                                                                                                  tatGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAg                                     <  1:4515257/150‑1 (MQ=255)
                                                                                                                     gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcc                                  >  2:2385809/1‑150 (MQ=255)
                                                                                                                     gtgtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCATCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcc                                  >  1:3327869/1‑150 (MQ=255)
                                                                                                                      tgtgAAAAGGCCTGAATCGAAGGGGCCGCGTTTCGGCCCCATCCGACCGGGCATACATCCAGCGTCTTTTGCTTCCAGCTGTCATTAACCGGCAGAATCAGGGTCAGGGCCATAGCCACAAATGCGCACAGGGCGATGCCCATCAGGccc                                 >  1:1447438/1‑150 (MQ=255)
                                                                                                                       gtgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGcccc                                <  2:4040707/150‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTTAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                               <  2:2616733/150‑1 (MQ=255)
                                                                                                                        tgAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCa                               <  2:3202523/150‑1 (MQ=255)
                                                                                                                           aaaGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTc                            >  2:4232527/1‑150 (MQ=255)
                                                                                                                             aGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTTTGCTTCCAGCTGTCATTAACCGGCAGAATCAGGTTAAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCAt                          >  1:4476997/1‑150 (MQ=255)
                                                                                                                               gCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAAAGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATcc                        >  1:4614969/1‑150 (MQ=255)
                                                                                                                                ccTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCg                       >  2:3581397/1‑150 (MQ=255)
                                                                                                                                ccTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCg                       >  2:4584598/1‑150 (MQ=255)
                                                                                                                                 cTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCAACCgg                      >  2:3207755/1‑150 (MQ=255)
                                                                                                                                 cTGAATCGAAGGGGCCGCTTTGCGGCCCCATCCGCCCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCCGCAGAATCAGGTTCAGGGCCATAGCCACAAAGGCGCACAGGGCGATGCCCGTCAGGCCCCAGTCATCCgg                      >  2:616527/1‑150 (MQ=255)
                                                                                                                                   gAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGc                    >  2:2302108/1‑150 (MQ=255)
                                                                                                                                         aaGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAc              >  1:2598064/1‑150 (MQ=255)
                                                                                                                                         aaGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAc              <  2:2593310/150‑1 (MQ=255)
                                                                                                                                         aaGGGGCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCCGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCgagtc              >  2:3838826/1‑146 (MQ=255)
                                                                                                                                          aGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCAcc             <  2:2265792/150‑1 (MQ=255)
                                                                                                                                           ggggCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCg            >  1:1160325/1‑150 (MQ=255)
                                                                                                                                           ggggCCGCAAAGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCg            >  1:654667/1‑150 (MQ=255)
                                                                                                                                            gggCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCCGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGc           >  2:1302022/1‑150 (MQ=255)
                                                                                                                                               ccGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGc        <  1:577378/150‑1 (MQ=255)
                                                                                                                                                 gCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCACTAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCg      >  2:46724/1‑150 (MQ=255)
                                                                                                                                                  cTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGa     <  1:389062/150‑1 (MQ=255)
                                                                                                                                                   tttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGAt    <  1:3228673/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATc   <  1:1139597/150‑1 (MQ=255)
                                                                                                                                                    ttGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATc   <  1:1139607/150‑1 (MQ=255)
                                                                                                                                                     tGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATCa  <  1:966073/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTGGAAGAGATAGGCATGCGTGGCAGCGAGCAGTTTGTCGAGGCGGGCGGCCAGGAGCTGGTGCTGGTGCCGTGCCTGAATGATCACCCGGAGTGGGTGAGGGTGTTGGCTGATATGTGTGAAAAGGCCTGAATCGAAGGGGCCGCTTTGCGGCCCCTTCCGACCGGTCATACATCCAGCTTCTTGTGCTTCCAGCTGTCATTACCCGGCAGAATCAGGTTCAGGGCAATAGCCACAATGGCGCACAGGGCGATGCCCTTCAGGCCCCAGTCATCCGGGCCGTCACCGCTGCCGATCA  >  NC_002947/863754‑864052

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