Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 495,027 T→G 8.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 8.0% 300.3 / 13.5 112intergenic (+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 (43/59);  new base G (3/6);  total (46/66)
Fisher's exact test for biased strand distribution p-value = 7.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

Only 100 of 115 total aligned reads displayed.

GAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTG  >  NC_000913/494879‑495176
                                                                                                                                                    |                                                                                                                                                     
gAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATc                                                                                                                                                      <  2:167884/150‑1 (MQ=255)
  gCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAc                                                                                                                                                    <  2:541985/150‑1 (MQ=255)
   cTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAcc                                                                                                                                                   <  1:956959/150‑1 (MQ=255)
    tttGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCAccg                                                                                                                                                  >  1:452476/1‑149 (MQ=255)
       gcgcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCACTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGct                                                                                                                                               >  2:940723/1‑150 (MQ=255)
        cgcgcAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGctc                                                                                                                                              >  1:1480034/1‑150 (MQ=255)
            cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCgt                                                                                                                                          <  1:55154/150‑2 (MQ=255)
            cAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGc                                                                                                                                          <  2:1318836/150‑1 (MQ=255)
              atatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCGttt                                                                                                                                        <  2:1654308/150‑1 (MQ=255)
              atatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGTGATCCTGCTCTCGttt                                                                                                                                        <  1:1226419/150‑1 (MQ=255)
               tatTACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGCGAATGGTCGACGGCACCACGTTGTCACACTCCCTTCCCGGGCGTCATAAGCTTCGTTTTTAAGCAAACGAGAGCAGGACCACCTGCTCTCGCTTg                                                                                                                                       <  1:242828/150‑1 (MQ=255)
               tatGACGTGGACGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTg                                                                                                                                       <  1:121230/150‑1 (MQ=255)
               tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTg                                                                                                                                       <  2:1541327/150‑1 (MQ=255)
               tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTg                                                                                                                                       <  1:2778/150‑1 (MQ=255)
               tatGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGGCGACGGCACCACGTTGTCACACTCCCTTTCCGGGCGTCATAAGATTCGTTTTTAAGCAAACTAGAGCAGGTGATCCTGCTCTCGTTTg                                                                                                                                       <  2:1108841/150‑1 (MQ=255)
                atGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGTAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCTCTCGTTTgc                                                                                                                                      <  2:260796/150‑2 (MQ=255)
                atGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAGCGAGAGCAGGTGATCCTGCtctcttttgc                                                                                                                                      <  2:415305/150‑7 (MQ=255)
                 tGACGTGGAAGCCAGCCGAATCGCTCATGGCGTGCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGTTTgct                                                                                                                                     <  2:1345959/150‑3 (MQ=255)
                   aCGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAat                                                                                                                                   <  1:901567/150‑1 (MQ=255)
                    cGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAatt                                                                                                                                  >  1:1534609/1‑150 (MQ=255)
                      tGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattat                                                                                                                                >  1:839434/1‑150 (MQ=255)
                       ggAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAattatt                                                                                                                               <  2:1560533/150‑1 (MQ=255)
                        gAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATtc                                                                                                                              <  1:1499685/150‑1 (MQ=255)
                            ccagccGAATCGCTCATGGCGTGCCGGTTGGCGGCGAGCTGGAAATGGGCGACGGCACCACGGTGTCACACTCCCTTGCCGGGCGTCATAAGATCCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTccc                                                                                                                          <  1:288363/150‑1 (MQ=255)
                             cagccGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGACGGGCGTCATAAGATTCGTTTTTAAGCAAACCAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCt                                                                                                                         >  1:73160/1‑150 (MQ=255)
                                  gAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTcc                                                                                                                    <  1:933104/150‑1 (MQ=255)
                                       gCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGGCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtc                                                                                                               <  2:165649/150‑1 (MQ=255)
                                         gcATGGCGTTCCGGGTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGGTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCAtctc                                                                                                             <  2:1121351/149‑1 (MQ=255)
                                          cATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCACAtctctc                                                                                                           >  1:46181/1‑150 (MQ=255)
                                             ggCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTccc                                                                                                         >  2:1346652/1‑150 (MQ=255)
                                               cGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCcac                                                                                                       >  1:1196931/1‑150 (MQ=255)
                                                gTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCcaca                                                                                                      <  1:1318835/150‑1 (MQ=255)
                                                  tCCGGTTGGCGGCGAGCTGGAAATGGGCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATc                                                                                                    <  1:460670/150‑1 (MQ=255)
                                                   ccGGTTGGCGGCGAGCTGGAAATGGGCGACGGCGCCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATcc                                                                                                   <  2:1726752/150‑1 (MQ=255)
                                                       ttGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtt                                                                                               <  2:1574182/150‑1 (MQ=255)
                                                         ggcggcGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtttt                                                                                             <  2:82713/150‑1 (MQ=255)
                                                         ggcggcGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGtttt                                                                                             <  2:1580001/150‑1 (MQ=255)
                                                          gcggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                            <  2:1728944/150‑1 (MQ=255)
                                                          gcggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGttttt                                                                                            >  1:388143/1‑150 (MQ=255)
                                                            ggcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCGCCCACATCCTGGTTTTaa                                                                                          >  1:940147/1‑150 (MQ=255)
                                                             gcgAGCTGGAATTGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAc                                                                                         <  1:445077/150‑1 (MQ=255)
                                                             gcgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGGTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCCCATCCTGTTTTTAAc                                                                                         >  2:77356/1‑150 (MQ=255)
                                                             gcgAGCTGGAAATGGTCGAAGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAc                                                                                         <  1:875592/150‑1 (MQ=255)
                                                              cgAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAAATATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAcc                                                                                        >  2:1424756/1‑150 (MQ=255)
                                                                       aaTGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATgg                                                                               <  1:816786/150‑1 (MQ=255)
                                                                       aaTGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATgg                                                                               <  2:29793/150‑1 (MQ=255)
                                                                        aTGGTCGACGGCACGACGTTGTCTCACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAATCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGc                                                                              <  2:1078497/150‑1 (MQ=255)
                                                                          ggTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCGCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCat                                                                            <  2:46305/150‑1 (MQ=255)
                                                                            tCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCatta                                                                          <  2:947332/150‑1 (MQ=255)
                                                                                 ggCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAg                                                                     <  1:372600/150‑1 (MQ=255)
                                                                                   caccacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGt                                                                   >  1:696836/1‑150 (MQ=255)
                                                                                      cacGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGa                                                                <  2:574961/150‑1 (MQ=255)
                                                                                         gttgttACACTCCCTTGCCGGGCGTCATAAGATTCTTTTTTAAGCATACGAGTGCAGGATCAACTGCTCTATCTTGATATTATTCTCGCTTGTCCCCATCTCTCGCACATCCTGTTTTTTACCTTAAATTGGCATTATTGAGGTAGACCt                                                             <  1:1042011/150‑1 (MQ=255)
                                                                                            gTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACa                                                          >  1:1316506/1‑150 (MQ=255)
                                                                                             tCACCCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACAt                                                         >  2:525516/1‑150 (MQ=255)
                                                                                              cacaCTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATg                                                        >  1:718100/1‑150 (MQ=255)
                                                                                               acacTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGa                                                       >  1:1607332/1‑150 (MQ=255)
                                                                                                  cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                    <  2:165374/150‑1 (MQ=255)
                                                                                                  cTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                    <  2:1211132/150‑1 (MQ=255)
                                                                                                  atCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTTTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAg                                                    <  2:696836/149‑1 (MQ=255)
                                                                                                   tCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAgg                                                   <  1:1434073/150‑1 (MQ=255)
                                                                                                     acTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGACAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCACCCACAGCCTGTTTGTACCCTTAAAAAGGCATTATTGACGTAGACATACATGaaactca                                                 >  2:300873/2‑146 (MQ=255)
                                                                                                      cTTGCCGGGCGTCATAATATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCATGTTTTTAAACTTAAAATTGCATTATTGAGTTAGACCTACATGAAAGGaaa                                                >  1:1621916/1‑148 (MQ=255)
                                                                                                       ttGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACaa                                               <  1:1000191/150‑1 (MQ=255)
                                                                                                       ttGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACaa                                               <  2:255747/150‑1 (MQ=255)
                                                                                                        tGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAg                                              <  2:207869/150‑1 (MQ=255)
                                                                                                             ggCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACt                                         <  2:1281083/150‑1 (MQ=255)
                                                                                                                  cATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACCTCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTgg                                    >  1:523174/1‑150 (MQ=255)
                                                                                                                   aTAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTATAACCTTAAAATGGCTTTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGt                                   >  2:1398046/1‑150 (MQ=255)
                                                                                                                    tAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACACGTGGtt                                  >  2:693095/1‑150 (MQ=255)
                                                                                                                      aGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAATCTCGTGGtttt                                <  1:940723/150‑1 (MQ=255)
                                                                                                                      aGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCACATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtc                               <  2:46181/150‑1 (MQ=255)
                                                                                                                       gATTCGTTTTTAAGCTAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCATTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCAATATTGAGGGAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtc                               <  2:879381/150‑1 (MQ=255)
                                                                                                                        aTTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATTGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTAGTGGTTTtca                              >  1:208892/1‑150 (MQ=255)
                                                                                                                         ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                             <  1:1030899/150‑1 (MQ=255)
                                                                                                                         ttCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGAATTATTGAGGAAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcag                             >  2:1327570/1‑150 (MQ=255)
                                                                                                                           cGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTtcagtc                           <  2:1345165/150‑1 (MQ=255)
                                                                                                                                 tAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCCCATCCTGTTTTTAACCTTAAAAGGGCATTATTGAGGTAGACCTACAGGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                     >  2:470024/1‑150 (MQ=255)
                                                                                                                                 tAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                     >  2:449029/1‑150 (MQ=255)
                                                                                                                                 tAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                     >  2:579066/1‑150 (MQ=255)
                                                                                                                                 atttcAAGCGAGAGCAGGTGATCCTGCTCTCGTTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGt                     >  1:1045780/5‑150 (MQ=39)
                                                                                                                                  aaGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTGCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTGTTCAGTAAGAAGTg                    >  1:341202/1‑150 (MQ=255)
                                                                                                                                   aGCAAACGAGAGCAGTATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGa                   <  2:58664/150‑1 (MQ=255)
                                                                                                                                   aGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCCCCTCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTAGTGGTTTTCAGTCAGAATTGa                   >  2:666674/1‑150 (MQ=255)
                                                                                                                                   aGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCGCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTCCATGAAAGGACAAGAAACTCGTGCTTTTCACTCAGAAGTGa                   >  1:805750/1‑150 (MQ=255)
                                                                                                                                     cAAGCGAGAGCAGGTGATCCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGaaa                 >  2:588350/1‑150 (MQ=255)
                                                                                                                                      aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAc                                                                                         <  2:1595879/77‑1 (MQ=255)
                                                                                                                                      aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAAc                                                                                         >  1:1595879/1‑77 (MQ=255)
                                                                                                                                      aaaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAAc                >  2:247119/1‑150 (MQ=255)
                                                                                                                                       agcGAGAGCAGGTGATCCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACa               >  1:1408197/3‑150 (MQ=255)
                                                                                                                                       aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACa               >  1:611468/1‑150 (MQ=255)
                                                                                                                                       aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACa               >  2:1349767/1‑150 (MQ=255)
                                                                                                                                       aaCGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCACCAACTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGCCCTACATGAAAGGACAAGAAACTCGTTGTTTTCAGTCAGAAGTGAAACa               >  2:881553/1‑150 (MQ=255)
                                                                                                                                            gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTc          <  1:726337/150‑1 (MQ=255)
                                                                                                                                            gagCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTc          <  2:1534609/150‑1 (MQ=255)
                                                                                                                                              gCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACTTCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTg        <  1:1361597/150‑1 (MQ=255)
                                                                                                                                               cAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGc       >  1:1657589/1‑150 (MQ=255)
                                                                                                                                               cAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGc       <  2:155458/150‑1 (MQ=255)
                                                                                                                                                    tCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACTTCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCtg  <  2:254088/150‑1 (MQ=255)
                                                                                                                                                    tCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCtg  >  1:1734103/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGCTTTGCGCGCAATATGACGTGGAAGCCAGCCGAATCGCTCATGGCGTTCCGGTTGGCGGCGAGCTGGAAATGGTCGACGGCACCACGTTGTCACACTCCCTTGCCGGGCGTCATAAGATTCGTTTTTAAGCAAACGAGAGCAGGATCACCTGCTCTCGCTTGAAATTATTCTCCCTTGTCCCCATCTCTCCCACATCCTGTTTTTAACCTTAAAATGGCATTATTGAGGTAGACCTACATGAAAGGACAAGAAACTCGTGGTTTTCAGTCAGAAGTGAAACAGCTTCTGCACCTG  >  NC_000913/494879‑495176

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