Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,755,639 G→A 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.8% 289.5 / ‑4.5 85intergenic (‑498/‑59)ydhZ/pykFuncharacterized protein/pyruvate kinase I
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (40/44);  total (41/44)
Fisher's exact test for biased strand distribution p-value = 4.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

TGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755543‑1755735
                                                                                                |                                                                                                
tGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAg                                                                                               <  1:2358747/100‑1 (MQ=255)
tGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAg                                                                                               <  1:690257/100‑1 (MQ=255)
tGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAg                                                                                               <  2:1613741/100‑1 (MQ=255)
   gaCACAATCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgc                                                                                            >  1:1523796/1‑100 (MQ=255)
     cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          >  1:2627722/1‑100 (MQ=255)
     cacaAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTgcgc                                                                                          >  2:1090099/1‑100 (MQ=255)
       caAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGccc                                                                                        <  2:1693716/100‑1 (MQ=255)
            acacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaaa                                                                                   <  1:2124211/100‑1 (MQ=255)
            acacaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGaaa                                                                                   <  1:1719445/100‑1 (MQ=255)
              acaTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGc                                                                                 <  1:969267/100‑1 (MQ=255)
               caTTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCa                                                                                >  2:672747/1‑100 (MQ=255)
                aTTCCTCTGCACGCTTTTTCTTTGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCaa                                                                               >  2:647384/1‑100 (MQ=255)
                 ttCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAg                                                                              >  1:2474932/1‑100 (MQ=255)
                   ccTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGtt                                                                            >  2:1061412/1‑100 (MQ=255)
                    ctctGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGttt                                                                           >  1:743058/1‑100 (MQ=255)
                       tGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtctc                                                                        <  2:2803162/100‑1 (MQ=255)
                       cttccgatctttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtctc                                                                        <  2:820849/91‑1 (MQ=255)
                        gCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTcc                                                                       >  2:65609/1‑100 (MQ=255)
                          aCGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                     <  2:1961268/100‑1 (MQ=255)
                          aCGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCa                                                                     <  2:2875493/100‑1 (MQ=255)
                           cGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCAt                                                                    >  1:631610/1‑100 (MQ=255)
                           cGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCAt                                                                    >  1:2058133/1‑100 (MQ=255)
                            gCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATc                                                                   >  2:1215098/1‑100 (MQ=255)
                             cTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATcc                                                                  <  1:1043533/100‑1 (MQ=255)
                             cTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATcc                                                                  <  2:2815955/100‑1 (MQ=255)
                                tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTGCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  2:1781326/100‑1 (MQ=255)
                                tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  1:1371632/100‑1 (MQ=255)
                                tttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               <  1:672749/100‑1 (MQ=255)
                                 ttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTtctcagatcggaag                                                              >  1:820850/1‑90 (MQ=255)
                                 ttCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtct                                                              >  2:283971/1‑100 (MQ=255)
                                  tCGATTTCACCTATACTTAGAGCGAGGCACCACAACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtc                                                               >  2:1477307/1‑98 (MQ=255)
                                  tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                             <  1:1470157/100‑1 (MQ=255)
                                  tCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTtctc                                                             <  1:1363454/100‑1 (MQ=255)
                                   cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            <  2:2773712/100‑1 (MQ=255)
                                   cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            <  1:1596225/100‑1 (MQ=255)
                                   cGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCa                                                            <  2:1053095/100‑1 (MQ=255)
                                      tGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACt                                                         <  1:24098/100‑1 (MQ=255)
                                      tGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACt                                                         >  1:2658925/1‑100 (MQ=255)
                                       gTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACtt                                                        >  2:1080449/1‑100 (MQ=255)
                                       gTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACtt                                                        >  2:971557/1‑100 (MQ=255)
                                         cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCTGATAGTGCGCCCAGAAAGCAAGTTTCTCCCATCATTCTCAACTTaa                                                      >  1:1431128/1‑100 (MQ=255)
                                         cACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaa                                                      >  1:441495/1‑100 (MQ=255)
                                          aCCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTaaa                                                     <  1:1304459/100‑1 (MQ=255)
                                          aCCTATCCTGAGAGCCAGGCACCAACCCTTTCGTAATGCCGGATTCGCTGTCCGACTGTGCGCTCAGTGAGCAAGTTTCTCCCAGCCTTCTTAACTTaaa                                                     >  2:236656/1‑100 (MQ=255)
                                           ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    >  1:1355025/1‑100 (MQ=255)
                                           ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  1:1275415/100‑1 (MQ=255)
                                           ccTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaag                                                    <  2:493977/100‑1 (MQ=255)
                                             tATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagac                                                  >  2:2487691/1‑100 (MQ=255)
                                                   tAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagac                                            <  1:1888294/100‑1 (MQ=255)
                                                   tAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagac                                            <  2:1551857/100‑1 (MQ=255)
                                                   tAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAaagactaagac                                            <  2:2466405/100‑1 (MQ=255)
                                                      agCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGt                                         <  2:826626/100‑1 (MQ=255)
                                                        cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                       <  1:519568/100‑1 (MQ=255)
                                                        cGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCa                                       >  2:2059377/1‑100 (MQ=255)
                                                           ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                    >  2:296652/1‑100 (MQ=255)
                                                           ggCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGa                                    <  2:2870879/100‑1 (MQ=255)
                                                            gCACCACCACTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAGGACTGTCATGaa                                   <  1:2371679/100‑1 (MQ=255)
                                                              accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                 >  1:2078841/1‑100 (MQ=255)
                                                              accaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaa                                 >  2:2885764/1‑100 (MQ=255)
                                                               ccaccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGaaaaa                                <  2:441279/100‑1 (MQ=255)
                                                                caccacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAg                               >  2:362369/1‑100 (MQ=255)
                                                                   cacTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAAACAGTCATAAAAAAGAcc                            <  2:1971371/100‑1 (MQ=255)
                                                                    acTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCa                           <  2:1104550/100‑1 (MQ=255)
                                                                     cTTTCGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCaa                          >  2:2014296/1‑100 (MQ=255)
                                                                         cGTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAAt                       <  1:1383320/99‑1 (MQ=255)
                                                                          gTAATACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACAAAAATTg                     >  2:2400584/1‑100 (MQ=255)
                                                                             aTACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGttt                  >  1:954370/1‑100 (MQ=255)
                                                                              tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 <  1:65609/100‑1 (MQ=255)
                                                                              tACCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTg                 >  1:650333/1‑100 (MQ=255)
                                                                               aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                >  2:2491533/1‑100 (MQ=255)
                                                                               aCCGGATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGc                <  2:1470465/100‑1 (MQ=255)
                                                                                  ggATTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCAcc             >  1:1145372/1‑100 (MQ=255)
                                                                                    aTTCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCAt           >  1:1118901/1‑100 (MQ=255)
                                                                                     ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc          <  1:438490/100‑1 (MQ=255)
                                                                                     ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc          <  1:647386/100‑1 (MQ=255)
                                                                                     ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc          >  2:591130/1‑100 (MQ=255)
                                                                                     ttCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATc          >  2:2226097/1‑100 (MQ=255)
                                                                                      tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCg         <  1:1090102/100‑1 (MQ=255)
                                                                                      tCGCTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCg         <  1:1061414/100‑1 (MQ=255)
                                                                                       cGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCgg        >  1:315995/1‑100 (MQ=255)
                                                                                         cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      <  2:2627703/100‑1 (MQ=255)
                                                                                         cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      >  2:1150475/1‑100 (MQ=255)
                                                                                         cTTTCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGAc      >  1:2480330/1‑100 (MQ=255)
                                                                                          tttCCGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCGTCGGAcc     <  2:1885842/100‑1 (MQ=255)
                                                                                             ccGACAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTACACCATCGGACCGaa  >  1:2847153/1‑100 (MQ=255)
                                                                                                |                                                                                                
TGAGACACAAGCACACATTCCTCTGCACGCTTTTTCGATGTCACCTATCCTTAGAGCGAGGCACCACCACTTTCGTAATACCGGATTCGCTTTCCGGCAGTGCGCCCAGAAAGCAAGTTTCTCCCATCCTTCTCAACTTAAAGACTAAGACTGTCATGAAAAAGACCAAAATTGTTTGCACCATCGGACCGAA  >  NC_000913/1755543‑1755735

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