Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,2460GA45.0% 174.8 / inf 292noncoding (2543/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (73/87);  new base A (34/97);  total (108/184)
Fisher's exact test for biased strand distribution p-value = 6.11e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.55e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 297 total aligned reads displayed.

GGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGG  >  NC_000913/3946147‑3946327
                                                                                                   |                                                                                 
 gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAgg                                                                                  <  1:989825/100‑1 (MQ=34)
   gATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGt                                                                                >  1:1113840/1‑100 (MQ=34)
     tAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAt                                                                              >  2:2457815/1‑100 (MQ=33)
      aaCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATg                                                                             <  1:2032826/100‑1 (MQ=34)
          ggCTGATACCGCCCAAGAGTTCATATCGCCGGCGGTGTTTGGCACCTCGATGTCGGCTCCTCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTg                                                                         <  2:1275331/100‑1 (MQ=11)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTg                                                                         >  2:3577905/1‑100 (MQ=34)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTg                                                                         <  1:3042646/100‑1 (MQ=34)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTg                                                                         <  2:502473/100‑1 (MQ=21)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTg                                                                         <  2:3256205/100‑1 (MQ=18)
           gCTGATACCGACCCAGAGTTCATATCGCCGGCGGTGTTTGGCGCCTCGATGTCGGCTCATCACCTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGt                                                                        <  2:417694/100‑1 (MQ=11)
              gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                     <  1:1468491/100‑1 (MQ=34)
               aTACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGc                                                                    >  2:3335700/1‑100 (MQ=34)
                  ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                 >  2:2020313/1‑100 (MQ=34)
                  ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                 >  2:2048695/1‑100 (MQ=34)
                   cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAtt                                                                <  2:4089410/100‑1 (MQ=34)
                    gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAttt                                                               >  1:819523/1‑100 (MQ=34)
                    gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAttt                                                               >  2:3831966/1‑100 (MQ=33)
                      ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaa                                                             <  2:2232738/100‑1 (MQ=18)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:1779282/100‑1 (MQ=33)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:211457/100‑1 (MQ=21)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:19841/100‑1 (MQ=21)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:152246/100‑1 (MQ=18)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:1603304/100‑1 (MQ=21)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:2478881/100‑1 (MQ=21)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:1119470/100‑1 (MQ=18)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:329375/100‑1 (MQ=21)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACACCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:165407/100‑1 (MQ=255)
                        aaGAGTTCATATCGACGGCGGTATTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                           <  2:2432715/100‑1 (MQ=18)
                        aaGAGGTCAGATCGACGGCGGTGTTTGGCACCTCGATGTGGGCTCATCCCATCCTGGAGCTGAAGTAGGTCCCAAAGGTCTGGCTGTTCGCCATTTACAg                                                           <  2:837405/100‑1 (MQ=11)
                          gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTg                                                         <  2:3096287/100‑1 (MQ=18)
                           agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTgg                                                        <  2:2281563/100‑1 (MQ=34)
                           agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTgg                                                        <  2:3013697/100‑1 (MQ=21)
                           agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTgg                                                        <  2:1395496/100‑1 (MQ=18)
                              tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAc                                                     <  2:2624721/100‑1 (MQ=33)
                              tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTAc                                                     <  2:2851445/100‑1 (MQ=21)
                               cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                    <  2:1295070/100‑1 (MQ=34)
                               cATACCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                    <  2:2517414/100‑1 (MQ=255)
                                   tCGATGGCGGTGTTTGGCGCCTCGATGTCGGGTCATCAGCTCCTGGAGCTGAGGTAGGTCCGAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAg                                                <  2:2746382/100‑1 (MQ=11)
                                   tCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAg                                                <  2:2674653/100‑1 (MQ=33)
                                   tCGACGGCGGTGTTTGGCACCTCGATGTCGGCGCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAg                                                >  2:906328/1‑100 (MQ=21)
                                    cGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAACGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                               >  1:4154176/1‑100 (MQ=21)
                                     gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                              <  1:2376480/100‑1 (MQ=34)
                                     gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                              >  1:4093536/1‑100 (MQ=21)
                                      tcggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                             >  1:2596592/2‑100 (MQ=16)
                                      aCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                             <  1:2901166/100‑1 (MQ=34)
                                       cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                            >  1:3613175/1‑100 (MQ=18)
                                       cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                            >  1:128142/1‑100 (MQ=21)
                                       cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                            >  1:2974810/1‑100 (MQ=21)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                           <  2:3728685/100‑1 (MQ=34)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                           >  1:156883/1‑100 (MQ=21)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                           >  1:2145612/1‑100 (MQ=21)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                           <  2:1848409/100‑1 (MQ=21)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATATAAAGTGGTACGCGAGCTggg                                           >  1:3651980/1‑100 (MQ=255)
                                          cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                         <  1:3578756/100‑1 (MQ=34)
                                          cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                         <  1:1353435/100‑1 (MQ=34)
                                          cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                         >  2:3225349/1‑100 (MQ=34)
                                          cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                         <  2:1953118/100‑1 (MQ=34)
                                           ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                        >  1:2496897/1‑100 (MQ=34)
                                           ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGCCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                        <  2:1172617/100‑1 (MQ=255)
                                            gtgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTa                                       <  1:3838872/100‑1 (MQ=34)
                                              gtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                     <  1:3280623/100‑1 (MQ=34)
                                              gtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                     <  2:2732914/100‑1 (MQ=18)
                                                   gCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                                >  1:4267828/1‑100 (MQ=34)
                                                    cACTTCGATGTCGGCTCATAACATCCTGGAGCTGAAGTAGGTCCAAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                               <  2:2552055/100‑1 (MQ=11)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                               <  1:1595068/100‑1 (MQ=34)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                               >  1:1690636/1‑100 (MQ=34)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                               <  2:540657/100‑1 (MQ=33)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                               <  2:696230/100‑1 (MQ=18)
                                                     aCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                              <  2:3773674/100‑1 (MQ=34)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                             >  1:2865903/1‑100 (MQ=21)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                             >  1:1815159/1‑100 (MQ=21)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                             >  1:582872/1‑100 (MQ=21)
                                                       cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                            <  2:2572990/100‑1 (MQ=33)
                                                        tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgan                           <  1:3513712/100‑2 (MQ=21)
                                                        tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgag                           <  1:1035148/100‑1 (MQ=33)
                                                         cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgaga                          <  1:2078797/100‑1 (MQ=34)
                                                          gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGAc                         <  2:1247507/100‑1 (MQ=21)
                                                            tGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAg                       >  1:4032919/1‑100 (MQ=34)
                                                            tGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAg                       <  1:4388500/100‑1 (MQ=34)
                                                            tGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAg                       <  1:47558/100‑1 (MQ=33)
                                                            tGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAg                       <  2:3449838/100‑1 (MQ=21)
                                                             gTCGGCGGATAACATCCTGAAGCTGAACTAGGTCTCAAGGGTATGGCTGTTCGCCATTAAAAGTGGTACTCGATCTTGGTTTAGAACGGCGTGAGACGGt                      <  2:176493/100‑1 (MQ=11)
                                                             ctCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGt                      >  1:3123372/2‑100 (MQ=18)
                                                               cGGCTGATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                    <  2:3193656/100‑1 (MQ=18)
                                                               cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                    >  1:1207556/1‑100 (MQ=34)
                                                               cGGCTCAGCACATCGTGGAGCGGAAGTAGGTCCCAAGGGTATTGCTGTTCGCCATATAGAGTCGTACGCGAGCAGGGTTCAGAACGTCGTGAAGCAGTTc                    <  2:4330945/100‑1 (MQ=11)
                                                                ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTAAGACAGTTCg                   >  1:3355570/1‑100 (MQ=21)
                                                                ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                   <  2:3591585/100‑1 (MQ=18)
                                                                ggCTCACCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAAGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGTACGTCGTGAGACAGTTCg                   >  2:4013280/1‑100 (MQ=255)
                                                                   tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc                <  2:105646/100‑1 (MQ=18)
                                                                     ntcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc              >  1:1908686/2‑100 (MQ=18)
                                                                     atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAACTGGTACGCGTGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc              <  2:4232371/100‑1 (MQ=16)
                                                                     atcaCATCCTGGAGCAGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGCGCTGGGGTTAGAACGTCGTGAGACAGTTCGGTccc              >  2:1495704/1‑100 (MQ=255)
                                                                        acaTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTGAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTAt           <  2:3922375/100‑1 (MQ=255)
                                                                         caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATc          <  2:2322193/100‑1 (MQ=34)
                                                                          aTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCt         >  1:4284128/1‑100 (MQ=21)
                                                                           tCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTg        >  2:765632/1‑100 (MQ=33)
                                                                            ccTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc       >  1:2717839/1‑100 (MQ=34)
                                                                            ccTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc       <  2:2955546/100‑1 (MQ=33)
                                                                               ggAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGt    <  1:897392/100‑1 (MQ=33)
                                                                                                   |                                                                                 
GGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTGG  >  NC_000913/3946147‑3946327

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: ‑