Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,755,639 G→A 100% intergenic (‑498/‑59) ydhZ ← / → pykF uncharacterized protein/pyruvate kinase I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,755,6390GA98.9% 315.2 / ‑4.7 93intergenic (‑498/‑59)ydhZ/pykFuncharacterized protein/pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (51/40);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755542‑1755735
                                                                                                 |                                                                                                
cTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCa                                                                                                <  1:1714648/100‑1 (MQ=255)
  gagaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGt                                                                                              >  1:1926090/1‑100 (MQ=255)
   agaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTg                                                                                             <  1:469557/100‑1 (MQ=255)
      cacaAGCACACATTCCTCTGCACGCTTTTTTGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          >  2:2384711/1‑100 (MQ=255)
       acaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGcc                                                                                         <  2:1903835/100‑1 (MQ=255)
        caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        <  1:691392/100‑1 (MQ=255)
        caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        >  1:1544106/1‑100 (MQ=255)
         aaGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        <  1:1390367/99‑1 (MQ=255)
          aGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAg                                                                                      <  2:2325810/100‑1 (MQ=255)
          aGCACACATTCCTCTGCACGCTTTTGCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCTTAATACCGGATTCGCTTTCCGACAGTGCGCCCAg                                                                                      <  2:1150939/100‑1 (MQ=255)
            cacacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaa                                                                                    <  2:2340217/100‑1 (MQ=255)
            cacacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaa                                                                                    <  2:2339695/100‑1 (MQ=255)
             acacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaaa                                                                                   >  1:1224916/1‑100 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCTCCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  <  2:1006572/100‑1 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  >  2:1696056/1‑100 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  <  1:836728/100‑1 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  >  1:404512/1‑100 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  >  2:2225062/1‑100 (MQ=255)
              cacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAg                                                                                  <  1:1411841/100‑1 (MQ=255)
                   tCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGt                                                                             <  1:1736985/100‑1 (MQ=255)
                    ccTCTGCACGCTTTTTTGATTTTACTTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            >  2:1884799/1‑100 (MQ=255)
                     ctctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGttt                                                                           >  1:739074/1‑100 (MQ=255)
                     ctctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGttt                                                                           >  2:1503948/1‑100 (MQ=255)
                      tctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtc                                                                          >  1:1189776/1‑100 (MQ=255)
                       ctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACCGTGCGCCCAGAAAACAAGTTtct                                                                         >  2:786080/1‑100 (MQ=255)
                         gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGCCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       >  1:1677279/1‑100 (MQ=255)
                          cACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       <  2:776350/99‑1 (MQ=255)
                           aCGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                     <  1:719310/100‑1 (MQ=255)
                            cGCTTTTTCGATATCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                     <  2:1125340/99‑1 (MQ=255)
                               tttttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCt                                                                 >  1:482669/1‑100 (MQ=255)
                                  ttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtct                                                              <  2:499118/100‑1 (MQ=255)
                                  ttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtct                                                              >  1:1219752/1‑100 (MQ=255)
                                   tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                             <  1:739809/100‑1 (MQ=255)
                                   tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTCCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                             >  2:2335172/1‑100 (MQ=255)
                                    cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            <  2:1022574/100‑1 (MQ=255)
                                    cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            >  1:1120463/1‑100 (MQ=255)
                                     gATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCaa                                                           >  2:2086355/1‑100 (MQ=255)
                                      aTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                          <  2:2342588/100‑1 (MQ=255)
                                         tCACCGATCCTTAGAGCGAGGCACCACCACTTTCGTATTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTa                                                       >  2:1420320/1‑100 (MQ=255)
                                          cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      <  1:901516/100‑1 (MQ=255)
                                           aCCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaaa                                                     <  1:1159412/100‑1 (MQ=255)
                                           aCCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaaa                                                     <  1:1820723/100‑1 (MQ=255)
                                           aCCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaaa                                                     >  1:1541851/1‑100 (MQ=255)
                                            ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    >  2:397201/1‑100 (MQ=255)
                                            ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    >  2:303227/1‑100 (MQ=255)
                                            ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  2:1592845/100‑1 (MQ=255)
                                                  cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  1:590330/1‑100 (MQ=255)
                                                  cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                              >  1:1418462/1‑100 (MQ=255)
                                                   ttAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaaga                                             <  2:676323/100‑1 (MQ=255)
                                                   ttAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaaga                                             >  2:1660134/1‑100 (MQ=255)
                                                    tAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagac                                            >  1:3665/1‑100 (MQ=255)
                                                      gagCGAGGCACCACCTCTTTCGTTATACCTGATTCGCTTTCCGTCAGTGCGCTCAGACAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagactn                                          <  2:1665229/100‑2 (MQ=255)
                                                        gCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTc                                        >  1:40822/1‑100 (MQ=255)
                                                         cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                       >  1:1664598/1‑100 (MQ=255)
                                                          gAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCAt                                      >  1:1751880/1‑100 (MQ=255)
                                                          gAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCAt                                      >  2:2103329/1‑100 (MQ=255)
                                                           aGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCACGTTTCCCCCATCCTTCTCAACTTAAAGACTAAGACTGGCATg                                     >  2:711219/1‑100 (MQ=255)
                                                             gCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaa                                   <  2:1194896/100‑1 (MQ=255)
                                                              caccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaa                                  <  2:1393403/100‑1 (MQ=255)
                                                               accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                 >  2:2387821/1‑100 (MQ=255)
                                                                 caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               <  1:2362330/100‑1 (MQ=255)
                                                                 caccacTTTCGTAATACCGGATTCGCTTTCCGACAGAGCGCCCCGAAAGCAAGATCCTCCCATCCTTCTCAACTTAAAGACTCAGACTGTCATGAAAAAg                               >  2:125010/1‑100 (MQ=255)
                                                                 caccacTTTCGTAATACCGGAGTCGCTTTCCGACAGTGAGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               <  2:1161459/100‑1 (MQ=255)
                                                                  accacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGa                              <  1:2117623/100‑1 (MQ=255)
                                                                   cctcTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCACCTTAAAGACTAAGACTGTCATGAAAAAGAc                             <  1:1112907/97‑1 (MQ=255)
                                                                   ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             >  2:2294068/1‑100 (MQ=255)
                                                                   ccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             >  2:206936/1‑100 (MQ=255)
                                                                   ccacTTTCGTAATACCGGATTCGCGTTCCGACAGGGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGAc                             >  2:1747130/1‑100 (MQ=255)
                                                                    cacaTTCATAATACCAGATTCGCTCTGCGACAGTGCGTTCAGAAAGCAAGTTTCTCCAATCCATCTCGACTTAAAGACTAAGTCCGTCATGAAAAAg                               <  2:109239/97‑1 (MQ=255)
                                                                       tttCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCaaa                         >  2:1927914/1‑100 (MQ=255)
                                                                          cGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAAtt                      >  1:203885/1‑100 (MQ=255)
                                                                          cGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAAtt                      >  1:1674745/1‑100 (MQ=255)
                                                                           gTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTg                     >  1:824949/1‑100 (MQ=255)
                                                                             aaTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGtt                   <  1:1880708/100‑1 (MQ=255)
                                                                              aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  <  1:1696069/100‑1 (MQ=255)
                                                                              aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  >  1:289196/1‑100 (MQ=255)
                                                                              aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  <  1:1538157/100‑1 (MQ=255)
                                                                               tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 >  1:1283592/1‑100 (MQ=255)
                                                                               tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 <  2:1762227/100‑1 (MQ=255)
                                                                               tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 >  1:1932596/1‑100 (MQ=255)
                                                                                 ccGGATTCGCTTTCCGACAGTGCGCCCATAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                <  2:1659215/99‑1 (MQ=255)
                                                                                 ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCCTGAAAAAGACCAAAATTGTTTGCa               >  2:2285036/1‑100 (MQ=255)
                                                                                  cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc              >  2:2327529/1‑100 (MQ=255)
                                                                                  cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc              <  2:1206073/100‑1 (MQ=255)
                                                                                     aTTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCAt           >  1:2341464/1‑100 (MQ=255)
                                                                                       tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCg         <  1:1884814/100‑1 (MQ=255)
                                                                                        cGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCgg        >  2:2118924/1‑100 (MQ=255)
                                                                                        cGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCgg        >  1:1370636/1‑100 (MQ=255)
                                                                                          cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      >  2:2328470/1‑100 (MQ=255)
                                                                                           tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc     <  1:240259/100‑1 (MQ=255)
                                                                                           tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc     >  2:896135/1‑100 (MQ=255)
                                                                                            ttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcca    >  1:2178122/1‑99 (MQ=255)
                                                                                              ccGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGaa  <  1:1827215/100‑1 (MQ=255)
                                                                                                 |                                                                                                
CTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755542‑1755735

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: