Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,670,627 T→G 5.1% intergenic (+50/‑45) rplS → / → PP_1466 50S ribosomal protein L19/acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,670,6270TG 5.1% 308.0 / 5.8 98intergenic (+50/‑45)rplS/PP_146650S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (43/50);  new base G (3/2);  total (46/52)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

GGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACC  >  NC_002947/1670481‑1670770
                                                                                                                                                  |                                                                                                                                               
ggCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCg                                                                                                                                               >  1:298999/1‑149 (MQ=255)
ggCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCg                                                                                                                                               >  2:219495/1‑149 (MQ=255)
 gCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTc                                                                                                                                                <  7:294188/147‑1 (MQ=255)
 gCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCgg                                                                                                                                              <  8:407699/149‑1 (MQ=255)
 gCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCg                                                                                                                                               <  3:348704/148‑1 (MQ=255)
  ccGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGATAAAAGCAGCCTTCggg                                                                                                                                             <  1:252515/149‑1 (MQ=255)
   cGTGAAGCGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGgc                                                                                                                                            >  3:174104/1‑149 (MQ=255)
    gTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGgct                                                                                                                                           >  1:312797/1‑149 (MQ=255)
            tcgtGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGGGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttttt                                                                                                                                   <  3:36608/149‑1 (MQ=255)
            tcgtGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTCTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGACTTCGGGCTGCtttttt                                                                                                                                   >  3:306452/1‑149 (MQ=255)
              gtggtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGc                                                                                                                                 >  2:205393/1‑149 (MQ=255)
              gtggtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGc                                                                                                                                 >  3:282056/1‑149 (MQ=255)
                 gtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGt                                                                                                                               <  1:110195/148‑1 (MQ=255)
                  tgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGttt                                                                                                                             >  5:389045/1‑149 (MQ=255)
                   gACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtttt                                                                                                                            <  7:106494/149‑1 (MQ=255)
                    aCGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtttt                                                                                                                            <  5:361332/148‑1 (MQ=255)
                      gTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtt                                                                                                                              >  4:350916/1‑144 (MQ=255)
                      gTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtt                                                                                                                              <  3:350916/144‑1 (MQ=255)
                      gTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                          >  5:144492/1‑148 (MQ=255)
                      gTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                          >  2:395420/1‑148 (MQ=255)
                      gTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGACTTCGTGCTGCTGTTTTGCGTTTTGa                                                                                                                          >  8:291002/1‑148 (MQ=255)
                        gCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGaacc                                                                                                                       >  3:269807/1‑147 (MQ=255)
                          gTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAccc                                                                                                                     >  2:1755/1‑149 (MQ=255)
                                ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                               >  6:430519/1‑149 (MQ=255)
                                    gCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCGTCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtgtg                                                                                                           <  4:219559/149‑1 (MQ=255)
                                     cTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGa                                                                                                          <  4:427276/149‑1 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGa                                                                                                          <  8:104396/146‑1 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGa                                                                                                          >  7:104396/1‑146 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                        <  3:147431/148‑1 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                        <  5:78815/148‑1 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                        <  5:430519/148‑1 (MQ=255)
                                        tactacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAt                                                                                                         <  4:282056/147‑1 (MQ=255)
                                          ctacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGcc                                                                                                     <  6:389045/149‑1 (MQ=255)
                                          ctacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGc                                                                                                      >  6:189556/1‑148 (MQ=255)
                                          ctacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTCGAAACCCCGAAAATGTGAtt                                                                                                        <  5:227406/146‑1 (MQ=255)
                                          ctacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTCGAAACCCCGAAAATGTGAtt                                                                                                        >  6:227406/1‑146 (MQ=255)
                                           tacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                        >  8:101524/1‑145 (MQ=255)
                                           tacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGAtt                                                                                                        <  7:101524/145‑1 (MQ=255)
                                             ccTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCAt                                                                                                   <  7:291002/148‑1 (MQ=255)
                                             ccTGCGCGACCTGTCCGGCAAAGCCGCACGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCAt                                                                                                   <  5:209477/148‑1 (MQ=255)
                                              cTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATga                                                                                                 >  2:262085/1‑149 (MQ=255)
                                               tGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATgac                                                                                                >  4:217790/1‑149 (MQ=255)
                                                gcggGACCTGTCCCGCTAAGCCGCTCGCATCAACGAAAAACTTTCCTGAGTGCTGTTTGCCCGGTGGCCCAGTCCGCCTAGCGAGAAAAAACCAGCCTTCTAGCTGCTTTTTTGCGTTTTGaaa                                                                                                                        <  7:66635/124‑1 (MQ=255)
                                                gcccgACCTCTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGACCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCGTCGCGCTGCTTTTTTGCGTTTTGaaa                                                                                                                        >  8:66635/4‑124 (MQ=255)
                                                  gcgACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtttt                                                                                                                            <  5:74934/118‑1 (MQ=255)
                                                  gcgACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtttt                                                                                                                            >  6:74934/1‑118 (MQ=255)
                                                      ccTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGAccc                                                                                           >  1:412107/1‑147 (MQ=255)
                                                       cTGTCCGGCAAAGCCGCTCGCATCAAGGAATAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGGGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGAccc                                                                                           >  1:37867/1‑146 (MQ=255)
                                                          tCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATgac                                                                                                <  7:349542/138‑1 (MQ=255)
                                                          tCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATgac                                                                                                >  8:349542/1‑138 (MQ=255)
                                                          tCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGAc                                                                                      >  6:107748/1‑148 (MQ=255)
                                                            cGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAg                                                                                   <  4:174104/149‑1 (MQ=255)
                                                                aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                               <  5:189556/149‑1 (MQ=255)
                                                                aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                               >  4:116439/1‑149 (MQ=255)
                                                                aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                               <  1:1755/149‑1 (MQ=255)
                                                                aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCAGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                               >  6:349954/1‑149 (MQ=255)
                                                                 aaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCttttttc                                                                                                                                  >  4:271084/1‑96 (MQ=255)
                                                                 aaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCttttttc                                                                                                                                  <  3:271084/97‑2 (MQ=255)
                                                                  aGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATc                                                                             <  7:302916/149‑1 (MQ=255)
                                                                  aGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATc                                                                             <  6:144492/149‑1 (MQ=255)
                                                                     cGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAg                                                                          >  4:68292/1‑149 (MQ=255)
                                                                       cTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAgcg                                                                        >  2:153842/1‑149 (MQ=255)
                                                                        tCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAgcgc                                                                       >  5:144703/1‑149 (MQ=255)
                                                                                  gAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAAc                                                             <  7:261572/149‑1 (MQ=255)
                                                                                      aaCTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTc                                                         <  5:417561/149‑1 (MQ=255)
                                                                                       aCTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCg                                                        <  2:412107/149‑1 (MQ=255)
                                                                                        cTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCgg                                                       <  1:219495/149‑1 (MQ=255)
                                                                                          gTCCTGAGTGCAGTTTGGCCGGTGGCCGAGGACGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCATCGCCGCACCGAACTGTCGGt                                                      <  8:355587/148‑1 (MQ=255)
                                                                                             cTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGAcc                                                  >  7:156340/1‑149 (MQ=255)
                                                                                                 gTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCcgcg                                              <  8:188817/149‑1 (MQ=255)
                                                                                                  tGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGt                                             >  6:357324/1‑149 (MQ=255)
                                                                                                    cAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTa                                           <  1:262085/149‑1 (MQ=255)
                                                                                                       tttGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTAcc                                         >  1:103447/1‑148 (MQ=255)
                                                                                                       tttGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTAcc                                         >  5:263675/1‑148 (MQ=255)
                                                                                                           cccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAgg                                     <  2:103447/148‑1 (MQ=255)
                                                                                                           cccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAgg                                     <  6:263675/148‑1 (MQ=255)
                                                                                                           cccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGc                                    <  1:395420/149‑1 (MQ=255)
                                                                                                            ccagTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                   <  5:349954/146‑1 (MQ=255)
                                                                                                            ccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                   >  8:399108/1‑149 (MQ=255)
                                                                                                            ccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                   >  3:104448/1‑149 (MQ=255)
                                                                                                            ccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCg                                   <  7:399108/149‑1 (MQ=255)
                                                                                                              ggTGGCCGAGGCCGCCTAGCTAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGgt                                 <  5:107748/149‑1 (MQ=255)
                                                                                                                 ggCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGCACCCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGtt                              >  5:321606/1‑149 (MQ=255)
                                                                                                                   ccGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                            <  3:68292/149‑1 (MQ=255)
                                                                                                                    cGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTccc                           <  1:153842/149‑1 (MQ=255)
                                                                                                                       ggCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCaa                        <  5:43189/149‑1 (MQ=255)
                                                                                                                                gCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccacc                 <  3:217790/147‑1 (MQ=255)
                                                                                                                                aacgcaaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaa               >  6:338031/6‑149 (MQ=255)
                                                                                                                                 cGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaac              >  3:322419/1‑149 (MQ=255)
                                                                                                                                 acgcaaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTc                             >  4:10492/5‑134 (MQ=39)
                                                                                                                                 acgcaaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTc                             <  3:10492/130‑1 (MQ=39)
                                                                                                                                   gcaaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacca            <  1:93895/147‑1 (MQ=255)
                                                                                                                                   gcaaaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacca            >  2:93895/3‑149 (MQ=255)
                                                                                                                                     aaaaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAATTCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacc          <  2:312797/149‑1 (MQ=255)
                                                                                                                                     aaaaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccac           >  7:102288/1‑148 (MQ=255)
                                                                                                                                      aaaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacca         <  8:102288/149‑1 (MQ=255)
                                                                                                                                         aaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCacaa      >  6:247462/1‑149 (MQ=255)
                                                                                                                                             agcCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACAcc  <  4:306452/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                               
GGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACACC  >  NC_002947/1670481‑1670770

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: