Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,670,628 C→A 6.2% intergenic (+51/‑44) 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,6280CA 6.2% 203.6 / 3.5 64intergenic (+51/‑44)rplS/PP_146650S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base C (32/28);  new base A (2/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

CCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACAC  >  NC_002947/1670483‑1670769
                                                                                                                                                 |                                                                                                                                             
ccGTGAAACGGCGGGGGGAAGGGCGTAAAGCCAAGCTGTACTACCCGCGCGACCTGTGCGGCAAAGCCGCTCGCATCAACGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCggg                                                                                                                                            <  5:71635/149‑1 (MQ=255)
  gTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGgct                                                                                                                                          >  3:47681/1‑149 (MQ=255)
  gTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGgct                                                                                                                                          <  2:81331/149‑1 (MQ=255)
        cgtcgtGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGGAGCCCGAAGGCTGCtttt                                                                                                                                    <  2:185630/149‑1 (MQ=255)
         gtcgtGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCCGAAGGCTGCttttt                                                                                                                                   <  6:119385/149‑1 (MQ=255)
          tcgtGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttttt                                                                                                                                  <  2:49495/149‑1 (MQ=255)
             tggtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCttttt                                                                                                                                   >  8:22491/1‑145 (MQ=255)
             tggtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCttttt                                                                                                                                   <  7:22491/145‑1 (MQ=255)
               gtgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGtt                                                                                                                             <  1:35314/149‑1 (MQ=255)
                tgACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGttt                                                                                                                            <  4:115166/149‑1 (MQ=255)
                       cGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAc                                                                                                                      <  1:164513/148‑1 (MQ=255)
                        gTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAccc                                                                                                                    >  2:39057/1‑149 (MQ=255)
                        gTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGATCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAccc                                                                                                                    >  1:78627/1‑149 (MQ=255)
                            aGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGa                                                                                                                 >  7:144327/1‑148 (MQ=255)
                              ccAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGaaaa                                                                                                              >  1:170626/1‑149 (MQ=255)
                                aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttt                                                                                                                                    >  8:313665/1‑125 (MQ=255)
                                aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttt                                                                                                                                    <  7:313665/125‑1 (MQ=255)
                                        ctacCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGACCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAAtggcgatag                                                                                                     >  3:211666/1‑141 (MQ=255)
                                          acCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCAt                                                                                                  >  4:47772/1‑149 (MQ=255)
                                             tGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATgac                                                                                               >  5:341643/1‑149 (MQ=255)
                                                        tCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGAc                                                                                     <  3:94642/148‑1 (MQ=255)
                                                          cGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCa                                                                                   <  5:73292/148‑1 (MQ=255)
                                                           ggCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAACAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAgg                                                                                 <  5:77243/149‑1 (MQ=255)
                                                              aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                              <  1:39057/149‑1 (MQ=255)
                                                              aaaGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGaaa                                                                              <  2:75250/149‑1 (MQ=255)
                                                                aGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATc                                                                            <  7:32789/149‑1 (MQ=255)
                                                                 gCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATcc                                                                           <  3:101571/149‑1 (MQ=255)
                                                                     cTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                         <  6:170511/99‑1 (MQ=255)
                                                                     cTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGa                                                                                                                         >  5:170511/1‑99 (MQ=255)
                                                                          aTCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCa                                                                  >  2:134287/1‑149 (MQ=255)
                                                                           tCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCTGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAcc                                                                                                                     <  1:301792/97‑1 (MQ=255)
                                                                           tCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAAcc                                                                                                                     >  2:301792/1‑97 (MQ=255)
                                                                               ggAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGaa                                                             <  2:206967/149‑1 (MQ=255)
                                                                                gAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAAc                                                            >  8:240273/1‑149 (MQ=255)
                                                                                  aaaaCTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTg                                                          >  4:205456/1‑149 (MQ=255)
                                                                                         tCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGt                                                     <  8:144327/147‑1 (MQ=255)
                                                                                             gAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCcg                                               >  4:216741/1‑149 (MQ=255)
                                                                                                tGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGt                                            >  5:301976/1‑149 (MQ=255)
                                                                                                tGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGt                                            >  8:100708/1‑149 (MQ=255)
                                                                                                         cccGGTGGCCGAGGCCGCCTCGCGAGAAAAAAGCAGCGTTCGGGCTGCGTTTTTGCGGTTTGCAACCCCGAAAAGGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTCCCAAGGc                                   >  5:108973/1‑149 (MQ=255)
                                                                                                         cccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGc                                   >  2:120532/1‑149 (MQ=255)
                                                                                                          ccGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGATCTGTCGGTGACCCGCGTTACCAAGGCg                                  >  7:323988/1‑149 (MQ=255)
                                                                                                            ggTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGgt                                >  3:340028/1‑149 (MQ=255)
                                                                                                             gTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGgtg                               <  3:47772/149‑1 (MQ=255)
                                                                                                                 ccGAGGCCGCCTAGCGAGAAAAAAGCAGCCTACGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                           >  1:225594/1‑149 (MQ=255)
                                                                                                                  cGAGGCCGCCTAGCGATAAAAAAGCAGCGTTCGGGCTGCTTTTTTGCGTTTTGAAACGCCGAAAATGTGATTGCCATGACGAGCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGa                                                   <  2:187030/124‑1 (MQ=255)
                                                                                                                  cGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTccc                          >  2:296668/1‑149 (MQ=255)
                                                                                                                  cGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTccc                          >  1:328015/1‑149 (MQ=255)
                                                                                                                   gAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcccc                         <  2:78627/149‑1 (MQ=255)
                                                                                                                   gAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTccc                          >  6:213310/1‑148 (MQ=255)
                                                                                                                    aGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCga                                                                                        >  4:187085/1‑85 (MQ=255)
                                                                                                                    aGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCga                                                                                        <  3:187085/85‑1 (MQ=255)
                                                                                                                      gccgccTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAAt                      <  7:240273/149‑1 (MQ=255)
                                                                                                                          ccTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                           >  6:294773/1‑140 (MQ=255)
                                                                                                                          ccTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                           <  5:294773/140‑1 (MQ=255)
                                                                                                                          ccTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATacca                  <  7:100708/149‑1 (MQ=255)
                                                                                                                            aaaacgcaaaaaaGCAGCCCGAAGGCTGCTTTTATGCGTTTTGAAACCCCGAAAATGTTAATGCCAAGACGAACCGAGACCAAGAAAAACCACGCCGCACCCAACAATAGGTGACACGCGTTACCAAAGCGGTGTTACCCAAAaacac                 >  6:268393/8‑148 (MQ=255)
                                                                                                                                gagaAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaac             >  6:55821/1‑148 (MQ=255)
                                                                                                                                 agaAAAAAGCAGCCTTCGGGCGGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTCCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaacca           <  7:332146/149‑1 (MQ=255)
                                                                                                                                    aaaaaGCAGCCCGAAGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacca        >  4:120550/1‑149 (MQ=255)
                                                                                                                                     aaaaGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATaccaccaaccacca        >  8:70339/1‑148 (MQ=255)
                                                                                                                                        agcagcCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCacaac    >  6:247985/1‑149 (MQ=255)
                                                                                                                                         gcagcCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCacaaca   <  4:47681/149‑1 (MQ=255)
                                                                                                                                          cagcCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCAAGCGCCGCACCGAACTGTCGGTGACCAGCGTTACCAAGGCGGTGTTACCCAATACCACCAACCACCACAacac  >  7:337186/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACAC  >  NC_002947/1670483‑1670769

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