Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*AC_0000913,543,5470AC38.5% 188.8 / 35.9 104Y102S (TAC→TCC) yiiLL‑rhamnose mutarotase
Reads supporting (aligned to +/- strand):  ref base A (56/8);  new base C (0/40);  total (56/48)
Fisher's exact test for biased strand distribution p-value = 3.79e-21
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 125 total aligned reads displayed.

TTGCCAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACATCTCCCG  >  AC_000091/3543458‑3543635
                                                                                         |                                                                                        
ttGCCAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTa                                                                                          <  2:2779638/90‑1 (MQ=255)
 tGCCAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTAc                                                                                         <  1:2379534/90‑1 (MQ=255)
   ccAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCt                                                                                       <  1:253801/90‑1 (MQ=255)
   ccAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCt                                                                                       >  1:1325043/1‑90 (MQ=255)
    cagcGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTg                                                                                      <  1:793848/87‑1 (MQ=255)
     aaCGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGc                                                                                     >  2:387751/1‑90 (MQ=255)
      aCGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGcc                                                                                    <  2:2464294/90‑1 (MQ=255)
      aCGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGcc                                                                                    <  2:111350/90‑1 (MQ=255)
       cGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCg                                                                                   >  1:3098835/1‑90 (MQ=255)
       cGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCg                                                                                   <  1:2154197/90‑1 (MQ=255)
        gTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGt                                                                                  >  1:3262151/1‑90 (MQ=255)
           ggtggAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAAt                                                                               >  1:990762/1‑90 (MQ=255)
           ggtggAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAAt                                                                               >  1:2151817/1‑90 (MQ=255)
            gtggAAATATATGACCGATGTTATGCCCGCTACCCCGGATAACAGCCGGGTGAGTGGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAAtt                                                                              <  1:934400/90‑1 (MQ=255)
             tggAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTg                                                                             >  2:3029059/1‑90 (MQ=255)
             tggAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTg                                                                             >  1:1324524/1‑90 (MQ=255)
               gAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGcc                                                                           <  1:1905826/90‑1 (MQ=255)
               gAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGcc                                                                           >  2:450130/1‑90 (MQ=255)
                aaaTATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGccc                                                                          <  1:3152932/90‑1 (MQ=255)
                aaaTATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGccc                                                                          >  2:389181/1‑90 (MQ=255)
                 aaTATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCa                                                                         >  2:3431185/1‑90 (MQ=255)
                 aaTATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGAGAGTAGCGAGCTGCAAGACGTGTTTTACCTCCCGTAATTGCCCa                                                                         >  2:1699081/1‑90 (MQ=255)
                  atatatGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCaa                                                                        >  2:3153026/1‑90 (MQ=255)
                     tatGACCGATGTTATGCCCGCTACCCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAAtgt                                                                     <  1:456137/90‑1 (MQ=255)
                     tatGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAAtgt                                                                     >  1:2235254/1‑90 (MQ=255)
                      atGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAAtgtg                                                                    >  1:2901478/1‑90 (MQ=255)
                      atGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAAtgtg                                                                    >  2:3313452/1‑90 (MQ=255)
                      atGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGAAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAAtgtg                                                                    >  2:2706791/1‑90 (MQ=255)
                       tGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGc                                                                   >  1:451638/1‑90 (MQ=255)
                       tGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGc                                                                   >  2:2510845/1‑90 (MQ=255)
                         aCCGATGTTATGCCCGCTACCCCGGATAACAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCa                                                                 <  1:882762/90‑1 (MQ=255)
                          ccGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCag                                                                >  1:1789512/1‑90 (MQ=255)
                              ttttATGCCCGCTACCCCGGATAACAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGAt                                                            <  1:1950691/88‑1 (MQ=255)
                              tGTTATGCCCGCTACCCCGGATAACAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGAt                                                            <  1:2728705/90‑1 (MQ=255)
                              tGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGAt                                                            >  2:2430208/1‑90 (MQ=255)
                                 tATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGcc                                                         >  1:1024059/1‑90 (MQ=255)
                                   tGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGCTGCCTc                                                       >  2:727984/1‑90 (MQ=255)
                                    gCCCGCTACCCCGGATAACAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGATTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTcc                                                      <  1:888213/90‑1 (MQ=255)
                                    gCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTcc                                                      >  2:288955/1‑90 (MQ=255)
                                        gCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGca                                                  >  1:2244733/1‑90 (MQ=255)
                                        gCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGca                                                  >  1:1843419/1‑90 (MQ=255)
                                         cTAACCCGGAAAACACCCGGGTGAGTAGCGAGCTGCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGcac                                                 <  2:1966113/90‑1 (MQ=255)
                                          tAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGcaca                                                >  1:2894324/1‑90 (MQ=255)
                                          tAACCCGGATAACAGCCCGGGGAGTAGCGAGCTGCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGcaca                                                <  1:3449198/90‑1 (MQ=255)
                                            ccccGGATAACAGCCGGGTGAGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAat                                              <  1:1707191/89‑1 (MQ=255)
                                             cccGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAata                                             <  1:1001549/90‑1 (MQ=255)
                                              ccGGATAACAGCCGGGTGAGTAGCGAGCGGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAatat                                            <  2:2847287/90‑1 (MQ=255)
                                               cGGATACAACCCGGGGGAGGAGGGAGCTGAAAAAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATAtt                                           <  2:282698/90‑1 (MQ=255)
                                               cGGAAAACAGCCGGGTGAGTGGGGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATAtt                                           <  1:815032/90‑1 (MQ=255)
                                                ggATAACAGCCGGGTGAGTAGGGAGCTGCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATAttt                                          <  2:3435356/90‑1 (MQ=255)
                                                ggATAACAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATAttt                                          <  1:595698/90‑1 (MQ=255)
                                                ggATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATAttt                                          >  1:2951038/1‑90 (MQ=255)
                                                  aTAACACCCGGGGGAGTAGCGAGGGGCAAGAATTGTTTCCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTa                                        <  2:2678772/90‑1 (MQ=255)
                                                   tAACCGCCGGGTGAGGAGCGAGCTGAAAGAGGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTaa                                       <  1:2800674/90‑1 (MQ=255)
                                                   aaaCAGCCCGGGAATTAGCGGTCTGCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTaa                                       <  2:3261491/89‑1 (MQ=255)
                                                    aaCAGCCTGGTGATTCGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTaaa                                      <  2:3211653/90‑1 (MQ=255)
                                                    aaCAGCCGGGTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTaaa                                      <  1:584294/90‑1 (MQ=255)
                                                     aCAGCCCGGCGGAAAGCGTCCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAt                                     <  2:2095738/90‑1 (MQ=255)
                                                      cAGCCGGGGGAGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAtt                                    <  1:3448898/90‑1 (MQ=255)
                                                      cAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAtt                                    >  2:3163747/1‑90 (MQ=255)
                                                      cAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAtt                                    >  2:3041422/1‑90 (MQ=255)
                                                       aGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAttt                                   >  2:1588345/1‑90 (MQ=255)
                                                       aGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAttt                                   >  2:1805854/1‑90 (MQ=255)
                                                       aGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAAttt                                   >  2:2481943/1‑90 (MQ=255)
                                                        gCCGGGTGAGTAGCGAGCTACAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTc                                  <  2:2945410/90‑1 (MQ=255)
                                                        gCCGGGTAAGTAGGGGGTTCCAAGAGTTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTc                                  <  1:3342317/90‑1 (MQ=255)
                                                        gCCCGGTGGGGAGCGACCGGCAAGAATTGTTTCACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTc                                  <  1:3326711/90‑1 (MQ=255)
                                                          cgggTGGGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTCAATTTCGa                                <  2:2821739/88‑1 (MQ=255)
                                                          ccGGTGAGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGa                                <  2:1613225/90‑1 (MQ=255)
                                                           gggTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGAt                               <  1:567531/89‑1 (MQ=255)
                                                           gggTGAGAAGCGAGTTTCAGGAAGTTTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGAt                               <  1:648348/89‑1 (MQ=255)
                                                           cGGGGAGTAGGGCGCTGTAAGAGTCTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGAt                               <  1:611637/90‑1 (MQ=255)
                                                            ggTGAGTAGCGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGAtt                              <  1:2694077/90‑1 (MQ=255)
                                                             gTGAGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGAttt                             <  1:2636403/90‑1 (MQ=255)
                                                               gggTAGCGAGCTTAAAAAGGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGa                           <  1:1921272/88‑1 (MQ=255)
                                                                aGTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGaa                          <  1:49098/90‑1 (MQ=255)
                                                                 gTAGGGAGCTGCAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGaaa                         <  2:397839/90‑1 (MQ=255)
                                                                 gTAGCGAGCTGAAAGAAGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGaaa                         <  1:1617827/90‑1 (MQ=255)
                                                                   aGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGaaaaa                       >  2:1513803/1‑90 (MQ=255)
                                                                      gAGCTGGAAAATTTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTa                    <  2:1657865/90‑1 (MQ=255)
                                                                      gAGCTGCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTa                    <  2:359145/90‑1 (MQ=255)
                                                                      gAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTa                    >  1:2700615/1‑90 (MQ=255)
                                                                       aGCTGCAAGAGGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAACTTCGATTTGAAAAATTaa                   <  1:474567/90‑1 (MQ=255)
                                                                       aGCTGCAAGAGGTGTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTaa                   <  2:691317/90‑1 (MQ=255)
                                                                         cTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAAc                 >  1:1151205/1‑90 (MQ=255)
                                                                          tGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACa                >  1:2662074/1‑90 (MQ=255)
                                                                          tGCAAGAAGGGTTTTACCTGCCGTAATTGCCCAATGTACCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACa                >  1:815247/1‑90 (MQ=255)
                                                                           gCAAGAAGTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACat               <  1:1042294/90‑1 (MQ=255)
                                                                           gCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACat               >  2:2603402/1‑90 (MQ=255)
                                                                            tgagatgTATTTTCCGTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACata              <  2:2983367/88‑1 (MQ=255)
                                                                             aagaagTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACatat             >  2:3166983/1‑90 (MQ=255)
                                                                             aagaagTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACatat             >  1:3348184/1‑90 (MQ=255)
                                                                                 agTTTTTTCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACa         <  1:1049254/90‑1 (MQ=255)
                                                                                   tgtTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtc       >  1:3431841/1‑90 (MQ=255)
                                                                                   tgtTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtc       >  2:2087268/1‑90 (MQ=255)
                                                                                    tttttGCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtct      <  1:1112802/89‑1 (MQ=255)
                                                                                    tttttCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtct      <  1:1476855/89‑1 (MQ=255)
                                                                                    tttttACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtct      <  2:549790/89‑1 (MQ=255)
                                                                                     ttttCCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACAtctc     <  2:1014931/90‑1 (MQ=255)
                                                                                       tcaCCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACATCTccc   <  2:1480142/88‑1 (MQ=255)
                                                                                         |                                                                                        
TTGCCAACGTTGGTGGAAATATATGACCGATGTTATGCCCGCTAACCCGGATAACAGCCCGGTGAGTAGCGAGCTGCAAGAAGTGTTTTACCTGCCGTAATTGCCCAATGTGCCAGAGATGCCTCCGGCACAATATTTTAAATTTCGATTTGAAAAATTAAACATATTACATCTCCCG  >  AC_000091/3543458‑3543635

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