Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 495,027 T→G 5.0% intergenic (+17/‑93) recR → / → htpG gap repair protein/protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913495,0270TG 5.0% 988.7 / 31.7 278intergenic (+17/‑93)recR/htpGgap repair protein/protein refolding molecular co‑chaperone Hsp90, Hsp70‑dependent; heat‑shock protein; ATPase
Reads supporting (aligned to +/- strand):  ref base T (131/132);  new base G (6/8);  total (138/140)
Fisher's exact test for biased strand distribution p-value = 7.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

Only 100 of 280 total aligned reads displayed.

CGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCT  >  NC_000913/494878‑495175
                                                                                                                                                     |                                                                                                                                                    
cGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGAGTCGTTTTTAAGCAAACGAGAGCAGGAt                                                                                                                                                      >  2:4691887/1‑150 (MQ=255)
 gAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATc                                                                                                                                                     <  2:121203/150‑1 (MQ=255)
   gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                   <  1:4008564/150‑1 (MQ=255)
   gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                   <  2:569945/150‑1 (MQ=255)
   gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                   <  1:2122914/150‑1 (MQ=255)
   gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                   <  1:2338643/150‑1 (MQ=255)
     tttGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCt                                                                                                                                                 <  2:1868722/150‑1 (MQ=255)
       tGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGc                                                                                                                                               <  1:2618339/150‑1 (MQ=255)
       tGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGc                                                                                                                                               <  2:4893858/150‑1 (MQ=255)
        gcgcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGct                                                                                                                                              <  1:1310107/150‑1 (MQ=255)
         cgcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGctc                                                                                                                                             <  1:2781192/150‑1 (MQ=255)
            gcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCg                                                                                                                                          <  2:4492277/150‑1 (MQ=255)
            gcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCg                                                                                                                                          >  1:3624064/1‑150 (MQ=255)
            gcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCg                                                                                                                                          <  1:4950158/150‑1 (MQ=255)
             cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGc                                                                                                                                         <  1:2649306/150‑1 (MQ=255)
             cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGc                                                                                                                                         <  2:2724004/150‑1 (MQ=255)
             cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGc                                                                                                                                         <  1:4646729/150‑1 (MQ=255)
                tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTg                                                                                                                                      <  1:5719537/150‑1 (MQ=255)
                  ggACGTGGAAGCCAGCCGAATCGCTCATGGCGGTCCGGTTGGCGGCGAGCGGGAAATGGTCGACGGCTCCACGTTGTCACACTCCCTTGCGGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaa                                                                                                                                    <  1:4527472/149‑1 (MQ=255)
                   gACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                   <  1:2904147/150‑1 (MQ=255)
                   gACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGaaa                                                                                                                                   <  1:1638401/150‑1 (MQ=255)
                      gTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatta                                                                                                                                <  2:133685/150‑1 (MQ=255)
                        ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattatt                                                                                                                              >  2:2957986/1‑150 (MQ=255)
                        ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattatt                                                                                                                              >  2:4125812/1‑150 (MQ=255)
                         gAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATtc                                                                                                                             <  1:3708653/150‑1 (MQ=255)
                           agccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATtctc                                                                                                                           <  1:5103024/150‑1 (MQ=255)
                            gccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTcc                                                                                                                          <  2:4610196/150‑1 (MQ=255)
                            gccagccGAATCGCTCATGGCGTACCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGAGTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCGCTCGCTGGTAATTATTCTcc                                                                                                                          >  1:3824613/1‑150 (MQ=255)
                             ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                         >  2:368817/1‑150 (MQ=255)
                             ccagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                         <  1:1713118/150‑1 (MQ=255)
                               agccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCtt                                                                                                                       <  2:5442961/150‑1 (MQ=255)
                                  cGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTc                                                                                                                    <  1:4435654/150‑1 (MQ=255)
                                   gAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTcc                                                                                                                   >  1:2202992/1‑150 (MQ=255)
                                   gAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTcc                                                                                                                   >  1:1049704/1‑150 (MQ=255)
                                     aTCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTcccc                                                                                                                 <  2:4354740/150‑1 (MQ=255)
                                      tCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCa                                                                                                                >  2:482561/1‑150 (MQ=255)
                                          tCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctc                                                                                                            >  2:3533368/1‑150 (MQ=255)
                                            aTGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACCTTGTAACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCACTCGCTTGAAATTATTCACCCATTTCCCCAtctctc                                                                                                          >  2:1247671/1‑150 (MQ=255)
                                             tGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTcc                                                                                                         >  1:2558076/1‑150 (MQ=255)
                                                   tCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATc                                                                                                   <  1:2418636/150‑1 (MQ=255)
                                                    ccGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATcc                                                                                                  >  1:1546881/1‑150 (MQ=255)
                                                       gTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGt                                                                                               >  1:5205462/1‑150 (MQ=255)
                                                         tGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttt                                                                                             >  1:5601334/1‑150 (MQ=255)
                                                           gcggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                           >  1:913706/1‑150 (MQ=255)
                                                             ggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTaa                                                                                         <  1:1694017/150‑1 (MQ=255)
                                                                gAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCt                                                                                      >  1:5799204/1‑150 (MQ=255)
                                                                    tGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaa                                                                                  >  2:3911462/1‑150 (MQ=255)
                                                                    tGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaa                                                                                  >  1:2014071/1‑150 (MQ=255)
                                                                    tGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCACTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaa                                                                                  >  1:4189007/1‑150 (MQ=255)
                                                                     ggAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTaaaa                                                                                 <  1:998196/150‑1 (MQ=255)
                                                                      gAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAAt                                                                                >  2:3968284/1‑150 (MQ=255)
                                                                        aaTGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATgg                                                                              <  1:5426146/150‑1 (MQ=255)
                                                                          tGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGGTGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTGCCACAACGTGTTTTTAACCTTAAAATGGCa                                                                            <  1:4897651/150‑1 (MQ=255)
                                                                               gACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCattatt                                                                       <  2:3385920/150‑1 (MQ=255)
                                                                                aCGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTg                                                                      <  1:676999/150‑1 (MQ=255)
                                                                                  ggCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAg                                                                    <  1:2833041/150‑1 (MQ=255)
                                                                                     accacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTa                                                                 >  2:2014464/1‑150 (MQ=255)
                                                                                     accacGGTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTa                                                                 <  2:4226888/150‑1 (MQ=255)
                                                                                       cacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGa                                                               >  2:335265/1‑150 (MQ=255)
                                                                                       cacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGa                                                               >  2:5628520/1‑150 (MQ=255)
                                                                                       cacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGa                                                               <  1:1084459/150‑1 (MQ=255)
                                                                                        acGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGAc                                                              <  1:4458691/150‑1 (MQ=255)
                                                                                                acacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGa                                                      >  2:5284765/1‑150 (MQ=255)
                                                                                                 cacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTCCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGAATTATTGAGGTAGACCTACATGaa                                                     >  1:3446905/1‑150 (MQ=255)
                                                                                                  acTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGaaa                                                    <  1:2957986/150‑1 (MQ=255)
                                                                                                      ccgtGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTAGAAATTATTCTCCCTTGTCCCCATCTCTCCCCCATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGAc                                                <  2:4850217/147‑1 (MQ=255)
                                                                                                      ccTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGAc                                                <  1:3345665/150‑1 (MQ=255)
                                                                                                       cTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACa                                               <  1:1999602/150‑1 (MQ=255)
                                                                                                        ttGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACaa                                              <  2:1794192/150‑1 (MQ=255)
                                                                                                        gtGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAACTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACaa                                              <  1:3172916/149‑1 (MQ=255)
                                                                                                                    aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                  >  2:1142642/1‑150 (MQ=255)
                                                                                                                    aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTAGTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                  >  2:2467893/1‑150 (MQ=255)
                                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                >  2:3641934/1‑150 (MQ=255)
                                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                >  1:5557474/1‑150 (MQ=255)
                                                                                                                      aaGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                >  1:3066513/1‑150 (MQ=255)
                                                                                                                        gATTCGTTTTTAAGCAATCGAGAGCATGATCTCCTGCTCTAGCTTGAAATTATTCGCCCTTGTCCCCATCTCTCCCACATCATGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACATACATGAAATGACAAGAAACTCGTGGtttta                              <  2:3058235/150‑2 (MQ=255)
                                                                                                                         aTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtca                             <  1:1537465/150‑1 (MQ=255)
                                                                                                                           tCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACACGAAACTCGTGGTTTtcagt                           >  1:3977212/1‑150 (MQ=255)
                                                                                                                           tCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagt                           >  1:472777/1‑150 (MQ=255)
                                                                                                                             gTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtca                         >  2:2907766/1‑150 (MQ=255)
                                                                                                                             gTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtca                         >  1:648256/1‑150 (MQ=255)
                                                                                                                              tttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtcag                        >  1:3481520/1‑150 (MQ=255)
                                                                                                                               ttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGa                       >  2:612028/1‑150 (MQ=255)
                                                                                                                               ttttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGa                       >  1:5341887/1‑150 (MQ=255)
                                                                                                                                 ttAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCACCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGGAGACCTACATGAAAGGCCACGACCCTCGTGGTTTTCAGTCAGAAg                     >  1:3276868/1‑150 (MQ=255)
                                                                                                                                   aaGCAAACGAGAGCAGGATCACCTGCTCTCGATTGAAATTATTCTACCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGTCATTATTGATGTAGACCTACATGAAAGTACAAGAAGATCGTGTTTTTCAGTCAGAtgtg                   >  1:3200829/1‑150 (MQ=255)
                                                                                                                                    aGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGGCAGAAGTGa                  >  2:4329162/1‑150 (MQ=255)
                                                                                                                                     gCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGaa                 <  1:2369879/150‑1 (MQ=255)
                                                                                                                                     gCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGaa                 <  1:3503452/150‑1 (MQ=255)
                                                                                                                                      cAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGTAACTCGTGGTTTTAAGTCAGAAGTGaaa                >  2:320622/1‑150 (MQ=255)
                                                                                                                                        agcGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGttt                                <  2:1683450/130‑1 (MQ=255)
                                                                                                                                         gcGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAg             >  2:3743194/2‑150 (MQ=255)
                                                                                                                                         aCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCTGAAGTGAAACAg             >  2:1467342/1‑150 (MQ=255)
                                                                                                                                         aCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAg             >  2:4377520/1‑150 (MQ=255)
                                                                                                                                         aCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAg             >  1:1450415/1‑150 (MQ=255)
                                                                                                                                         aCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAg             >  1:1174164/1‑150 (MQ=255)
                                                                                                                                          cGAGAGCAGGATCACCTGCTCTCGCTTGTAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGc            >  1:610206/1‑150 (MQ=255)
                                                                                                                                                cAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGc      >  1:4685678/1‑150 (MQ=255)
                                                                                                                                                    aTCACCTGCTCTCGCTTGAATTTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCt  <  2:5348583/150‑1 (MQ=255)
                                                                                                                                                    aTCACCGGCTCGCGCTTGAAATTCGTCGCCCTGGTCCCCATCTCTCCCCCATCCTGTTTTGAACCTTAACATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCt  <  2:2361829/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCT  >  NC_000913/494878‑495175

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: