Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,220,357 G→C 6.3% L545V (CTG→GTG)  bglX ← beta‑D‑glucoside glucohydrolase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,220,3570GC 6.3% 395.7 / 13.5 111L545V (CTG→GTG) bglXbeta‑D‑glucoside glucohydrolase, periplasmic
Reads supporting (aligned to +/- strand):  ref base G (53/51);  new base C (4/3);  total (57/54)
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.70e-01

Only 100 of 121 total aligned reads displayed.

GACGGGTTGTAATCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTTTCGCCGTC  >  NC_000913/2220215‑2220499
                                                                                                                                              |                                                                                                                                              
gACGGGTTGTAATCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGta                                                                                                                                           >  4:93292/1‑147 (MQ=255)
 aCGGGTTGTAATCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGttt                                                                                                                                          >  1:215562/1‑148 (MQ=255)
    ggTTGTAATCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTAcc                                                                                                                                       <  3:926852/148‑1 (MQ=255)
            tCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGTAAACCGCTgg                                                                                                                                  <  2:1040440/145‑1 (MQ=255)
            tCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGTAAACCGCTgg                                                                                                                                  >  1:1040440/1‑145 (MQ=255)
             cGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCCCCAGCGCCAGCGGACGCCCGTTCATCAGCaccagcaccagc                                                                                                                                               >  4:979634/1‑131 (MQ=255)
                 aaaCAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCgg                                                                                                                                             <  2:66022/129‑1 (MQ=255)
                 aaaCAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCgg                                                                                                                                             >  1:66022/1‑129 (MQ=255)
                  aaCAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCgg                                                                                                                                             >  1:1112976/1‑128 (MQ=255)
                    cAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCgg                                                                                                                                     >  1:88185/1‑134 (MQ=255)
                            cGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCaa                                                                                                               >  4:1143857/1‑148 (MQ=255)
                               cAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCg                                                                                                                                      <  4:797881/122‑1 (MQ=255)
                               cAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCg                                                                                                                                      >  3:797881/1‑122 (MQ=255)
                                           cGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTg                                                                                                <  4:493028/148‑1 (MQ=255)
                                             ccTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGt                                                                                                                                            >  4:174260/1‑102 (MQ=255)
                                                 cAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTc                                                                                                            >  1:276716/1‑130 (MQ=255)
                                                 cAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTc                                                                                                            <  2:276716/130‑1 (MQ=255)
                                                    tCCCCGCAAAGCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGt                                                                                       >  3:747937/1‑148 (MQ=255)
                                                     ccccGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTg                                                                                      <  3:1143857/148‑1 (MQ=255)
                                                       ccGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGc                                                                                                                                 >  1:584895/1‑103 (MQ=255)
                                                        cGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGTAAACCGCTg                                                                                                                                   <  2:845696/100‑1 (MQ=255)
                                                        cGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGTAAACCGCTg                                                                                                                                   >  1:845696/1‑100 (MQ=255)
                                                         gCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGtaaa                                                                                                                                         >  1:1021169/1‑90 (MQ=255)
                                                         gCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCACCGGtaaa                                                                                                                                         <  2:1021169/93‑4 (MQ=255)
                                                              ccAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGcgcc                                                                                                                       >  4:327692/1‑106 (MQ=255)
                                                              ccAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGcgcc                                                                                                                       <  3:327692/106‑1 (MQ=255)
                                                                aGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGtt                                                                                                                                           <  1:238337/84‑1 (MQ=255)
                                                                aGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGtt                                                                                                                                           >  2:238337/1‑84 (MQ=255)
                                                                aGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCg                                                                           >  2:372049/1‑148 (MQ=255)
                                                                  gTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCttt                                                                                                 <  4:844081/124‑1 (MQ=255)
                                                                      ccAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGa                                                                     <  1:1007783/148‑1 (MQ=255)
                                                                        aGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAgg                                                                   >  3:475030/1‑148 (MQ=255)
                                                                         gAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGc                                                                  >  1:1066602/1‑148 (MQ=255)
                                                                           aTCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCtt                                                                <  2:1066602/148‑1 (MQ=255)
                                                                           aTCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCtt                                                                >  4:910061/1‑148 (MQ=255)
                                                                           aTCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCtt                                                                <  1:372049/148‑1 (MQ=255)
                                                                               cATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTg                                                            <  4:475030/148‑1 (MQ=255)
                                                                                aTCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCaccagcaccagc                                                                                                                                               <  3:1115614/64‑1 (MQ=255)
                                                                                aTCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCaccagcaccagc                                                                                                                                               >  4:1115614/1‑64 (MQ=255)
                                                                                 tCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCtt                                                                                                                               >  1:447430/1‑79 (MQ=255)
                                                                                   aGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACg                                                                                                         <  4:594516/99‑1 (MQ=255)
                                                                                        gctgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGc                                                                                                    >  2:41515/1‑99 (MQ=255)
                                                                                        gctgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGc                                                                                                    <  1:41515/99‑1 (MQ=255)
                                                                                        gctgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGa                                                                                     >  1:84718/1‑114 (MQ=255)
                                                                                        gctgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGa                                                                                     <  2:84718/114‑1 (MQ=255)
                                                                                        gctgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATccc                                                   <  1:1001251/148‑1 (MQ=255)
                                                                                          tgATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCt                                                 <  2:215562/148‑1 (MQ=255)
                                                                                            aTCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATcc                                                    <  2:373867/143‑1 (MQ=255)
                                                                                            aTCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCtgt                                               >  3:270589/1‑148 (MQ=255)
                                                                                             tcttctTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCgg                                                                                                                                             >  1:51458/1‑53 (MQ=255)
                                                                                             tcttctTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCgg                                                                                                                                             <  2:51458/53‑1 (MQ=255)
                                                                                             tcttctTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTg                                                            >  4:417837/1‑134 (MQ=255)
                                                                                             tcttctTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTg                                                            <  3:417837/134‑1 (MQ=255)
                                                                                              cttcttTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTg                                                            <  2:734628/133‑1 (MQ=255)
                                                                                                tcttTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCtgtg                                              <  2:1090823/145‑1 (MQ=255)
                                                                                                tcttTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCtt                                           >  3:132519/1‑148 (MQ=255)
                                                                                                tcttTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCtt                                           >  3:271971/1‑148 (MQ=255)
                                                                                                 cttTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTc                                          <  3:910061/148‑1 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCaa                                                                                                                   >  2:860261/1‑73 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCaa                                                                                                                   <  1:860261/73‑1 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTgg                                                                      >  2:353443/1‑118 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTgg                                                                      <  1:353443/118‑1 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCt                                                                 <  4:510799/123‑1 (MQ=255)
                                                                                                   ttCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCt                                                                 >  3:510799/1‑123 (MQ=255)
                                                                                                       ccagcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCttt                                                                                                                              <  2:249765/58‑1 (MQ=255)
                                                                                                       ccagcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGAc                                    >  1:56146/1‑148 (MQ=255)
                                                                                                        cagcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCa                                                                                                                    <  2:1106268/67‑1 (MQ=255)
                                                                                                        cagcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCa                                                                                                                    >  1:1106268/1‑67 (MQ=255)
                                                                                                        cagcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACt                                   <  4:132519/148‑1 (MQ=255)
                                                                                                         agcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCttt                                                                                                                              >  2:480597/1‑56 (MQ=255)
                                                                                                         agcgccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCttt                                                                                                                              <  1:480597/56‑1 (MQ=255)
                                                                                                            gccagcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAg                               <  2:56146/148‑1 (MQ=255)
                                                                                                               agcgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGcca                            <  1:291394/148‑1 (MQ=255)
                                                                                                                 cgGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCa                                                                                                                >  3:413812/1‑62 (MQ=255)
                                                                                                                     cGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGa                                     <  2:217071/133‑1 (MQ=255)
                                                                                                                     cGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGa                                     >  1:217071/1‑133 (MQ=255)
                                                                                                                      gCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCCGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCcaccaccac                     <  4:747937/148‑1 (MQ=255)
                                                                                                                       cccGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGata                                                                                   <  3:987900/85‑1 (MQ=255)
                                                                                                                       cccGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCAt                                                      >  4:413546/1‑114 (MQ=255)
                                                                                                                       cccGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCAt                                                      <  3:413546/114‑1 (MQ=255)
                                                                                                                         cGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATa                                                                                         >  4:345772/1‑77 (MQ=255)
                                                                                                                         cGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATc                  >  1:316829/1‑148 (MQ=255)
                                                                                                                         cGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATc                  >  1:705540/1‑148 (MQ=255)
                                                                                                                         cGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATc                  >  4:887445/1‑148 (MQ=255)
                                                                                                                          gTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCa                 <  4:270589/148‑1 (MQ=255)
                                                                                                                           ttCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCg                                                                           <  3:157557/89‑1 (MQ=255)
                                                                                                                               tcaGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGa                                                                                            >  4:861485/1‑68 (MQ=255)
                                                                                                                               tcaGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGa                                                                                            <  3:861485/68‑1 (MQ=255)
                                                                                                                                cagcaccagcacCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACAt                   >  1:530408/1‑140 (MQ=255)
                                                                                                                                cagcaccagcacCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACAt                   <  2:530408/140‑1 (MQ=255)
                                                                                                                                cagcaccagcacCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGt           <  2:559181/148‑1 (MQ=255)
                                                                                                                                  gcaccagcaccaGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGatat                                                                                  <  3:881929/75‑1 (MQ=255)
                                                                                                                                  gcaccagcaccaGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGttt         <  3:887445/148‑1 (MQ=255)
                                                                                                                                  gcaccagcaccaGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTCCTTCACCGACTACAGCCACCACCACATCAGATTGttt         >  3:740343/1‑148 (MQ=255)
                                                                                                                                   caccagcaccagCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTTTc        >  1:782135/1‑148 (MQ=255)
                                                                                                                                    accagcaccagcGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTTTCg       <  3:645390/148‑1 (MQ=255)
                                                                                                                                         caccagcGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTTTCGCCGTc  >  3:185660/1‑148 (MQ=255)
                                                                                                                                         caccagcGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTGTCGCCGTc  >  2:970010/1‑148 (MQ=255)
                                                                                                                                          agcagcGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGa               >  4:761141/3‑134 (MQ=255)
                                                                                                                                          accagcGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGa               <  3:761141/134‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                              
GACGGGTTGTAATCGCCAAACAATACATCGGCAATTGCATTACCGCCTTCAGTCCCCGCAAACCAGGTTTCCAGAATCGCATCAGCCTGCTGATCTTCTTTCACCAGCGCCAGCGGACGCCCGTTCATCAGCACCAGCACCAGCGGTTTACCGGTGGCTTTCAGCGCCGCAATCAAGTCACGTTGGCTTTGCGGAATAGTGATATCGGTCCGGCTGGAGGCTTCGTGCGCCATCCCCTGTGCTTCACCGACTACAGCCACCACCACATCAGATTGTTTCGCCGTC  >  NC_000913/2220215‑2220499

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: