Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,2460GA27.3% 464.8 / 187.1 275noncoding (2543/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (89/111);  new base A (23/52);  total (112/163)
Fisher's exact test for biased strand distribution p-value = 3.97e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

Only 100 of 280 total aligned reads displayed.

GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGT  >  NC_000913/3946148‑3946325
                                                                                                  |                                                                               
gggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAg                                                                                 >  1:2448996/1‑99 (MQ=33)
 ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAgg                                                                                <  1:116990/99‑1 (MQ=34)
    tAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTa                                                                             <  2:1332245/99‑1 (MQ=33)
     aaCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAt                                                                            <  1:5277560/99‑1 (MQ=34)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                          <  1:656731/99‑1 (MQ=34)
       cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATgg                                                                          <  1:2476020/99‑1 (MQ=34)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGc                                                                         >  1:451491/1‑99 (MQ=33)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGc                                                                         >  2:3027995/1‑99 (MQ=34)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGc                                                                         <  2:3816080/99‑1 (MQ=21)
        aGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGc                                                                         <  2:167301/99‑1 (MQ=21)
         ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCt                                                                        <  2:1268847/99‑1 (MQ=21)
         ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGACGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCt                                                                        <  2:4448444/99‑1 (MQ=255)
         ggCTGATACCGCCCAAGAGTGCATATCGACAGCGGTGTTTGGCACCTCGAGGTCGGCTCATCGCATCCTGGAGCTGAAGTAGGTCCCAAAGGTACGGCt                                                                        <  2:1182179/99‑1 (MQ=11)
            tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGtt                                                                     >  1:3300803/1‑99 (MQ=33)
               tACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGc                                                                  <  1:4318120/99‑1 (MQ=34)
                aCCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGcc                                                                 <  1:4459655/99‑1 (MQ=34)
                  cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                               >  2:1703419/1‑99 (MQ=33)
                  cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                               <  2:3698198/99‑1 (MQ=34)
                   gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAtt                                                              <  2:2198412/99‑1 (MQ=34)
                     ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTa                                                            >  1:3839841/1‑99 (MQ=34)
                     ccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACAACCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTa                                                            >  1:413184/1‑99 (MQ=18)
                      cAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaa                                                           >  2:1029103/1‑99 (MQ=34)
                      cAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaa                                                           >  2:4712955/1‑99 (MQ=34)
                      cAAGAGTTCATATCGACGGCGGGGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTGTGGCTGTTCGCCATTTag                                                           >  2:194025/1‑98 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTTGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:562827/99‑1 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                          <  1:2697928/99‑1 (MQ=34)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                          <  2:3482467/99‑1 (MQ=34)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaa                                                          <  2:4638753/99‑1 (MQ=34)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGTCATTTaaa                                                          <  1:1088518/99‑1 (MQ=255)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:4058218/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:1892271/99‑1 (MQ=18)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:2955477/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:2089866/99‑1 (MQ=21)
                       aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGAAGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTaaa                                                          <  2:717548/99‑1 (MQ=255)
                        agagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTaaga                                                         >  2:5028998/1‑97 (MQ=21)
                        agagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                         <  1:4936624/99‑1 (MQ=34)
                        agagTTCATACCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                         <  2:918295/99‑1 (MQ=255)
                         gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGt                                                        <  2:5352502/99‑1 (MQ=34)
                          agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTg                                                       <  1:919550/99‑1 (MQ=34)
                           gTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGCgg                                                      <  2:3173393/99‑1 (MQ=21)
                            ttCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGt                                                     <  1:994454/99‑1 (MQ=34)
                             tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                    <  1:1727185/99‑1 (MQ=34)
                             tCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                    <  1:2949018/99‑1 (MQ=33)
                             tCATATCGACGGCGGTGTTTGGCACCTAGATGTCGGATCATCACATCCTGGAGCTGAAGGAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTa                                                    <  2:1573581/99‑1 (MQ=11)
                              cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTAGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAc                                                   <  1:5119279/99‑1 (MQ=18)
                               atatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                  >  1:3920031/1‑99 (MQ=34)
                                tatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCCAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgc                                                 <  2:5389362/99‑1 (MQ=21)
                                 atCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                >  2:2308337/1‑99 (MQ=33)
                                  tCGACGGCGGGGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGa                                               >  2:1450651/1‑99 (MQ=18)
                                    gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                             <  1:3792174/99‑1 (MQ=34)
                                    gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGc                                             >  1:3521606/1‑99 (MQ=34)
                                      cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                           >  2:3249502/1‑99 (MQ=33)
                                      cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTg                                           >  2:2892449/1‑99 (MQ=34)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                          <  1:1083951/99‑1 (MQ=34)
                                       ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                          <  1:2102920/99‑1 (MQ=33)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                         >  1:1985837/1‑99 (MQ=33)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                         >  1:5490780/1‑99 (MQ=21)
                                        gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                         >  1:5532796/1‑99 (MQ=21)
                                         cggTGTTTGGCACCTCGATGTCGGCGCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                        <  2:4516275/99‑1 (MQ=18)
                                          ggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                       <  2:261523/99‑1 (MQ=21)
                                          ggTGTTCGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGtt                                       <  2:3435638/99‑1 (MQ=255)
                                           gtgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                      <  1:4714371/99‑1 (MQ=34)
                                           gtgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGttt                                      >  1:2759584/1‑99 (MQ=34)
                                            tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCCGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTCAAGTGGTACGCGACCTGGGTTTa                                     >  2:3654191/1‑99 (MQ=255)
                                            ggtTTGGCAGCTCGATATCGACTCATCACATCCTTCAGCTAAAGTAGCTCCCAAGGGTATCGCTGCACGCCATTTAAAGTCGTACGCGAGCTGGGTTTa                                     <  2:1107114/98‑1 (MQ=11)
                                               ttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                  <  1:3727672/99‑1 (MQ=34)
                                                tGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAc                                 <  2:2778079/99‑1 (MQ=34)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                              <  1:2363223/99‑1 (MQ=34)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                              <  1:1489233/99‑1 (MQ=33)
                                                   cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                              <  2:3939510/99‑1 (MQ=21)
                                                    aCCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                             <  2:3394543/99‑1 (MQ=33)
                                                       tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                          >  2:4364415/1‑99 (MQ=34)
                                                       tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                          <  2:3755366/99‑1 (MQ=33)
                                                         gATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgaga                        <  1:5158220/99‑1 (MQ=34)
                                                             tCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGCCTGTTAGCCATTTAAAGTGGAACGCGGGCTGGGTTTAGCGCGTAGTgcgatggt                    >  1:1396530/1‑95 (MQ=11)
                                                              cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGtt                   <  1:1292980/99‑1 (MQ=34)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                 >  1:4977898/1‑99 (MQ=34)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                 <  2:3088992/99‑1 (MQ=33)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                 <  2:3136797/99‑1 (MQ=34)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                 <  1:5298855/99‑1 (MQ=34)
                                                                gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                 >  1:2071475/1‑99 (MQ=21)
                                                                 cTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCgg                >  1:735523/1‑99 (MQ=34)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt               <  1:2260169/99‑1 (MQ=33)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt               >  1:3668727/1‑99 (MQ=34)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt               >  1:777104/1‑99 (MQ=33)
                                                                  tcatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGt               >  2:2033226/1‑99 (MQ=34)
                                                                   catcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc              >  2:2652135/1‑99 (MQ=34)
                                                                   catcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc              >  1:322833/1‑99 (MQ=34)
                                                                   catcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc              <  1:1465697/99‑1 (MQ=33)
                                                                   catcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc              >  1:5554554/1‑99 (MQ=33)
                                                                    atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc             >  2:1125558/1‑99 (MQ=33)
                                                                     tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc            >  2:1831751/1‑99 (MQ=33)
                                                                      cacaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt           >  2:4829154/1‑99 (MQ=34)
                                                                       acaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTa          <  2:4764926/99‑1 (MQ=33)
                                                                        caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTAt         <  2:2368112/99‑1 (MQ=33)
                                                                          tCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGGCGTGAGACAGTTCGGTCCCTATCt       >  2:4569518/1‑99 (MQ=18)
                                                                             tGGAGCTGACGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAAGTCGTGAGACAGTTCGGTCCCTATCTGcc    >  2:2526935/1‑99 (MQ=16)
                                                                              ggAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCg   >  1:2323935/1‑99 (MQ=33)
                                                                              ggAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCg   <  1:3077610/99‑1 (MQ=33)
                                                                               gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGt  <  2:5163933/99‑1 (MQ=32)
                                                                                                  |                                                                               
GGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGT  >  NC_000913/3946148‑3946325

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: