Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,2460GA34.1% 237.5 / 173.3 213noncoding (2543/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (69/70);  new base A (19/53);  total (88/125)
Fisher's exact test for biased strand distribution p-value = 1.21e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 224 total aligned reads displayed.

GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946148‑3946326
                                                                                                  |                                                                                 
 ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTTGCACCTCGATGTCGGCTAATCACATCCTGGAGCTGAAGTAGGTCCCAAgg                                                                                  <  1:1533395/99‑1 (MQ=16)
 ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAgg                                                                                  <  1:1509867/99‑1 (MQ=34)
 ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAgg                                                                                  <  1:2399931/99‑1 (MQ=34)
     aaCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAt                                                                              <  2:2153664/99‑1 (MQ=33)
     aaCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAt                                                                              >  1:5016732/1‑99 (MQ=33)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                            <  2:1891723/99‑1 (MQ=33)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATgg                                                                            <  2:4513265/99‑1 (MQ=21)
            tGATACCGCCCAAGAGTTCATCTAGACGGCGGTATTAGGCACCTCGATATCGGCTCAACACATCATGGAGCAGCTGCAGGCCCGAAGGGTACGGCAGtt                                                                      >  2:4360720/1‑99 (MQ=2)
            tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGtt                                                                      <  2:3190784/99‑1 (MQ=33)
            tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGtt                                                                      >  1:2564625/1‑99 (MQ=34)
             gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTc                                                                     >  1:2550696/1‑99 (MQ=34)
             gATACCGCCCAAGAGTTCATAGCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTc                                                                     <  2:4413446/99‑1 (MQ=18)
                aCCGCCCAAGGGTTCATATCGACGGTGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGTTATGGCTGTTCGct                                                                  <  2:1131738/99‑2 (MQ=11)
                aCCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGcc                                                                  <  1:1688418/99‑1 (MQ=33)
                 ccGCCCAAGAGTTCATATCTACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCa                                                                 <  2:4601907/99‑1 (MQ=21)
                 ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCa                                                                 >  2:3672674/1‑99 (MQ=33)
                    cccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAGGTAGGTCCCAAGGGTATGGCTGTTCGCCAttt                                                              >  2:1110361/1‑99 (MQ=18)
                     ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGCATGGCTGTTCGCCATTTa                                                             >  1:4946381/1‑99 (MQ=18)
                      cAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaa                                                            >  2:3612412/1‑99 (MQ=34)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACTTCGATGTCGGCTCAGCACATCCTGGAGCTGAAGCAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:4771619/99‑1 (MQ=11)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGCTGTCGGCTCATCACATCCTGCAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:4910966/99‑1 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                           >  2:5131802/1‑99 (MQ=34)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  1:5278815/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:1421496/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:4067246/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:2677816/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  1:5158117/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:179311/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:403800/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                           <  2:4795129/99‑1 (MQ=21)
                       aaGAGTTCATAACCACGGCGGTGTGTGGATACTTGATATCGGTTTATCGCATCCTGGAGCTGAAGTAGGGCCCAAGTGTTTGACTGTTCGCCATTTaaa                                                           <  2:5161439/99‑1 (MQ=2)
                         gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGt                                                         >  1:913428/1‑99 (MQ=33)
                            ctCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGt                                                      <  2:2446853/98‑1 (MQ=18)
                             tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                     <  2:1842627/99‑1 (MQ=33)
                              cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAc                                                    >  2:646592/1‑99 (MQ=34)
                               atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                   >  1:2898445/1‑99 (MQ=33)
                               atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                   >  1:1343019/1‑99 (MQ=34)
                               atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                   >  2:1254892/1‑99 (MQ=34)
                                tatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTAcgc                                                  <  2:2973918/99‑1 (MQ=18)
                                tatCGACGGCGGGGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAATTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgc                                                  >  2:653944/1‑99 (MQ=255)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 >  1:4135436/1‑99 (MQ=33)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 >  1:1254259/1‑99 (MQ=33)
                                  tCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGa                                                <  1:4409611/99‑1 (MQ=33)
                                  tCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGa                                                >  2:705932/1‑99 (MQ=34)
                                    gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                              >  2:3874111/1‑99 (MQ=34)
                                    gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                              >  1:4062734/1‑99 (MQ=21)
                                    gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                              >  1:194340/1‑99 (MQ=21)
                                     aCGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGAGGTACGCGAGCt                                             <  2:3282960/99‑1 (MQ=21)
                                      cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                            <  1:4207481/99‑1 (MQ=33)
                                      cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                            >  1:983822/1‑99 (MQ=21)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGTAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  1:349360/99‑1 (MQ=18)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  2:2457497/99‑1 (MQ=34)
                                       ggcggTGTGTGGCACCTGGAGGTGGGGTCATGACATGCTGGAGGGGGAGTGGGTCCCGGGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           <  2:3949763/99‑1 (MQ=2)
                                        gcggTGTTTGGCACCTCGATGTGGGCTCATTAAATCCTGCAGCTGAAGTGGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACCCGAGCTggg                                          >  1:1077733/1‑99 (MQ=11)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  1:3882896/1‑99 (MQ=21)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCAATTTAAAGTGGTACGCGAGCTggg                                          >  1:1208499/1‑99 (MQ=255)
                                         cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                         >  2:4266017/1‑99 (MQ=34)
                                         cggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                         <  2:506594/99‑1 (MQ=18)
                                          ggTGTTTGGCACCTCGATGTCGGCTGATCACATCCTGGAGCTGAAGGAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  1:4895312/99‑1 (MQ=16)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  1:4276622/99‑1 (MQ=33)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  2:3766240/99‑1 (MQ=33)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        >  1:2476180/1‑99 (MQ=18)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        >  1:773345/1‑99 (MQ=21)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                        <  2:3200931/99‑1 (MQ=18)
                                          gatgtTTTGCACCTCGATGACGGCTCATCACATCCTGAAGCTCAAGTACGTCGCAAGGCTATGGCTGGATCGCTTTTAAAGTCGCACGCGAACTAGGtc                                        <  2:891528/97‑2 (MQ=2)
                                            tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTa                                      <  2:1552160/99‑1 (MQ=34)
                                             gtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAg                                     <  2:5107285/99‑1 (MQ=33)
                                              tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                    >  2:3729098/1‑99 (MQ=34)
                                                tGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAc                                  <  1:1006203/99‑1 (MQ=33)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGTCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                               <  2:2848066/99‑1 (MQ=255)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                               <  2:4975525/99‑1 (MQ=21)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGATGAAGTAGGTCCAACAAGTCTGGCTGTTCGCCATTTAAAGAGGCCAGAGAGCTGGGTTAAGAGcgtc                               <  2:5180693/99‑1 (MQ=2)
                                                     ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                             <  2:5016877/99‑1 (MQ=33)
                                                     ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                             <  2:3136564/99‑1 (MQ=18)
                                                      cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            <  2:4682777/99‑1 (MQ=33)
                                                      cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            <  1:3548880/99‑1 (MQ=33)
                                                       tggATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCGCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                           <  2:3086554/97‑1 (MQ=255)
                                                       tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                           <  2:2729281/99‑1 (MQ=33)
                                                        cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgag                          >  1:4549484/1‑99 (MQ=33)
                                                         gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTATCACGTCGTgaga                         >  2:3555805/1‑99 (MQ=255)
                                                          aTGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTCGGTTTACAACGTCGTGAGAc                        <  2:3652133/99‑1 (MQ=16)
                                                             tCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGCCTGTTCGCCATTTAAAGTGGTACGCGCGCTGGGTTTAGAACGTCGTGGGACAGt                     >  1:4106701/1‑99 (MQ=11)
                                                              cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGtt                    >  1:663798/1‑99 (MQ=34)
                                                              cGGCTCAGCACATCCTGAAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGtt                    <  1:4852527/99‑1 (MQ=16)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  1:577615/1‑99 (MQ=34)
                                                               ggCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   <  2:1579418/99‑1 (MQ=18)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  >  2:3505452/1‑99 (MQ=33)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  >  2:750378/1‑99 (MQ=34)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  <  2:1879816/99‑1 (MQ=21)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt                <  2:1500037/99‑1 (MQ=18)
                                                                  tcatcaCATCCTGAAGCTAAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGGGCGTCGTGAGACAGTTCGGt                <  2:3349963/99‑1 (MQ=11)
                                                                    atcaCATCCTGGAGCTGAAGTAGGTCTCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              <  1:5164853/99‑1 (MQ=18)
                                                                    atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              >  1:3831487/1‑99 (MQ=34)
                                                                     tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc             >  1:325957/1‑99 (MQ=34)
                                                                      cacaTCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCCTTTAAAGTGGTACGCGAGCTGGGTTTAGACCGTAGTGCGGCAATTCGGtaact            <  2:1281068/99‑5 (MQ=11)
                                                                            cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc      >  2:4847611/1‑99 (MQ=34)
                                                                            cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc      >  1:2204684/1‑99 (MQ=33)
                                                                            cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc      <  1:1860706/99‑1 (MQ=34)
                                                                             tGGAGCGGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGATGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGcc     >  1:2094730/1‑99 (MQ=16)
                                                                               gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGt   <  1:3170887/99‑1 (MQ=32)
                                                                                                  |                                                                                 
GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946148‑3946326

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