Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,532,022 T→G 6.0% W621G (TGG→GGG)  topB → DNA topoisomerase III

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,532,0220TG 6.0% 195.8 / 3.1 67W621G (TGG→GGG) topBDNA topoisomerase III
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base G (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.41e-01

GGAGCAGGCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTAACA  >  NC_002947/4531874‑4532169
                                                                                                                                                     |                                                                                                                                                   
ggAGCAGGCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGACCGGGGCGGCCCCGCCCt                                                                                                                                                     <  7:36383/149‑1 (MQ=255)
 gAGCAGGCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTg                                                                                                                                                    <  3:175316/149‑1 (MQ=255)
    cAGGCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGa                                                                                                                                                  >  1:365242/1‑148 (MQ=255)
      ggCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGaaa                                                                                                                                                >  6:246705/1‑148 (MQ=255)
       gcgcTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGGCGGGGCGGCCCCGCCCTGGaaaaa                                                                                                                                              <  2:365242/149‑1 (MQ=255)
        cgcTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAg                                                                                                                                             <  5:351209/149‑1 (MQ=255)
                       aaaGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGc                                                                                                                              >  2:70108/1‑149 (MQ=255)
                         aGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCgg                                                                                                                            <  1:282501/149‑1 (MQ=255)
                             gcgAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGTTGGCggg                                                                                                                           >  4:163607/1‑146 (MQ=255)
                                 gaTGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGc                                                                                                                     <  5:246705/148‑1 (MQ=255)
                                     ccccTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGc                                                                                                                >  6:299215/1‑149 (MQ=255)
                                      cccTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCg                                                                                                               >  8:234377/1‑149 (MQ=255)
                                      cccTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCGGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGc                                                                                                                >  1:190745/1‑148 (MQ=255)
                                      cccTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCGGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGc                                                                                                                >  1:184819/1‑148 (MQ=255)
                                       ccTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGa                                                                                                              >  1:208476/1‑149 (MQ=255)
                                        cTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAATCGTCGCGGTGGCGGGAAGAGCAAAGCGAc                                                                                                             <  1:230939/149‑1 (MQ=255)
                                           gAAGAGTTCGTCGCCCGGCAGTCAGTGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGTGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGAc                                                                                                             <  1:97359/147‑1 (MQ=255)
                                             agagTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTgg                                                                                                          >  1:141874/1‑147 (MQ=255)
                                             agagTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGaa                                                                                                        <  1:70108/149‑1 (MQ=255)
                                              gagTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAg                                                                                                       >  8:49044/1‑149 (MQ=255)
                                                gTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACt                                                                                                            >  8:10097/1‑142 (MQ=255)
                                                gTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACt                                                                                                            <  7:10097/142‑1 (MQ=255)
                                                 ttCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGc                                                                                                      >  2:350461/1‑147 (MQ=255)
                                                 ttCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTAGCAGCGCCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTCGTAAAAGAAACGTCGCGGTGGCGGGAAGAGCACAGCGCCTCGAAGc                                                                                                      >  6:247913/1‑147 (MQ=255)
                                                    gtcgCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGTGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGc                                                                                                 >  6:398102/1‑149 (MQ=255)
                                                      cgCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGa                                                                                               >  1:115633/1‑149 (MQ=255)
                                                        cccGGCAGTCAGCGTGTATGGGCATGCTGGTCGAGCGTTCAAGCGGCATGCGCATTACAATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCGGTAAAAATAAACGTCGCGGTGGCGGTAACAGCAAAGCGACTGGAAGCAAGGCGAcc                                                                                             <  2:184819/149‑1 (MQ=255)
                                                        cccGGCAGTCAGCGTGGATGGGCAAGCTGGTCTCGGGTTGCAGCGGCAGGCGCATGACCATCAGCGGCCCAGCGGCCGGGGCGGCCCCGCCCGGTAAAAACAAACGTCGCGGTGGCGGGAATAGCAAAGCGACTGGAAGCAAGGCGAcc                                                                                             <  2:190745/149‑1 (MQ=255)
                                                         ccGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGAAAGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAACAACCGTCGCGGGGGCGGGAAGAGCAAAGCGACTGGAAGAAAGGCGAccc                                                                                            >  2:238576/1‑149 (MQ=255)
                                                            gCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCcaa                                                                                         <  7:49044/149‑1 (MQ=255)
                                                               gTCAGCGTGGATGGACAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCcaagca                                                                                      >  4:76768/1‑149 (MQ=255)
                                                               gTCAGCGTGGATGGACAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCcaagc                                                                                       >  4:82215/1‑148 (MQ=255)
                                                                tCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCcaagcaa                                                                                     <  3:177998/149‑1 (MQ=255)
                                                                  aGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCcaagcaag                                                                                    <  2:141874/148‑1 (MQ=255)
                                                                          ggggCAAGCTGGTCGAGCGGTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGGGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAg                                                                           <  2:367538/148‑1 (MQ=255)
                                                                             gCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGgcg                                                                        >  5:154572/1‑149 (MQ=255)
                                                                              cAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGgcgc                                                                       <  8:35213/149‑1 (MQ=255)
                                                                               aaGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGgcgc                                                                       >  7:189535/1‑148 (MQ=255)
                                                                                 gCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTcgcg                                                                                                                                   <  7:206942/86‑1 (MQ=255)
                                                                                 gCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTcgcg                                                                                                                                   >  8:206942/1‑86 (MQ=255)
                                                                                     gTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGaaaa                                                                  <  5:299215/147‑1 (MQ=255)
                                                                                         aGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCGAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGa                                                            <  3:387420/149‑1 (MQ=255)
                                                                                          gCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAg                                                           <  6:154572/149‑1 (MQ=255)
                                                                                                  gCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTg                                                   >  4:196933/1‑149 (MQ=255)
                                                                                                    ggCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTa                                                 <  1:350461/149‑1 (MQ=255)
                                                                                                     gCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAg                                                <  3:163607/149‑1 (MQ=255)
                                                                                                           gcATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCgg                                          <  7:234377/149‑1 (MQ=255)
                                                                                                             aTGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGt                                        <  2:208476/149‑1 (MQ=255)
                                                                                                             aTGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGt                                        >  3:202541/1‑149 (MQ=255)
                                                                                                              tGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGTAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGtt                                       <  2:151164/149‑1 (MQ=255)
                                                                                                               gACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGtt                                       >  8:192059/1‑148 (MQ=255)
                                                                                                                aCCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTc                                      <  4:202541/148‑1 (MQ=255)
                                                                                                                aCCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCa                                     >  6:133052/1‑149 (MQ=255)
                                                                                                                     cAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCcgc                                >  1:334343/1‑149 (MQ=255)
                                                                                                                        cGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTGAACCcgctcc                             <  2:324801/149‑4 (MQ=255)
                                                                                                                            ccAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAgg                         >  5:369266/1‑149 (MQ=255)
                                                                                                                                cGGCCGGGGCGGCCCCGCCCTGGAAAAAGACACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAg                     <  5:247913/149‑1 (MQ=255)
                                                                                                                                 ggCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGc                    <  4:144691/149‑1 (MQ=255)
                                                                                                                                   ccGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGa                                              >  3:365083/1‑121 (MQ=255)
                                                                                                                                   ccGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGa                                              <  4:365083/121‑1 (MQ=255)
                                                                                                                                       ggcggcCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGgc              <  7:192059/149‑1 (MQ=255)
                                                                                                                                              cccGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGATGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCCCCCGCGAAGAGGCCCGCCCCGGCAGcacc        >  8:373690/1‑148 (MQ=255)
                                                                                                                                                cGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCTGGCAGCACAt       >  3:317708/1‑148 (MQ=255)
                                                                                                                                                  cccTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTaa    <  6:369266/149‑1 (MQ=255)
                                                                                                                                                  cccTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTaa    <  8:189535/149‑1 (MQ=255)
                                                                                                                                                   ccTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTaa    >  3:381950/1‑148 (MQ=255)
                                                                                                                                                   ccTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTAAc   <  3:76768/149‑1 (MQ=255)
                                                                                                                                                   ccTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTAAc   <  3:82215/149‑1 (MQ=255)
                                                                                                                                                     tGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTAACa  >  3:215998/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGAGCAGGCGCTGGACATGGTGCAAAGCGGCGAGATGCCCCTGGAAGAGTTCGTCGCCCGGCAGTCAGCGTGGATGGGCAAGCTGGTCGAGCGTTGCAGCGGCATGCGCATGACCATCAGCGGGCCAGCGGCCGGGGCGGCCCCGCCCTGGAAAAAGAAACGTCGCGGTGGCGGGAAGAGCAAAGCGACTGGAAGCAAGGCGACCCCAAGCAAGCCACGGCAGGCGCGAAAAAAGCGAGCCCCCCTGTAGGAGCGGGTTCACCCGCGAAGAGGCCAGCCCAGGCAGCACATTAACA  >  NC_002947/4531874‑4532169

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