Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,796,091 A→G intergenic (+26/+21) ygaP → / ← stpA hypothetical protein/StpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,796,0910AG97.7% 396.1 / ‑0.1 129intergenic (+26/+21)ygaP/stpAhypothetical protein/StpA
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (55/71);  total (57/72)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.87e-01

Only 100 of 132 total aligned reads displayed.

AATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCAGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTCGGTGTACGGC  >  NC_000913/2796002‑2796179
                                                                                         |                                                                                        
aaTCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAgcggcg                                                                                       <  1:3291293/93‑1 (MQ=255)
aaTCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAgcggcg                                                                                       >  2:1284322/1‑93 (MQ=255)
  tCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGAc                                                                                     >  1:2286614/1‑93 (MQ=255)
  tCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGAc                                                                                     <  2:3122076/93‑1 (MQ=255)
   cAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACa                                                                                    >  2:410323/1‑93 (MQ=255)
     gCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATc                                                                                  <  1:1249205/93‑1 (MQ=255)
     gCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATc                                                                                  <  1:1341027/93‑1 (MQ=255)
      cGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATcc                                                                                 <  1:1259019/93‑1 (MQ=255)
      cGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATcc                                                                                 <  1:2639731/93‑1 (MQ=255)
         tttttGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGc                                                                              <  1:860112/93‑1 (MQ=255)
          ttttGTGGAAGGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGcc                                                                             <  1:2166478/93‑1 (MQ=255)
           tttGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCt                                                                            >  1:2562801/1‑93 (MQ=255)
           ttagtgGAGTGGCAAGGTTGTTAGATAAGATGCAGGAGAACGAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCt                                                                            <  2:2443269/90‑1 (MQ=255)
            tggtgGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTc                                                                           <  2:443903/91‑1 (MQ=255)
             tgtgGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCa                                                                          >  2:693965/1‑93 (MQ=255)
              gtgGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAg                                                                         <  1:4087012/93‑1 (MQ=255)
               tgGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGt                                                                        <  2:1569411/93‑1 (MQ=255)
                    tGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTa                                                                   <  2:3334104/93‑1 (MQ=255)
                      gCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGa                                                                 <  2:2695342/93‑1 (MQ=255)
                       gaaGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGAt                                                                <  1:1688098/92‑1 (MQ=255)
                       cAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGAt                                                                >  2:2489015/1‑93 (MQ=255)
                        aaGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATc                                                               >  2:3665287/1‑93 (MQ=255)
                        aaGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATc                                                               >  2:3665288/1‑93 (MQ=255)
                        aaGGTCGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATc                                                               <  2:1660073/93‑1 (MQ=255)
                         aGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATcg                                                              >  1:4136115/1‑92 (MQ=255)
                         aGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCa                                                              <  1:3482219/93‑1 (MQ=255)
                         aGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCa                                                              <  2:1690747/93‑1 (MQ=255)
                         aGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCa                                                              <  2:2003208/93‑1 (MQ=255)
                         aGGGTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCa                                                              <  1:3087697/93‑1 (MQ=255)
                          ggTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAg                                                             >  1:2821619/1‑93 (MQ=255)
                          ggTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAg                                                             >  1:73148/1‑93 (MQ=255)
                          ggTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAg                                                             <  2:2045769/93‑1 (MQ=255)
                           gttgttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAgg                                                            <  1:2787604/93‑1 (MQ=255)
                           gttgttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAgg                                                            <  1:2672065/93‑1 (MQ=255)
                            ttgttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGa                                                           <  2:4658280/93‑1 (MQ=255)
                             tgttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGaa                                                          >  1:269096/1‑93 (MQ=255)
                              gttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTATTTAGATCAGGaaa                                                         >  2:269820/1‑93 (MQ=255)
                              gttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGaaa                                                         >  1:1809580/1‑93 (MQ=255)
                              gttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGaaa                                                         <  2:1205612/93‑1 (MQ=255)
                              gttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGaaa                                                         <  2:781556/93‑1 (MQ=255)
                               ttAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAt                                                        <  1:362729/93‑1 (MQ=255)
                                tAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtc                                                       <  1:3584757/93‑1 (MQ=255)
                                  gATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgt                                                     <  1:2489121/93‑1 (MQ=255)
                                  gATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgt                                                     >  2:306215/1‑93 (MQ=255)
                                  gATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgt                                                     <  1:3554912/93‑1 (MQ=255)
                                  gATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgt                                                     >  1:2803963/1‑93 (MQ=255)
                                  gATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgc                                                     >  1:2803964/1‑92 (MQ=255)
                                   aTAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgtc                                                    <  2:3913961/93‑1 (MQ=255)
                                   aTAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgtc                                                    <  1:1820968/93‑1 (MQ=255)
                                   aTAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAAtcgtc                                                    <  1:3851266/93‑1 (MQ=255)
                                     aaGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCga                                                  >  1:4477470/1‑93 (MQ=255)
                                     aaGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCga                                                  <  1:2823552/93‑1 (MQ=255)
                                     aaGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCga                                                  <  2:2754348/93‑1 (MQ=255)
                                      aGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgag                                                 >  2:194560/1‑93 (MQ=255)
                                      aGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgag                                                 >  1:2835125/1‑93 (MQ=255)
                                       gATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgaga                                                <  1:1325668/93‑1 (MQ=255)
                                        aTGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgagag                                               <  2:4029982/93‑1 (MQ=255)
                                        aTGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgagag                                               >  1:4388362/1‑93 (MQ=255)
                                         tcccGTGGATCCAACGAGCGTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCgataga                                              <  1:4458125/91‑1 (MQ=255)
                                          gCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGAt                                             <  2:1553366/93‑1 (MQ=255)
                                           ccGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGAtt                                            <  1:3635839/93‑1 (MQ=255)
                                            cGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGAttt                                           >  1:657179/1‑93 (MQ=255)
                                            cGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGAttt                                           >  1:1932379/1‑93 (MQ=255)
                                             gTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTa                                          >  2:1603160/1‑93 (MQ=255)
                                               ggAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCGCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTAcc                                        <  1:708518/93‑1 (MQ=255)
                                               ggAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTAcc                                        >  1:3665531/1‑93 (MQ=255)
                                                gAACCAACGAGCTTGAGAAGCGACGCCGGAAGCGCACTAGCGCCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCt                                       >  2:1463257/1‑93 (MQ=39)
                                                 aaCCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCtt                                      >  2:2548002/1‑93 (MQ=255)
                                                 aaCCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCATCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCtt                                      <  2:3835883/93‑1 (MQ=255)
                                                  aCCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTc                                     <  1:150940/93‑1 (MQ=255)
                                                    cAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTg                                   >  1:2153625/1‑93 (MQ=255)
                                                     aaCGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGc                                  >  1:629825/1‑93 (MQ=255)
                                                     aaCGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGc                                  >  1:525489/1‑93 (MQ=255)
                                                      aCGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGcc                                 >  2:2471195/1‑93 (MQ=255)
                                                       cGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCa                                <  2:3518219/93‑1 (MQ=255)
                                                        gAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAg                               >  2:2963897/1‑93 (MQ=255)
                                                         aGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAgc                              <  1:3619937/93‑1 (MQ=255)
                                                         aGCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAgc                              >  1:1035156/1‑93 (MQ=255)
                                                          gCTTGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAgcg                             >  1:28203/1‑93 (MQ=255)
                                                            ttGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGcc                           >  1:653488/1‑93 (MQ=255)
                                                             tGAGAAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCt                          <  2:1322455/93‑1 (MQ=255)
                                                                gaAGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATCTACCTTCTGCCAGCGCCTGAg                       >  2:1058300/1‑93 (MQ=255)
                                                                 aaGCGACGCCGGACGCGCCCTAGCGGCTACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGc                      >  2:473897/1‑93 (MQ=255)
                                                                  aGCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCa                     >  1:3919628/1‑93 (MQ=255)
                                                                   gCGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCaa                    <  1:4577429/93‑1 (MQ=255)
                                                                    cGACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAAt                   >  1:2196453/1‑93 (MQ=255)
                                                                     gACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAAtt                  <  1:4037309/93‑1 (MQ=255)
                                                                     gACGCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAAtt                  <  2:4682924/93‑1 (MQ=255)
                                                                        gCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGc               <  2:1866634/93‑1 (MQ=255)
                                                                        gCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGc               <  1:2106248/93‑1 (MQ=255)
                                                                        gCCGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGc               <  2:4484116/93‑1 (MQ=255)
                                                                         ccGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCt              <  2:2067458/93‑1 (MQ=255)
                                                                         ccGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCt              <  2:2163287/93‑1 (MQ=255)
                                                                         ccGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCt              <  1:2148034/93‑1 (MQ=255)
                                                                         ccGGACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCt              >  1:2474224/1‑93 (MQ=255)
                                                                           ggACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTc            <  2:783208/93‑1 (MQ=255)
                                                                           ggACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCACCGCCTGAGCAATTGGCTTc            >  2:4256832/1‑93 (MQ=255)
                                                                            gACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTCg           >  1:3380523/1‑93 (MQ=255)
                                                                            gACGCGCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTCg           >  2:1099392/1‑93 (MQ=255)
                                                                              cgcgCCCTAGCGGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTCGgt         >  2:3357158/1‑93 (MQ=255)
                                                                                         |                                                                                        
AATCAGCGGTTTTTGTGGAATGGCAAGGTTGTTAGATAAGATGCCGTGGAACCAACGAGCTTGAGAAGCGACGCCGGACGCGCCCTAGCAGCGACATCCGGCCTCAGTAATTAGATCAGGAAATCGTCGAGAGATTTACCTTCTGCCAGCGCCTGAGCAATTGGCTTCGGTGTACGGC  >  NC_000913/2796002‑2796179

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