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,6390GA92.3% 251.5 / 9.8 91intergenic (‑498/‑59)ydhZ/pykFuncharacterized protein/pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base A (29/55);  total (30/61)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

GCTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGA  >  NC_000913/1755541‑1755734
                                                                                                  |                                                                                               
gCTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGc                                                                                                <  1:695843/100‑1 (MQ=255)
  tGAGACACAAGCACACATTCCTCTGCCCGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAg                                                                                              <  1:813080/100‑1 (MQ=255)
  tGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAg                                                                                              <  1:2111788/100‑1 (MQ=255)
     gaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTg                                                                                            <  2:2846471/99‑1 (MQ=255)
     gaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgc                                                                                           <  2:2637645/100‑1 (MQ=255)
     gaCACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgc                                                                                           >  2:91553/1‑100 (MQ=255)
       cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                         >  1:81835/1‑100 (MQ=255)
       cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                         <  2:2754535/100‑1 (MQ=255)
        acaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGcc                                                                                        >  2:1396967/1‑100 (MQ=255)
        acaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGcc                                                                                        <  2:592412/100‑1 (MQ=255)
         caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTTTCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                       <  2:924053/100‑1 (MQ=255)
         caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGccc                                                                                       <  2:799673/100‑1 (MQ=255)
         caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                       >  2:62914/1‑100 (MQ=255)
         caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                       <  2:2575663/100‑1 (MQ=255)
           aGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACCTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAg                                                                                     <  1:1789096/100‑1 (MQ=255)
            gCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGa                                                                                    <  1:2413691/100‑1 (MQ=255)
               tcttccgatctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAg                                                                                 <  1:325306/92‑1 (MQ=255)
                     ccTCTGCACGCTACTTCTATGTCACCAATCCTTACAGCAATGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                           <  1:229851/100‑1 (MQ=255)
                       tctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtc                                                                         <  2:1658025/100‑1 (MQ=255)
                       tctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtc                                                                         <  1:1090390/100‑1 (MQ=255)
                          gCACGCTTTTTCGATGTCACCTCTCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                      <  2:1237548/100‑1 (MQ=255)
                           cACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTccc                                                                     <  1:281252/100‑1 (MQ=255)
                           cACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTccc                                                                     >  2:353959/1‑100 (MQ=255)
                            aCGCTTTTTCTATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                    <  1:1910947/100‑1 (MQ=255)
                            aCGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                    <  2:367146/100‑1 (MQ=255)
                               cTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATcc                                                                 >  2:729794/1‑100 (MQ=255)
                                 ttttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCtt                                                               >  1:1031653/1‑100 (MQ=255)
                                  tatCGATGGCACCTATTCTTATAGGGAGGAACCACCACTTTCGTAATACCGAATTCGCTTTACGAGTCGGCGTCCAGAAAGCAAGTGTCTCTCATccat                                                               >  1:218127/3‑97 (MQ=255)
                                   ttCGATGTCACCTGTCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtct                                                             >  2:736371/1‑100 (MQ=255)
                                    tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                            <  2:81836/100‑1 (MQ=255)
                                       aTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                         <  2:2440981/100‑1 (MQ=255)
                                       aTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAAc                                                         <  1:1836356/100‑1 (MQ=255)
                                         gTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACtt                                                       >  1:1829762/1‑100 (MQ=255)
                                           cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                     >  1:130545/1‑100 (MQ=255)
                                           cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                     >  2:514074/1‑100 (MQ=255)
                                             ccTATCCTTAGACCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                   >  2:1911103/1‑100 (MQ=255)
                                               tATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCTCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagac                                                 <  1:1589048/100‑1 (MQ=255)
                                               tATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagac                                                 <  2:2584894/100‑1 (MQ=255)
                                                 tCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                               <  2:2098088/100‑1 (MQ=255)
                                                 tCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                               >  1:1042965/1‑100 (MQ=255)
                                                 tCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                               <  1:196609/100‑1 (MQ=255)
                                                 tCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagacta                                               >  1:814746/1‑100 (MQ=255)
                                                   cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                             >  1:2325825/1‑100 (MQ=255)
                                                   cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                             >  2:264594/1‑100 (MQ=255)
                                                   cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                             >  2:2777115/1‑100 (MQ=255)
                                                   cTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaag                                             >  2:1457012/1‑100 (MQ=255)
                                                    ttAGATCGAGGCACCACTACTTTCGTAACATCGGATTCGCTTTCCGACACTGCGCGCAGAAAGCACGTTTCTCCCATCCTTCTCAACTCAatgactaag                                             <  1:1750611/99‑1 (MQ=255)
                                                    ttAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaaga                                            >  2:1719414/1‑100 (MQ=255)
                                                      agagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagact                                          <  1:264592/100‑1 (MQ=255)
                                                       gagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTg                                         <  2:2832190/100‑1 (MQ=255)
                                                       gagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTg                                         >  2:2587217/1‑100 (MQ=255)
                                                       gagCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTg                                         <  2:1413832/100‑1 (MQ=255)
                                                        agCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGt                                        <  1:1757797/100‑1 (MQ=255)
                                                         gCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTc                                       <  1:1457009/100‑1 (MQ=255)
                                                          cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                      >  1:342973/1‑100 (MQ=255)
                                                          cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                      <  2:1168763/100‑1 (MQ=255)
                                                           gAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTCTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCAt                                     <  1:1786213/100‑1 (MQ=255)
                                                           gAGGCACCACCACTTGCGTAATACCGGATTCGCCTTCCGACAGTGCCCCCAGAAATCAAGTTTCTCCCATCCTTCCCAACTTAAAGACTAAGACTGTCAt                                     <  1:2381638/100‑1 (MQ=255)
                                                            aGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATg                                    <  1:2165353/100‑1 (MQ=255)
                                                             ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATg                                    <  1:2417432/99‑1 (MQ=255)
                                                             ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                   <  1:2128563/100‑1 (MQ=255)
                                                             ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                   <  2:1114170/100‑1 (MQ=255)
                                                             ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                   <  1:211651/100‑1 (MQ=255)
                                                               caccaccacTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaa                                 >  2:1299666/1‑100 (MQ=255)
                                                                accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                <  2:151082/100‑1 (MQ=255)
                                                                accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACCGTCATGaaaa                                <  2:282773/100‑1 (MQ=255)
                                                                 ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                               <  2:304293/100‑1 (MQ=255)
                                                                 ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                               <  2:1748211/100‑1 (MQ=255)
                                                                 ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                               <  2:1747975/100‑1 (MQ=255)
                                                                 ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                               >  1:2215795/1‑100 (MQ=255)
                                                                  naccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                               >  2:1239407/2‑99 (MQ=255)
                                                                   accacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGa                             <  2:1784447/100‑1 (MQ=255)
                                                                    ccacGTTCGTAATACCTCATTCGCCTCCTGACAGTGCTCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGCCTCAGACTGTCATGAAAAAGa                             <  1:2445002/99‑1 (MQ=255)
                                                                          tCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAAt                      <  2:2834638/100‑1 (MQ=255)
                                                                               aTACCGGATTCGCTTTCCTACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                 <  2:2033323/100‑1 (MQ=255)
                                                                               aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGtct                 <  2:1957613/100‑3 (MQ=255)
                                                                                tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                >  2:307218/1‑100 (MQ=255)
                                                                                tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                <  1:354517/100‑1 (MQ=255)
                                                                                 aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc               <  1:1244110/100‑1 (MQ=255)
                                                                                 aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc               <  1:1217060/100‑1 (MQ=255)
                                                                                  ccGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCa              <  2:1875089/100‑1 (MQ=255)
                                                                                   cGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAc             <  2:2251745/100‑1 (MQ=255)
                                                                                      aTTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCAt          >  1:2268465/1‑100 (MQ=255)
                                                                                       ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc         <  2:2169813/100‑1 (MQ=255)
                                                                                           cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAATGACCAAAATTGTTTGCACCATCGCAc     >  1:2799137/1‑100 (MQ=255)
                                                                                           cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc     >  1:885102/1‑100 (MQ=255)
                                                                                            tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc    <  2:1359637/100‑1 (MQ=255)
                                                                                            tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAcc    <  1:1396965/100‑1 (MQ=255)
                                                                                             ttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCg   <  2:2783313/100‑1 (MQ=255)
                                                                                              tCCGACAGTGTGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGa  <  2:1041148/100‑1 (MQ=255)
                                                                                              tCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGa  >  2:2529937/1‑100 (MQ=255)
                                                                                                  |                                                                                               
GCTGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGA  >  NC_000913/1755541‑1755734

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