Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,2460GA31.8% 361.8 / 196.5 258noncoding (2543/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (84/92);  new base A (20/62);  total (104/154)
Fisher's exact test for biased strand distribution p-value = 3.80e-04
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 261 total aligned reads displayed.

GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946148‑3946326
                                                                                                  |                                                                                
gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAg                                                                                  <  2:1971073/99‑1 (MQ=34)
  gATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAggg                                                                                <  2:38556/99‑1 (MQ=34)
      aCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATg                                                                            >  1:2172539/1‑99 (MQ=33)
      aCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATg                                                                            <  2:1271810/99‑1 (MQ=21)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                           <  1:685535/99‑1 (MQ=34)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                           <  2:590310/99‑1 (MQ=33)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                           <  1:3819727/99‑1 (MQ=34)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATgg                                                                           <  2:1752650/99‑1 (MQ=21)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGCCCCAAAGGTATgg                                                                           <  2:2348855/99‑1 (MQ=255)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGc                                                                          <  1:599088/99‑1 (MQ=34)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGc                                                                          <  2:713205/99‑1 (MQ=21)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGc                                                                          <  2:1348708/99‑1 (MQ=18)
          gCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTg                                                                        >  1:1546601/1‑99 (MQ=34)
          gCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGGAGGTCCCAAGGGTATGGCTg                                                                        >  1:1136419/1‑99 (MQ=18)
           cTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCAGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGt                                                                       <  2:3724504/99‑1 (MQ=21)
             gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTc                                                                     <  2:2271997/99‑1 (MQ=33)
             gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTc                                                                     <  2:794432/99‑1 (MQ=21)
             gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTc                                                                     <  2:1595871/99‑1 (MQ=18)
              aTACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                    >  2:2728555/1‑99 (MQ=33)
               tACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGc                                                                   <  2:348147/99‑1 (MQ=21)
               tACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGGTCGc                                                                   <  2:4018897/99‑1 (MQ=255)
                aCCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCACATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGcc                                                                  <  1:1098596/99‑1 (MQ=18)
                 ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCa                                                                 <  2:3058862/99‑1 (MQ=34)
                 ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCa                                                                 <  2:3187058/99‑1 (MQ=21)
                 ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCa                                                                 <  2:2649441/99‑1 (MQ=21)
                  cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                >  1:1824094/1‑99 (MQ=34)
                  cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                <  1:2976694/99‑1 (MQ=34)
                  cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTAAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                >  2:3629421/1‑99 (MQ=21)
                      cAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCGTGTAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaa                                                            >  2:4077166/1‑99 (MQ=16)
                       aaGTGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:329166/99‑1 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGCTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:4205287/99‑1 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACTTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:1503552/99‑1 (MQ=16)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:2027103/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:1692032/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  1:2729290/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:3101103/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGCGCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                           <  2:2575662/99‑1 (MQ=16)
                       aaGAGTTCATATCGACAGCGGTGTTGGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGAAGGTCCCAAGGGTATGGCTGTTCGCCATTGaaa                                                           <  2:1257008/99‑1 (MQ=11)
                        agagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:2140092/99‑1 (MQ=34)
                        agaattCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                          <  2:2215591/95‑1 (MQ=255)
                         gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGt                                                         <  2:3699319/99‑1 (MQ=33)
                         gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGt                                                         <  2:2094602/99‑1 (MQ=18)
                             tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                     >  1:1994/1‑99 (MQ=33)
                             tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                     >  1:3132621/1‑99 (MQ=34)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 >  1:1436971/1‑99 (MQ=33)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 <  1:3161291/99‑1 (MQ=33)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 >  1:3925572/1‑99 (MQ=34)
                                   cGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAg                                               >  1:4126709/1‑99 (MQ=34)
                                     aCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                             >  1:2154296/1‑99 (MQ=18)
                                     aCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                             >  1:2315166/1‑99 (MQ=21)
                                     aCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                             >  1:918990/1‑99 (MQ=18)
                                      cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATGTAAAGTGGTACGCGAGCTg                                            <  1:3385511/99‑1 (MQ=21)
                                       ggcggTGTTTGGCACCTTGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           >  1:1657458/1‑99 (MQ=255)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  1:2625694/99‑1 (MQ=34)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  1:4330421/99‑1 (MQ=34)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  2:3804838/99‑1 (MQ=34)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          <  2:1564324/99‑1 (MQ=34)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  1:28672/1‑99 (MQ=18)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          <  2:232428/99‑1 (MQ=21)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  1:3504514/1‑99 (MQ=21)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  1:1264575/1‑99 (MQ=21)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  1:2553441/1‑99 (MQ=21)
                                         cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGGAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                         >  1:1726675/1‑99 (MQ=255)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  2:1073969/99‑1 (MQ=33)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  1:77977/99‑1 (MQ=33)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  1:253920/99‑1 (MQ=34)
                                            tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTa                                      >  1:2827177/1‑99 (MQ=34)
                                             gtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTACGTCCCATGGGTATGGCTGTTCGCCATTAAAAGTGATACGGGAGCTGGGTTTAg                                     <  2:3993012/99‑1 (MQ=11)
                                              tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                    <  1:4238482/99‑1 (MQ=33)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                               <  1:2485395/99‑1 (MQ=18)
                                                    aCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                              >  1:4195810/1‑99 (MQ=34)
                                                    aCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                              >  1:3351601/1‑99 (MQ=33)
                                                     ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTCTGGCTGTTCGCCATATAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                             <  2:82746/99‑1 (MQ=16)
                                                     ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                             <  2:2525944/99‑1 (MQ=34)
                                                      cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            >  2:2889219/1‑99 (MQ=34)
                                                       tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                           >  1:1892913/1‑99 (MQ=34)
                                                         gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgaga                         <  2:2598574/99‑1 (MQ=34)
                                                         gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgaga                         <  2:2511860/99‑1 (MQ=34)
                                                          aTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGAc                        <  1:548905/99‑1 (MQ=34)
                                                          aTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGAc                        <  2:11320/99‑1 (MQ=18)
                                                            gTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAg                      >  1:3703267/1‑99 (MQ=34)
                                                             gcGGCTCACCGCCTCCTGGAGCTCTAGTAAGTCCCAGGGATATTGCTATTTGCCATTGAAAGTGGTACGCGAGCTTTGTTTATAACGTCGTCAGACAGt                     <  2:1933045/98‑1 (MQ=2)
                                                               ggCTCATCACATCCTGGAGCTGGAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   <  2:2985112/99‑1 (MQ=21)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  1:3240898/1‑99 (MQ=34)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  1:4170777/1‑99 (MQ=34)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  1:3287565/1‑99 (MQ=34)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCGTTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   <  2:383513/99‑1 (MQ=255)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  >  2:628489/1‑99 (MQ=33)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  <  2:3257233/99‑1 (MQ=21)
                                                                 cTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCACCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCgg                 <  2:700043/99‑1 (MQ=21)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt                <  1:3781922/99‑1 (MQ=34)
                                                                    atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              >  2:1196395/1‑99 (MQ=34)
                                                                     tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc             <  1:2064589/99‑1 (MQ=34)
                                                                      cacaTCCTGGAGCTGAAGTAGGTCGCAAGGGTATGGCTGTTCGCCATCTAAAGTGGTACGCGAGCTGGGATTAGAACGTCGTGAGACAGGTCGGTCCCt            >  2:955012/1‑99 (MQ=11)
                                                                       acaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTa           >  1:2959569/1‑99 (MQ=34)
                                                                        caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTAt          >  2:3909937/1‑99 (MQ=33)
                                                                         aTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATc         >  1:4286114/1‑99 (MQ=33)
                                                                          tCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCt        >  1:1594301/1‑99 (MQ=34)
                                                                          tCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCt        <  2:539432/99‑1 (MQ=18)
                                                                               gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCTGTCCCTATCTGCCg  <  2:4129887/99‑2 (MQ=17)
                                                                                                  |                                                                                
GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946148‑3946326

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