Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,946,246 G→A 28.1% noncoding (2543/2904 nt) rrlC → 23S ribosomal RNA of rrnC operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,946,2460GA28.1% 242.5 / 97.8 149noncoding (2543/2904 nt)rrlC23S ribosomal RNA of rrnC operon
Reads supporting (aligned to +/- strand):  ref base G (55/52);  new base A (12/30);  total (67/82)
Fisher's exact test for biased strand distribution p-value = 1.68e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

Only 100 of 157 total aligned reads displayed.

GGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946147‑3946326
                                                                                                   |                                                                                
ggggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAg                                                                                  <  1:1119320/100‑1 (MQ=32)
  ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAggg                                                                                <  1:2263555/100‑1 (MQ=33)
  ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAggg                                                                                <  1:1319074/100‑1 (MQ=33)
  ggATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAggg                                                                                <  1:721238/100‑1 (MQ=32)
    aTAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTa                                                                              >  1:609060/1‑100 (MQ=32)
     tAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAt                                                                             >  1:1935249/1‑100 (MQ=32)
      aaCAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATg                                                                            <  1:160579/100‑1 (MQ=32)
        cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGc                                                                          >  1:2842046/1‑100 (MQ=33)
        cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGc                                                                          <  2:1845640/100‑1 (MQ=33)
        cAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGc                                                                          <  1:1326753/100‑1 (MQ=18)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTg                                                                        <  1:2155055/100‑1 (MQ=18)
          ggCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTg                                                                        <  1:423535/100‑1 (MQ=17)
           gCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCATATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGt                                                                       <  1:1547264/100‑1 (MQ=255)
             tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTc                                                                     <  1:1212030/100‑1 (MQ=32)
             tGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTc                                                                     <  1:256756/100‑1 (MQ=32)
              gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                    <  2:1690979/100‑1 (MQ=33)
              gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                    <  2:2441735/100‑1 (MQ=33)
              gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                    >  2:2905135/1‑100 (MQ=32)
              gATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCg                                                                    <  1:2882182/100‑1 (MQ=33)
                tACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGcc                                                                  >  1:665901/1‑100 (MQ=32)
                tACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTAAGGCTGTTCg                                                                    >  1:270406/1‑98 (MQ=17)
                  ccGCCCAAGAGTTTATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCa                                                                 <  1:1983175/99‑1 (MQ=16)
                  ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                >  2:1594360/1‑100 (MQ=32)
                  ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                >  1:1457433/1‑100 (MQ=32)
                  ccGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAt                                                                >  1:305756/1‑100 (MQ=32)
                   cGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAtt                                                               <  1:1921803/100‑1 (MQ=32)
                    gCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCAttt                                                              >  2:735494/1‑100 (MQ=33)
                     cccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTa                                                             <  1:23676/100‑1 (MQ=33)
                     cccAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTa                                                             >  1:2123496/1‑100 (MQ=33)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                          <  2:930124/100‑1 (MQ=32)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:251094/100‑1 (MQ=33)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:2634221/100‑1 (MQ=17)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:505470/100‑1 (MQ=18)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:1026805/100‑1 (MQ=17)
                        aaGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAg                                                          <  1:1565238/100‑1 (MQ=18)
                         agagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGt                                                         <  2:261471/100‑1 (MQ=33)
                         agagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGGAGGTCCCAAATGTATGGCTGTTCGCCATTTAAAg                                                          <  1:2857370/99‑1 (MQ=255)
                         agagTTCATATCGACGGCGGTGTGTGGCACCTCGATGTCGGCTCATCACAGCCTAGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAga                                                         <  1:277177/100‑2 (MQ=11)
                          gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTg                                                        <  2:2430193/100‑1 (MQ=32)
                          gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTg                                                        <  1:1241490/100‑1 (MQ=17)
                          gagTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTg                                                        <  1:759385/100‑1 (MQ=18)
                          gagTTCATATCGACGGCGGTGTTTGGCACCCCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTg                                                        <  1:1756037/100‑1 (MQ=255)
                           agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTgg                                                       <  1:1077243/100‑1 (MQ=33)
                           agTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTgg                                                       <  1:1366550/100‑1 (MQ=33)
                            gTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGt                                                      >  2:2311639/1‑100 (MQ=33)
                               cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                   <  1:1231098/100‑1 (MQ=33)
                               cATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTAcg                                                   >  1:2899822/1‑100 (MQ=33)
                                 tatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTATAGTAGTAcgcg                                                 <  1:1525131/100‑1 (MQ=16)
                                 tatCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTAcgcg                                                 <  1:2513278/100‑1 (MQ=17)
                                     gACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCt                                             <  1:1680123/100‑1 (MQ=32)
                                       cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           >  2:972331/1‑100 (MQ=33)
                                       cggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTgg                                           >  2:2284032/1‑100 (MQ=18)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          <  1:2166990/100‑1 (MQ=33)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          <  1:236263/100‑1 (MQ=17)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  2:901503/1‑100 (MQ=18)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  2:2757231/1‑100 (MQ=18)
                                        ggcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTggg                                          >  2:2454127/1‑100 (MQ=17)
                                         gcggTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGt                                         >  2:423550/1‑100 (MQ=32)
                                             tgtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAg                                     <  1:1524780/100‑1 (MQ=33)
                                              gtTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGa                                    <  2:1398694/100‑1 (MQ=33)
                                               tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                   <  1:398920/100‑1 (MQ=32)
                                               tttGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGaa                                   <  2:1642659/100‑1 (MQ=32)
                                                   gCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                               >  2:1449109/1‑100 (MQ=33)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtc                               >  1:1822039/1‑99 (MQ=33)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGCATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                              >  2:510379/1‑100 (MQ=17)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                              <  1:568290/100‑1 (MQ=17)
                                                    cACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcg                              <  1:1012273/100‑1 (MQ=17)
                                                     aCCTCGATGTCGGCTCATCACATCATGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAcgtcgt                             <  2:2779989/100‑1 (MQ=18)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            <  1:991485/100‑1 (MQ=32)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            <  2:2804254/100‑1 (MQ=32)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            >  2:2534985/1‑100 (MQ=33)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAAAAcgtcgt                             >  2:2009732/1‑99 (MQ=17)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            >  2:1072064/1‑100 (MQ=17)
                                                      ccTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTg                            <  1:1950799/100‑1 (MQ=18)
                                                       cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                           >  2:2511935/1‑100 (MQ=33)
                                                       cTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTga                           >  2:2836799/1‑100 (MQ=33)
                                                        tCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgag                          >  2:270422/1‑100 (MQ=33)
                                                         cGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTgaga                         >  2:1129567/1‑100 (MQ=32)
                                                               cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   <  1:742448/100‑1 (MQ=33)
                                                               cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  2:1209972/1‑100 (MQ=33)
                                                               cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGATATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   <  2:2157273/100‑1 (MQ=18)
                                                               cGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTc                   >  2:561152/1‑100 (MQ=18)
                                                                ggCTCATCACATCCTGGAGCTGCAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCg                  >  1:1689394/1‑100 (MQ=18)
                                                                 gCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCgg                 >  2:2317219/1‑100 (MQ=33)
                                                                    catcaCATCCTGGAGCTGAAGTAGGTCCCTAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGATCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              <  1:1747617/100‑1 (MQ=16)
                                                                    catcaCATCCTGGAGCTGAAGTAGGTCCCAAAGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              <  1:909418/100‑1 (MQ=17)
                                                                    aatcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTcc              >  1:716761/2‑100 (MQ=18)
                                                                     atcaCATCCTGGAGCTGAAGTAGGTCCGAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTc               >  1:2601213/1‑98 (MQ=17)
                                                                     atcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTccc             <  2:716758/100‑1 (MQ=32)
                                                                      tcaCATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt            >  1:624073/1‑100 (MQ=33)
                                                                       cacaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTa           <  2:1551275/100‑1 (MQ=33)
                                                                        acaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCt            >  2:2604762/1‑98 (MQ=33)
                                                                        acaTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTAt          >  2:2001796/1‑100 (MQ=33)
                                                                        acaTCCTGGAGCTGAAGTAGGTCACAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAAGGTCGTGAGACAGTTCGGTCCCTAt          >  1:1903723/1‑100 (MQ=12)
                                                                         caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATc         >  2:1323760/1‑100 (MQ=33)
                                                                         caTCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATc         <  2:1001260/100‑1 (MQ=32)
                                                                            ccTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGc      >  2:309273/1‑100 (MQ=33)
                                                                             cTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGcc     >  2:1210632/1‑100 (MQ=32)
                                                                              tGGAGCTGAAGTAGGTCCCAAGGGTCTGGCTGTTCGCCATATAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGcc     >  1:1827223/1‑99 (MQ=12)
                                                                                gAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTg  <  1:2317208/100‑1 (MQ=31)
                                                                                                   |                                                                                
GGGGATAACAGGCTGATACCGCCCAAGAGTTCATATCGACGGCGGTGTTTGGCACCTCGATGTCGGCTCATCACATCCTGGAGCTGAAGTAGGTCCCAAGGGTATGGCTGTTCGCCATTTAAAGTGGTACGCGAGCTGGGTTTAGAACGTCGTGAGACAGTTCGGTCCCTATCTGCCGTG  >  NC_000913/3946147‑3946326

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: