Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,831 C→T 15.3% intergenic (+23/+17) tonB → / ← yciA membrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,8310CT15.3% 216.8 / 22.5 94intergenic (+23/+17)tonB/yciAmembrane spanning protein in TonB‑ExbB‑ExbD transport complex/acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base C (30/48);  new base T (8/7);  total (38/56)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGT  >  NC_000913/1311736‑1311929
                                                                                               |                                                                                                  
gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                >  1:2795014/1‑100 (MQ=255)
gCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGa                                                                                                >  1:1031399/1‑100 (MQ=255)
 ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCTGAAAGTCAAAAGCCTCCGACCGGAg                                                                                               <  2:220829/100‑1 (MQ=255)
 ccGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                               <  1:2066777/100‑1 (MQ=255)
   gggTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAg                                                                                               <  1:743091/98‑1 (MQ=255)
      tAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCttt                                                                                          <  2:2018346/100‑1 (MQ=255)
        tccgatctaGTGGGATTGTGGTGAATATCCTGTTTACAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                        <  1:2135903/92‑1 (MQ=255)
        aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTg                                                                                        <  2:1083785/100‑1 (MQ=255)
        aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                        <  2:2522355/100‑1 (MQ=255)
        aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                        <  2:168506/100‑1 (MQ=255)
        aGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTg                                                                                        <  1:2086763/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGa                                                                                       <  1:281111/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                       <  1:463227/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                       <  2:1105950/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                       <  2:2210031/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGa                                                                                       <  1:1076466/100‑1 (MQ=255)
         gCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGACTTTTGa                                                                                       >  1:806767/1‑100 (MQ=255)
          ccAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACAGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                      <  1:495876/100‑1 (MQ=255)
          ccAGGCAGTGGGATTGTGGTGAATATACTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                      <  2:1695057/100‑1 (MQ=255)
           cAGGCAGTGGGATTGTGGTGAATATCCTGTTTATAATTAACGGCTCCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACt                                                                                     <  1:2351189/100‑1 (MQ=38)
           cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGAc                                                                                      <  1:1213517/99‑1 (MQ=255)
           cAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACt                                                                                     <  2:1978233/100‑1 (MQ=255)
            aGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACt                                                                                     <  2:1462952/99‑1 (MQ=255)
            aGGCAGTGGGATTGTAGTGAATATCCTGTGTAAAATTAACGGCACCACCGACATTGAGTAAGCAGAAAGTCAAAAGTTTCCGGTCGGAGGCTTTTGACt                                                                                     <  2:1392220/99‑1 (MQ=37)
             ggCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTAt                                                                                   <  2:1031392/100‑1 (MQ=255)
              gCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGGTCGGAGGCTTTTGACtttc                                                                                  <  1:88304/100‑4 (MQ=38)
                aGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGagatcgga                                                                                >  2:2135880/1‑93 (MQ=255)
                  tGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTc                                                                              <  2:1957016/100‑1 (MQ=255)
                   gtgATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCa                                                                             <  2:437834/98‑1 (MQ=255)
                   gggATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCa                                                                             <  2:226143/100‑1 (MQ=255)
                     gATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAc                                                                           <  2:2583724/100‑1 (MQ=255)
                      aTTGTAGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACa                                                                          <  2:2698048/100‑1 (MQ=255)
                          tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTa                                                                      <  2:242869/100‑1 (MQ=255)
                          tggtgAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAAGCTTTTGACTATTACTCAACAGGTa                                                                      >  1:1158218/1‑100 (MQ=255)
                                  atCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGa                                                              <  1:2820346/100‑1 (MQ=255)
                                  atCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGa                                                              >  1:2259549/1‑100 (MQ=255)
                                    ccgatctAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAgg                                                            <  1:2911755/94‑1 (MQ=255)
                                     cTGTTTAAAATTAACGGTACCACCGAAATTTAGTAAGTAGAAAGTCAAAAGCCTCCGAACGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGt                                                           <  1:203559/100‑1 (MQ=39)
                                          tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCTACCGGAGGCTTTTGACTATTACTCAACAGAAAAGGCGCGAGGTTTTcc                                                      >  1:1834950/1‑100 (MQ=255)
                                          tAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGAAAGGCGCGAGGTTTTcc                                                      >  1:1417320/1‑100 (MQ=255)
                                           aaaaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAggagatcgg                                                     >  2:2911723/1‑93 (MQ=255)
                                             aaTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTc                                                   <  2:2555873/100‑1 (MQ=255)
                                              aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                  <  1:2841497/100‑1 (MQ=255)
                                              aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                  <  2:1502560/100‑1 (MQ=255)
                                              aTTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAAAAGGTAAGGCGCGAGGTTTTACtt                                                    >  2:3014566/1‑98 (MQ=255)
                                               ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGTTAAGGCGCGAGGTTTTCCTTCAg                                                 <  2:693007/100‑1 (MQ=255)
                                               ttAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCa                                                  >  2:1268852/1‑99 (MQ=255)
                                                tAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAgg                                                <  2:626172/100‑1 (MQ=255)
                                                     gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCaa                                           >  2:2002001/1‑100 (MQ=255)
                                                     gCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCa                                            <  2:1107582/99‑1 (MQ=255)
                                                      caccacCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTAAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAc                                          >  1:1820429/1‑100 (MQ=255)
                                                       accaccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAAcc                                         <  2:1595798/100‑1 (MQ=255)
                                                         caccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACcgc                                       <  2:2547499/100‑1 (MQ=255)
                                                         caccGAAATTCAGTAAGCAGAAACTCAAAAGTCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTACTTCAGGATCaa                                           <  2:566208/96‑1 (MQ=255)
                                                           ccGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGATTTCCTTCAGGATCAACCGCGa                                     >  1:813938/1‑100 (MQ=255)
                                                             gagaTTCAGTAAGCAGAAAATCAAAAGCCTCCGACCTGACGCTTTTGACTATTACTCAACAGGTAACACGCGAGGCTTTCCTTCAGGATCAACCGCGACa                                   >  2:1359587/4‑100 (MQ=255)
                                                             gAAATTCAGTAGGCAGAAAGTCAAAAGCCGGCGACCGGAGGCTTTTGACTATTACTCAGCAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACa                                   <  2:2353435/100‑1 (MQ=255)
                                                               aaTTCAGTAAGCAGAATGTCAAAAGCCTCCGACCGGAGGCTTTTGAATATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATa                                 <  2:1765638/100‑1 (MQ=255)
                                                                aTTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATAc                                >  1:445212/1‑100 (MQ=255)
                                                                 ttCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACt                               >  2:1671730/1‑100 (MQ=255)
                                                                      tAAGCAGAAATTCAATAGCCTCCGATCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGAGCAACCGCGACATACTTaaa                           <  1:262600/99‑1 (MQ=255)
                                                                      tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                          >  1:830164/1‑100 (MQ=255)
                                                                      tAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaat                          <  2:1085619/100‑1 (MQ=255)
                                                                        accAGAGAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGTGAGGTTTTCCGTCAGGATCATCCGCGACATACTTAaataa                        <  2:70070/98‑1 (MQ=39)
                                                                        aGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaata                         >  2:2403588/1‑99 (MQ=255)
                                                                         gCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTTACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                       <  2:1595351/100‑1 (MQ=255)
                                                                         gCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCAAGGTTTTCCTTCAGGATCAACCGCGACATACTTAaataat                       >  1:501990/1‑100 (MQ=255)
                                                                           aGAAAGTCAAAAGCCTCCGACCGGAGTCTTTTGATTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGc                     <  1:272544/100‑1 (MQ=255)
                                                                            gAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCt                    >  1:2483986/1‑100 (MQ=255)
                                                                             aaaGTCAAAAGGCGCCGATCGGAGGCTTTTGACTAGTATTCAACAGTTAAGTCGCGAGGGTCTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt                   <  2:1659594/100‑1 (MQ=255)
                                                                             aaaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCtt                   <  2:2727771/100‑1 (MQ=255)
                                                                              aaGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTc                  <  2:2234246/100‑1 (MQ=255)
                                                                                gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTg                >  2:712820/1‑100 (MQ=255)
                                                                                gTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAAGACCAACCAAGACAGACATAAATAATGCTTCTg                >  1:1111009/1‑100 (MQ=255)
                                                                                 tCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGt               >  2:1436852/1‑100 (MQ=255)
                                                                                  cAAAAGCCTCCGGTCGGAGGCTTTTGACTTTTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc              >  1:251379/1‑100 (MQ=38)
                                                                                  cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc              >  2:1513030/1‑100 (MQ=255)
                                                                                  cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc              >  2:1893930/1‑100 (MQ=255)
                                                                                  cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTc              >  2:2629286/1‑100 (MQ=255)
                                                                                  cAAAAGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACAGCGACATACTTAAATAATGATTCTGTc              >  1:2984322/1‑100 (MQ=38)
                                                                                  cAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCCCGACATACTTAAATAATGCTTCt                 >  1:2280834/1‑97 (MQ=255)
                                                                                   aaaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg             >  2:1850456/1‑100 (MQ=255)
                                                                                   aaaaGCCTCCGGTCGGAGGATTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTATg                >  2:114291/1‑97 (MQ=39)
                                                                                   aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg             <  1:1889964/100‑1 (MQ=255)
                                                                                   aaaaGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCg             >  1:2918585/1‑100 (MQ=255)
                                                                                    aaaGCCTCCGGTCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTccc            >  2:2395819/1‑98 (MQ=255)
                                                                                      aGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCtt          <  2:1959657/100‑1 (MQ=255)
                                                                                       gCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCttt         >  1:621633/1‑100 (MQ=255)
                                                                                        ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTta        <  2:505456/100‑1 (MQ=255)
                                                                                        ccTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGAATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTta        >  1:309493/1‑100 (MQ=255)
                                                                                         cTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTtat       <  2:2058225/100‑1 (MQ=255)
                                                                                          nccGACCGGAGGCTTTTGACGATTACACAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTACTGTCGCTTta        >  2:1848555/2‑98 (MQ=255)
                                                                                           ccGACCGGAGGCTTTTGAATATTACTCAACAGGTAAAGCACGAGGTTTTCATTCAGGATCAACCGCGACATAATTAAATAATGCTTCTGTCGCTTtat       >  2:1408443/1‑98 (MQ=255)
                                                                                              aCCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGt  <  1:1212177/100‑1 (MQ=255)
                                                                                               |                                                                                                  
GCCGGGTAAGCCAGGCAGTGGGATTGTGGTGAATATCCTGTTTAAAATTAACGGCACCACCGAAATTCAGTAAGCAGAAAGTCAAAAGCCTCCGACCGGAGGCTTTTGACTATTACTCAACAGGTAAGGCGCGAGGTTTTCCTTCAGGATCAACCGCGACATACTTAAATAATGCTTCTGTCGCTTTATAGCGT  >  NC_000913/1311736‑1311929

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