Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,946,790 T→G 18.2% G180G (GGT→GGG yebB → DUF830 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,946,7900TG18.2% 212.7 / 51.5 105G180G (GGT→GGGyebBDUF830 family protein
Reads supporting (aligned to +/- strand):  ref base T (36/50);  new base G (9/10);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 7.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.37e-02

Only 100 of 107 total aligned reads displayed.

CCGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAG  >  NC_000913/1946692‑1946884
                                                                                                  |                                                                                              
ccGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGtt                                                                                               <  2:3575935/100‑1 (MQ=255)
 cGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGttt                                                                                              >  1:3388245/1‑100 (MQ=255)
  gAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTg                                                                                             <  2:1133247/100‑1 (MQ=255)
   aaTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTgg                                                                                            <  1:1061835/100‑1 (MQ=255)
    aTGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGgt                                                                                           <  1:557357/100‑1 (MQ=255)
      gCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGgtgt                                                                                         >  2:3112543/1‑100 (MQ=255)
       cAAAACTCACTTTCTGGAAATTCTGTTGCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGtt                                                                                        <  2:3500186/100‑1 (MQ=255)
       cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGtt                                                                                        <  1:3112562/100‑1 (MQ=255)
       cAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGtt                                                                                        <  1:594002/100‑1 (MQ=255)
         aaaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGa                                                                                      >  1:1910172/1‑100 (MQ=255)
          aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGAt                                                                                     <  1:2305045/100‑1 (MQ=255)
          aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGAt                                                                                     <  2:2433735/100‑1 (MQ=255)
          aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGAt                                                                                     >  1:1012516/1‑100 (MQ=255)
          aaCTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATGCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGAt                                                                                     <  2:663530/100‑1 (MQ=255)
           aCTCACTTTCTGGAAATGCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCAGGCCAGACAGTTTGTGGCATCATCCGGGTTTGGTGCTGAt                                                                                     <  2:3937332/99‑1 (MQ=255)
            cTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTc                                                                                   >  2:1223439/1‑100 (MQ=255)
             tCACTTTCTGGAGATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCa                                                                                  <  1:1195142/100‑1 (MQ=255)
             tCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCa                                                                                  <  1:3619763/100‑1 (MQ=255)
             tCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCa                                                                                  >  1:1168805/1‑100 (MQ=255)
              cACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCAt                                                                                 >  1:1322811/1‑100 (MQ=255)
              cACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCAt                                                                                 >  1:1322411/1‑100 (MQ=255)
              cACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCAt                                                                                 >  1:2964744/1‑100 (MQ=255)
                cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGc                                                                               <  1:32272/100‑1 (MQ=255)
                cTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGc                                                                               >  2:767759/1‑100 (MQ=255)
                 tttCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCg                                                                              <  2:1783104/100‑1 (MQ=255)
                   tCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCgga                                                                            <  1:1318839/100‑1 (MQ=255)
                    cTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCggag                                                                           >  2:3105434/1‑100 (MQ=255)
                     tGGAAATTTTGGTTCGTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCTCGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggagg                                                                          <  2:4057784/100‑1 (MQ=255)
                     tGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggagg                                                                          <  1:1812811/100‑1 (MQ=255)
                      ggAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTGGATTCATGCggaggg                                                                         >  2:1663237/1‑100 (MQ=255)
                        aaaTTTTGCTTCCTGGGTTTTATTCCGTAGGTGCGTAAAACCGTCCCGCCAGCCAGTTTGTGGCATAATCCGGGTTTGTTGTTGATTCAGGCggaggt                                                                         <  2:1335465/98‑2 (MQ=255)
                        aaaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggagggag                                                                       >  1:1544652/1‑100 (MQ=255)
                        aaaTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCggagggag                                                                       >  2:1913732/1‑100 (MQ=255)
                          aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTg                                                                     <  2:2427754/100‑1 (MQ=255)
                          aTTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTg                                                                     <  1:1866293/100‑1 (MQ=255)
                           ttCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTgg                                                                    <  1:3558208/100‑1 (MQ=255)
                                gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCCGCCAGTTTGTGGCATCATCCGGCTTTGGGGTTGATTCATGCGGAGGGAGTGGAAAc                                                                >  1:635166/1‑99 (MQ=255)
                                gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                               <  1:1226651/100‑1 (MQ=255)
                                gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                               <  2:67068/100‑1 (MQ=255)
                                gTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACg                                                               <  1:457212/100‑1 (MQ=255)
                                  tCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                             <  1:3477298/100‑1 (MQ=255)
                                   ccTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCt                                                            <  2:2324205/100‑1 (MQ=255)
                                    cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                           >  1:3659502/1‑100 (MQ=255)
                                    cTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTc                                                           >  2:3873003/1‑100 (MQ=255)
                                     tGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCa                                                          >  1:1389543/1‑100 (MQ=255)
                                       ggTTGTATTCCGTGGGAGCGTAAAGCCGTCACGCCAGCCAGTTTGGGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGcc                                                             <  2:2241707/95‑1 (MQ=255)
                                        gTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGcc                                                       <  2:3839823/100‑1 (MQ=255)
                                         ttCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCg                                                      <  2:1644325/100‑1 (MQ=255)
                                         ttCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGCGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCg                                                      <  2:2209337/100‑1 (MQ=255)
                                             aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  >  2:1851486/1‑100 (MQ=255)
                                             aTTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAAc                                                  <  1:2324392/100‑1 (MQ=255)
                                              ttCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACt                                                 <  1:1495446/100‑1 (MQ=255)
                                                ccGTGGGAGCGTAAAACCTTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATTCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGa                                               >  2:401536/1‑100 (MQ=255)
                                                 cGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGAc                                              >  2:1343664/1‑100 (MQ=255)
                                                    gggAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGa                                           <  2:3274593/100‑1 (MQ=255)
                                                    gggAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGa                                           <  2:2189662/100‑1 (MQ=255)
                                                     ggAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAg                                          >  2:4057151/1‑100 (MQ=255)
                                                      gAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                         >  1:588465/1‑100 (MQ=255)
                                                      gAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                         >  1:3828982/1‑100 (MQ=255)
                                                       aGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAgg                                         <  1:1266673/99‑1 (MQ=255)
                                                       aGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGc                                        <  1:48187/100‑1 (MQ=255)
                                                         cGTAAAACCGTCACGCCAGCCAGTTTGTTGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                      <  1:1225488/100‑1 (MQ=255)
                                                         cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                      >  1:4052566/1‑100 (MQ=255)
                                                         cGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCgg                                      <  2:2364238/100‑1 (MQ=255)
                                                          gTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGt                                     >  1:1718612/1‑100 (MQ=255)
                                                           tAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGta                                    <  2:1851289/100‑1 (MQ=255)
                                                             aaaCCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGtata                                  <  1:1343673/100‑1 (MQ=255)
                                                                ccGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACt                               >  2:748970/1‑100 (MQ=255)
                                                                 cGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAAc                                <  2:2962315/98‑1 (MQ=255)
                                                                 cGTCACGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACtt                              <  2:2841595/100‑1 (MQ=255)
                                                                     aCGCCAGCCAGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACg                          <  2:2929321/100‑1 (MQ=255)
                                                                      cGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGc                         >  2:4029439/1‑100 (MQ=255)
                                                                      cGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGc                         <  2:2663247/100‑1 (MQ=255)
                                                                       gccagccaGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                        >  1:3671778/1‑100 (MQ=255)
                                                                       gccagccaGTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCa                        >  1:2410276/1‑100 (MQ=255)
                                                                         cagccagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGt                      >  1:2625880/1‑100 (MQ=255)
                                                                         cagccagGTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAg                       <  2:182001/99‑1 (MQ=255)
                                                                           gccagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCg                    >  1:3120824/1‑100 (MQ=255)
                                                                           gccagTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCg                    >  1:606626/1‑100 (MQ=255)
                                                                              agTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCTGTCGccc                 <  1:2167632/100‑1 (MQ=255)
                                                                              agTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGccc                 <  1:3871861/100‑1 (MQ=255)
                                                                               gTTTGTGGCATCATCCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCct                >  2:260247/1‑100 (MQ=255)
                                                                                tttGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctc               <  1:4104826/100‑1 (MQ=255)
                                                                                 ttGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCctct              <  2:1923152/100‑1 (MQ=255)
                                                                                   gtgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            <  1:505263/100‑1 (MQ=255)
                                                                                   gtgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            >  1:3330880/1‑100 (MQ=255)
                                                                                   gtgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            >  2:2085015/1‑100 (MQ=255)
                                                                                   gtgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCg            <  1:88913/100‑1 (MQ=255)
                                                                                    tgGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGc           <  2:3403194/100‑1 (MQ=255)
                                                                                     ggCATCGTCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGcc          >  1:2165267/1‑100 (MQ=255)
                                                                                        atcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg       <  2:3454100/100‑1 (MQ=255)
                                                                                        atcaacCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAgg       <  2:3665758/100‑1 (MQ=255)
                                                                                         tcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt      >  1:2400306/1‑100 (MQ=255)
                                                                                         tcatcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGt      >  2:2698280/1‑100 (MQ=255)
                                                                                          catcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCGCAGCCGGAGCTGAGCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt     >  2:1482572/1‑100 (MQ=255)
                                                                                          catcCGGGGTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGtt     >  1:1298929/1‑100 (MQ=255)
                                                                                           atcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGCGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc    >  1:1208916/1‑100 (MQ=255)
                                                                                           atcCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTc    >  1:2174644/1‑100 (MQ=255)
                                                                                             ccGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAg  <  2:194772/100‑1 (MQ=255)
                                                                                                 gTTTGGTGTTGATTCATGCGGAGGGAGGGGAAACGCCTCAGCCGGGACTGCGCGAGGCGGGAAAACTTAACCCACTCGCCCGCTGGACAGGTTc    >  2:3699079/1‑94 (MQ=255)
                                                                                                  |                                                                                              
CCGAATGCAAAACTCACTTTCTGGAAATTCTGGTTCCTGGGTTCTATTCCGTGGGAGCGTAAAACCGTCACGCCAGCCAGTTTGTGGCATCATCCGGGTTTGGTGTTGATTCATGCGGAGGGAGTGGAAACGCCTCAGCCGGAACTGACCGAGGCGGTATAACTTAACGCAGTCGCCCTCTCGCCAGGTTCAG  >  NC_000913/1946692‑1946884

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: