Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,465,236 C→T 100% P315P (CCG→CCA treB ← trehalose‑specific PTS enzyme: IIB and IIC component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,465,2360CT97.8% 295.5 / ‑2.2 95P315P (CCG→CCAtreBtrehalose‑specific PTS enzyme: IIB and IIC component
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (56/35);  total (59/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.24e-01

ACAGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCAT  >  NC_000913/4465139‑4465331
                                                                                                 |                                                                                               
aCAGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTgg                                                                                               <  1:1499244/100‑1 (MQ=255)
 cAGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTggg                                                                                              <  1:1881434/100‑1 (MQ=255)
 cAGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTggg                                                                                              >  2:1138355/1‑100 (MQ=255)
  aGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTgggg                                                                                             <  1:1499438/100‑1 (MQ=255)
  aGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTgggg                                                                                             <  2:372104/100‑1 (MQ=255)
    cgcAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCg                                                                                           >  1:393060/1‑100 (MQ=255)
    cgatcccAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCg                                                                                           <  2:1319936/94‑1 (MQ=255)
      cAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTa                                                                                         >  2:1005269/1‑100 (MQ=255)
      cAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTa                                                                                         >  2:1779239/1‑100 (MQ=255)
      cAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTa                                                                                         >  1:1836364/1‑100 (MQ=255)
      cAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTa                                                                                         >  2:1244418/1‑100 (MQ=255)
        aTCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACa                                                                                       >  1:2066874/1‑100 (MQ=255)
          cAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCgagatcg                                                                                     >  1:1319943/1‑94 (MQ=255)
           aGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGa                                                                                    >  2:106398/1‑100 (MQ=255)
           aGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGa                                                                                    >  2:871053/1‑100 (MQ=255)
           aGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGa                                                                                    >  2:1995947/1‑100 (MQ=255)
              ggCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAg                                                                                  >  2:2013427/1‑99 (MQ=255)
              ggCCACACTGGCGTACCACCCATGCTTTGAATCATCTACAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGaa                                                                                   >  2:1448799/1‑98 (MQ=255)
                 cacaCTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCan                                                                              <  2:2245938/100‑2 (MQ=255)
                   cacTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAAc                                                                            >  1:1046950/1‑100 (MQ=255)
                      tGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAAt                                                                         <  2:2325481/100‑1 (MQ=255)
                      tGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAAt                                                                         <  2:273345/100‑1 (MQ=255)
                       ggCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAAtg                                                                        >  2:2013492/1‑100 (MQ=255)
                       ggCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAAtg                                                                        <  2:841725/100‑1 (MQ=255)
                             ccaccCATGCTTTGATTCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTgcgc                                                                  <  1:2345987/100‑1 (MQ=255)
                              caccCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGcc                                                                 <  2:1095977/100‑1 (MQ=255)
                              caccCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGcc                                                                 >  1:1994000/1‑100 (MQ=255)
                               accCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCa                                                                >  1:1031604/1‑100 (MQ=255)
                                cccATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCaa                                                               >  2:680112/1‑100 (MQ=255)
                                cccATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCaa                                                               <  2:2193926/100‑1 (MQ=255)
                                cccATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCa                                                                >  2:1974204/1‑99 (MQ=255)
                                  cATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATc                                                             <  1:751259/100‑1 (MQ=255)
                                     gCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGa                                                          <  2:1263469/100‑1 (MQ=255)
                                     gCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGa                                                          >  2:342456/1‑100 (MQ=255)
                                      cTTTAAAGCATCAGCAAATCACTGGCAAGCGTGGTCTGGTGTACACCGGTGTTCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAg                                                         <  2:97081/100‑1 (MQ=255)
                                       tttGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGc                                                        >  1:1285432/1‑100 (MQ=255)
                                          gAATCATCTGCAAATCAATAGCAATCGTGGTCTAGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCACACAATGCTGCGCCAATCGGAGCa                                                       >  2:2120626/1‑98 (MQ=255)
                                           aaTCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAg                                                    <  2:553577/100‑1 (MQ=255)
                                            atcatcTGCAAATCAATAGCAGGCGTGGTCTGGTGTACCCCGGTGATCACCAGAGGGGCGTAAAGGAAGCCACACAATGCTGCGCCAATCgggtcaa                                                      >  2:778106/1‑92 (MQ=255)
                                            atcatcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAg                                                    <  2:721623/99‑1 (MQ=255)
                                              catcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTg                                                 >  2:2095262/1‑100 (MQ=255)
                                              catcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTg                                                 >  2:2080076/1‑100 (MQ=255)
                                              catcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTg                                                 >  1:1707115/1‑100 (MQ=255)
                                              catcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTg                                                 >  2:1196642/1‑100 (MQ=255)
                                              catcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTg                                                 <  1:2077707/100‑1 (MQ=255)
                                                tcTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGcc                                               <  1:2117207/100‑1 (MQ=255)
                                                 cTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCg                                              >  1:1968952/1‑100 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCTAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             >  1:447195/1‑100 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             >  1:329813/1‑100 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             <  2:2336209/100‑1 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             <  1:1594471/100‑1 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             <  2:1646591/100‑1 (MQ=255)
                                                  tGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCgg                                             >  1:946327/1‑100 (MQ=255)
                                                   gCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCg                                              <  2:330687/98‑1 (MQ=255)
                                                     aaaTCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGtca                                          >  2:2279558/1‑100 (MQ=255)
                                                     aaaTCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGtca                                          <  2:393058/100‑1 (MQ=255)
                                                     aaaTCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGtca                                          <  2:476037/100‑1 (MQ=255)
                                                       atgaaTAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGGGGGGCGTACAGGGAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCt                                                  >  2:1005172/4‑90 (MQ=255)
                                                       aTCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGtcatc                                        >  1:1306713/1‑100 (MQ=255)
                                                         cAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAg                                      <  2:1262601/100‑1 (MQ=255)
                                                         cAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAg                                      <  1:186507/100‑1 (MQ=255)
                                                          aaTAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAgg                                     >  1:2279568/1‑100 (MQ=255)
                                                            tAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGTGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                                   <  2:1455612/100‑1 (MQ=255)
                                                            tAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                                   >  2:275127/1‑100 (MQ=255)
                                                            tAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                                   <  1:145844/100‑1 (MQ=255)
                                                            tAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAAGGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTg                                   >  2:2371844/1‑100 (MQ=255)
                                                             aGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGa                                  >  2:1458463/1‑100 (MQ=255)
                                                               cgaGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACg                                >  1:1616153/3‑100 (MQ=255)
                                                               cAAGCGTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACg                                <  2:1480289/100‑1 (MQ=255)
                                                                    gTGGTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCg                           <  1:1531652/100‑1 (MQ=255)
                                                                      ggTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCa                         <  1:1465333/100‑1 (MQ=255)
                                                                       gTCTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCaa                        >  1:2175367/1‑100 (MQ=255)
                                                                         cTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAg                      <  1:414727/100‑1 (MQ=255)
                                                                         cTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAg                      <  1:2080093/100‑1 (MQ=255)
                                                                         cTGGTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAg                      <  1:560389/100‑1 (MQ=255)
                                                                           ggTGTACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGc                    >  2:1893581/1‑100 (MQ=255)
                                                                             gccccactggTGATCACCGGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCCGAACAAAGCCGCCGGTCATCAGGTGACGTACCGCAAAGGCaa                  >  2:1672183/9‑100 (MQ=255)
                                                                               tACACCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCTTCAGGTGACGTACCGCAAAGGCAACg                <  1:1455292/100‑1 (MQ=255)
                                                                                 cacCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGCCGTACCGCAACGGCAACGcc              >  2:351904/1‑100 (MQ=255)
                                                                                 cacCGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGcc              >  2:1227675/1‑100 (MQ=255)
                                                                                   ccGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCa             >  1:1271382/1‑99 (MQ=255)
                                                                                   ccGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCAt            >  2:1069945/1‑100 (MQ=255)
                                                                                   ccGGTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCAt            >  2:1290135/1‑100 (MQ=255)
                                                                                     ggTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCg          >  2:233141/1‑100 (MQ=255)
                                                                                     ggTGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCg          >  1:67422/1‑100 (MQ=255)
                                                                                       tGATCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGcc        >  1:1072481/1‑100 (MQ=255)
                                                                                          tCACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGCCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAAt     >  2:2217991/1‑100 (MQ=255)
                                                                                           cACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAAt     >  2:1092672/1‑99 (MQ=255)
                                                                                           cACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATc    >  1:1317141/1‑100 (MQ=255)
                                                                                           cACCAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATc    >  1:2113207/1‑100 (MQ=255)
                                                                                            aCCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCa   >  1:2300251/1‑100 (MQ=255)
                                                                                             ccAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCAt  >  1:416799/1‑100 (MQ=255)
                                                                                             ccAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCAt  >  2:610735/1‑100 (MQ=255)
                                                                                             ccAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACCTACCGCAAAGGCAACTCCATCGCCAATCAt  >  1:2011850/1‑100 (MQ=255)
                                                                                             ccAGTGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGACATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCAt  >  1:709661/1‑100 (MQ=255)
                                                                                                 |                                                                                               
ACAGCGCAATCAGCGGCCACACTGGCGTACCACCCATGCTTTGAATCATCTGCAAATCAATAGCAAGCGTGGTCTGGTGTACACCGGTGATCACCAGCGGGGCGTACAGGAAGCCAAACAATGCTGCGCCAATCGGAGCAAAGCTGCCGGTCATCAGGTGACGTACCGCAAAGGCAACGCCATCGCCAATCAT  >  NC_000913/4465139‑4465331

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