Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,670,626 T→C 6.2% intergenic (+49/‑46) 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,6260TC 6.2% 191.1 / 2.3 64intergenic (+49/‑46)rplS/PP_146650S ribosomal protein L19/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base C (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 = 8.45e-01

CTGGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACAC  >  NC_002947/1670479‑1670769
                                                                                                                                                   |                                                                                                                                               
cTGGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCtt                                                                                                                                                >  1:206967/1‑149 (MQ=255)
    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                                                                                                                                    <  7:313665/125‑1 (MQ=255)
                                    aaGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCtttt                                                                                                                                    >  8:313665/1‑125 (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                          >  1:328015/1‑149 (MQ=255)
                                                                                                                      cGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTccc                          >  2:296668/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                           <  5:294773/140‑1 (MQ=255)
                                                                                                                              ccTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTcc                           >  6:294773/1‑140 (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)
                                                                                                                                                   |                                                                                                                                               
CTGGCCGTGAAACGTCGTGGTGACGTGCGTAAAGCCAAGCTGTACTACCTGCGCGACCTGTCCGGCAAAGCCGCTCGCATCAAGGAAAAACTGTCCTGAGTGCAGTTTGCCCGGTGGCCGAGGCCGCCTAGCGAGAAAAAAGCAGCCTTCGGGCTGCTTTTTTGCGTTTTGAAACCCCGAAAATGTGATTGCCATGACGACCCGAGACCAGGAAATCCAGCGCCGCACCGAACTGTCGGTGACCCGCGTTACCAAGGCGGTGTTCCCCAATACCACCAACCACCACAACAC  >  NC_002947/1670479‑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: ‑