Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 495,027 T→G 8.6% intergenic (+17/‑93) recR → / → htpG DNA repair protein RecR/chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913495,0270TG 8.6% 436.4 / 29.6 155intergenic (+17/‑93)recR/htpGDNA repair protein RecR/chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base T (66/74);  new base G (6/7);  total (73/82)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

Only 100 of 155 total aligned reads displayed.

GCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCT  >  NC_000913/494881‑495175
                                                                                                                                                  |                                                                                                                                                    
gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                   <  1:1812447/150‑1 (MQ=255)
  tttGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCt                                                                                                                                                 >  1:1174430/1‑150 (MQ=255)
  tttGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGGCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCt                                                                                                                                                 >  1:996490/1‑150 (MQ=255)
      cgcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGctc                                                                                                                                             <  1:2075246/150‑1 (MQ=255)
       gcgcAATATGACTTGGAAGCCTGCCGAATCGCTCTTGGGGTTCCTTTTGGCGTAGAGATTGAAATTGTCGACTGCACCACTTTGTCACACTCCCTTGACGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGTATCACCTGctat                                                                                                                                            <  2:1930742/150‑3 (MQ=255)
       gcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGctct                                                                                                                                            >  2:1262443/1‑150 (MQ=255)
          cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGc                                                                                                                                         <  1:251262/150‑1 (MQ=255)
           aaTATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCgtt                                                                                                                                        <  2:311130/150‑3 (MQ=255)
             tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTg                                                                                                                                      <  2:195835/150‑1 (MQ=255)
               tGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaa                                                                                                                                    <  2:1191214/150‑1 (MQ=255)
               tGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaa                                                                                                                                    >  1:2359097/1‑150 (MQ=255)
                gACGTGGAAGCCAGCCGAATCGCTCATTGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                   <  2:387881/150‑1 (MQ=255)
                gACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgctt                                                                                                                                   <  1:90918/150‑4 (MQ=255)
                gACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                   <  1:2386776/150‑1 (MQ=255)
                 aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACTAGAGCAGGATCACCTGCTCTCGCTTGAAat                                                                                                                                  <  2:1006782/150‑1 (MQ=255)
                 aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAat                                                                                                                                  <  2:40006/150‑1 (MQ=255)
                  cgaggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatt                                                                                                                                 >  1:2078011/4‑150 (MQ=255)
                   gTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGGCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgcttaaa                                                                                                                                <  1:1164510/150‑7 (MQ=38)
                    tGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                               <  1:2449224/150‑1 (MQ=255)
                    tGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                               >  1:2461708/1‑150 (MQ=255)
                        agccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATtctc                                                                                                                           >  1:705137/1‑150 (MQ=255)
                         gccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTcc                                                                                                                          >  2:947136/1‑150 (MQ=255)
                          ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                         >  1:239185/1‑150 (MQ=255)
                          ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                         <  2:1265577/150‑1 (MQ=255)
                           cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGACATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCt                                                                                                                        >  1:2294529/1‑150 (MQ=255)
                             gccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTg                                                                                                                      >  2:583252/1‑150 (MQ=255)
                                   tCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCa                                                                                                                <  2:1614139/150‑1 (MQ=255)
                                   tCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCa                                                                                                                <  1:1170222/150‑1 (MQ=255)
                                     gCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                              <  2:1878515/150‑1 (MQ=255)
                                     gCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                              >  2:880975/1‑150 (MQ=255)
                                         aTGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctctc                                                                                                          <  2:136997/150‑1 (MQ=255)
                                           ggCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATGAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTccc                                                                                                        >  1:1593280/1‑150 (MQ=255)
                                              gTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCcaca                                                                                                     <  1:665556/150‑1 (MQ=255)
                                               ttCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACAt                                                                                                    >  1:1639158/1‑150 (MQ=255)
                                                tCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATc                                                                                                   <  1:2134097/150‑1 (MQ=255)
                                                         cggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCCCACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTa                                                                                          <  1:342810/150‑1 (MQ=255)
                                                         cggcgAGCGGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTa                                                                                          <  2:1314200/150‑1 (MQ=255)
                                                            cgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAcc                                                                                       <  1:1081835/150‑1 (MQ=255)
                                                               gCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTa                                                                                    >  1:1953385/1‑150 (MQ=255)
                                                                 tGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGCTCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaa                                                                                  >  1:1031529/1‑150 (MQ=255)
                                                                 tGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaa                                                                                  >  1:2323736/1‑150 (MQ=255)
                                                                   gAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAAt                                                                                <  2:1446713/150‑1 (MQ=255)
                                                                       tGGTCGACTGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                            <  2:1776124/150‑1 (MQ=255)
                                                                       tGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCa                                                                            >  1:1009032/1‑150 (MQ=255)
                                                                        ggTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCat                                                                           >  1:1864886/1‑150 (MQ=255)
                                                                         gTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCatt                                                                          >  1:2118724/1‑150 (MQ=255)
                                                                            ccccccaccacGCCCTCACACTCCCTTGCCGGGCCTCATAACATTCGTGTTTACCCCAACCCGACCAGGCCCACCTGCTCTCGCTTGACATTATTCTCCCTTCTCCCCATCTCTCCCACATCCTGATTTTAACCTTAAAATGGCattatt                                                                       <  2:1079358/145‑1 (MQ=255)
                                                                             aCGGCACCACGTTGTCACACTCCCTTGCCTGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTg                                                                      <  2:256895/150‑1 (MQ=255)
                                                                              cGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGa                                                                     <  1:1762282/150‑1 (MQ=255)
                                                                              cGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGa                                                                     <  2:713533/150‑1 (MQ=255)
                                                                              cGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGa                                                                     <  2:705137/150‑1 (MQ=255)
                                                                                gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                   >  1:1101015/1‑150 (MQ=255)
                                                                                gCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                   >  1:41360/1‑150 (MQ=255)
                                                                                gCACCACGTTGTCACAATCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAgg                                                                   <  2:442284/150‑1 (MQ=255)
                                                                                      cGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGAcc                                                             <  2:1850101/150‑1 (MQ=255)
                                                                                            cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                       <  2:2399607/150‑1 (MQ=255)
                                                                                            cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                       <  2:1181635/150‑1 (MQ=255)
                                                                                            cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTCTCCCCATCTCTCACACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                       >  1:329470/1‑150 (MQ=255)
                                                                                              cacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaa                                                     >  1:146662/1‑150 (MQ=255)
                                                                                              cacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaa                                                     <  2:391444/150‑1 (MQ=255)
                                                                                                       gCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTACCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGa                                            >  1:2173003/1‑150 (MQ=255)
                                                                                                        ccGGGCGTCATAAGATTCAATTTTAAGCAAACGAGAGAAGGAAAACCTGATCTCGCTTTAAATTAATATAACTTGTCCACATCTCTACCACATCCTGTTTTTAACCTTAAAATGGCATTATTGATGTAGACCTAAATCAAAGGAAAAGaa                                           >  2:704678/1‑150 (MQ=255)
                                                                                                         cGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTCCTCTAGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACCTCCTGTTATTAACCTTAAAATGGCATTATTGAGGTAGAACTACATGAAAGGACAAGaaa                                          >  2:2123756/1‑150 (MQ=255)
                                                                                                             cGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCg                                      <  2:146662/150‑1 (MQ=255)
                                                                                                               tCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                    >  2:72187/1‑150 (MQ=255)
                                                                                                               tCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                    >  1:685377/1‑150 (MQ=255)
                                                                                                               tCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                    <  2:1831102/150‑1 (MQ=255)
                                                                                                               tCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTg                                    <  1:1260247/150‑1 (MQ=255)
                                                                                                                cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTgg                                   >  1:1861477/1‑150 (MQ=255)
                                                                                                                 aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                  >  1:272562/1‑150 (MQ=255)
                                                                                                                 aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                  >  1:655002/1‑150 (MQ=255)
                                                                                                                     gATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtc                              <  2:2212962/150‑1 (MQ=255)
                                                                                                                      aTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                             >  1:1509291/1‑150 (MQ=255)
                                                                                                                      aTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                             >  2:1966181/1‑150 (MQ=255)
                                                                                                                       ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                            >  2:236285/1‑150 (MQ=255)
                                                                                                                       ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                            >  1:1794565/1‑150 (MQ=255)
                                                                                                                         cGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtc                          >  2:568722/1‑150 (MQ=255)
                                                                                                                         cGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtc                          <  2:2016501/150‑1 (MQ=255)
                                                                                                                         cGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtc                          >  2:168711/1‑150 (MQ=255)
                                                                                                                           tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtcag                        >  2:25203/1‑150 (MQ=255)
                                                                                                                              aatttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACATACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAg                     >  2:1886029/6‑150 (MQ=38)
                                                                                                                               atttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                    >  2:1808755/5‑150 (MQ=39)
                                                                                                                               atttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                    >  2:1413576/5‑150 (MQ=39)
                                                                                                                                aaGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTg                   >  1:576409/1‑150 (MQ=255)
                                                                                                                                 ttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                  >  1:319629/3‑150 (MQ=255)
                                                                                                                                 aGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAACTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                  >  1:1839344/1‑150 (MQ=255)
                                                                                                                                  gCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGaa                 >  1:1422000/1‑150 (MQ=255)
                                                                                                                                    aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAAc               <  1:1597288/150‑1 (MQ=255)
                                                                                                                                     agcGAGAGCAGGTGATCCTGCTCTTGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                            <  2:879657/134‑1 (MQ=255)
                                                                                                                                     agcGAGAGCAGGTGATCCTGCTCTTGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                            >  1:879657/3‑136 (MQ=255)
                                                                                                                                     aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGCCCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACa              >  1:986139/1‑150 (MQ=255)
                                                                                                                                            gCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTg       >  1:1743663/1‑150 (MQ=255)
                                                                                                                                             cAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGc      <  1:1659063/150‑1 (MQ=255)
                                                                                                                                              aGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCa     <  2:310472/150‑1 (MQ=255)
                                                                                                                                              aGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCa     >  1:1985972/1‑150 (MQ=255)
                                                                                                                                              aGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCa     <  2:1173276/150‑1 (MQ=255)
                                                                                                                                               ggATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagt                           <  2:2297038/127‑1 (MQ=255)
                                                                                                                                               ggATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagt                           >  1:2297038/1‑127 (MQ=255)
                                                                                                                                               ggATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAc    >  2:47458/1‑150 (MQ=255)
                                                                                                                                                 aTCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCAccg  >  2:2047166/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCT  >  NC_000913/494881‑495175

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