Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,381,402 T→G 10.2% intergenic (+31/+56) oprF → / ← cobA porin F/uroporphyrinogen‑III C‑methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,381,4020TG10.2% 156.9 / 9.9 59intergenic (+31/+56)oprF/cobAporin F/uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (21/32);  new base G (4/2);  total (25/34)
Fisher's exact test for biased strand distribution p-value = 3.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.23e-01

CAGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2381254‑2381549
                                                                                                                                                    |                                                                                                                                                   
cAGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCt                                                                                                                                                     <  7:74622/149‑1 (MQ=255)
 aGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTc                                                                                                                                                    <  3:115746/149‑1 (MQ=255)
        ttGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg                                                                                                                                              <  2:117600/148‑1 (MQ=255)
        ttGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg                                                                                                                                              <  2:101363/148‑1 (MQ=255)
           aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                          <  5:100761/149‑1 (MQ=255)
            cTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt                                                                                                                                         <  4:335147/149‑1 (MQ=255)
               ggTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttttt                                                                                                                                        <  3:166262/147‑1 (MQ=255)
                    gCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGACGTGCAGTCAACCGTCGCGTTGAAGCTCACGTAGAAGCCCAGGCCAAGTAATAGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTATTCTTTGcc                                                                                                                                 >  4:284912/1‑149 (MQ=255)
                      tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTg                                                                                                                               >  7:212100/1‑149 (MQ=255)
                       aCGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                              >  1:30491/1‑149 (MQ=255)
                        cGGCGAAACCCGTCCGGTTTCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                              <  3:225939/148‑1 (MQ=255)
                        cGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                              <  5:191066/148‑1 (MQ=255)
                         ggCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGt                                                                                                                            <  2:109926/149‑1 (MQ=255)
                         ggCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                            <  8:340505/133‑1 (MQ=255)
                         ggCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                            >  7:340505/1‑133 (MQ=255)
                           cGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGttt                                                                                                                          >  3:199291/1‑149 (MQ=255)
                              aaCCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTGGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                        <  1:221429/148‑1 (MQ=255)
                                 ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg                                                                                                                    >  3:20190/1‑149 (MQ=255)
                                         ttGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                            <  4:199291/149‑1 (MQ=255)
                                                 aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc                                                                                                    >  6:302755/1‑149 (MQ=255)
                                                    gCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCt                                                                                                  >  2:249137/1‑148 (MQ=255)
                                                      cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg                                                                                               <  6:19440/149‑1 (MQ=255)
                                                      cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg                                                                                               <  8:50725/149‑1 (MQ=255)
                                                       aCCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCt                                                                                                                                      <  4:347676/109‑1 (MQ=255)
                                                       aCCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCt                                                                                                                                      >  3:347676/1‑109 (MQ=255)
                                                          gaataaGGCCGTGCTGTCAACCTTCGCGTTGAATCTCATGTAGAACCCCAGGCCAATTAATTGGTTTTAAGCTTCATGAAAAACCCGGCCTCTGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCa                                                                                           <  7:135438/149‑1 (MQ=255)
                                                                  ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcga                                                                                   <  5:302755/149‑1 (MQ=255)
                                                                    gTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTTGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgatt                                                                                 <  6:371164/149‑1 (MQ=255)
                                                                      gCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg                                                                                          <  3:145500/138‑1 (MQ=255)
                                                                      gCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg                                                                                          >  4:145500/1‑138 (MQ=255)
                                                                         gTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTg                                                                                                                                   <  6:128925/94‑1 (MQ=255)
                                                                         gTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTg                                                                                                                                   >  5:128925/1‑94 (MQ=255)
                                                                               cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCAGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGccc                                                                      <  2:30491/149‑1 (MQ=255)
                                                                                 tCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCa                                                                    <  8:212100/149‑1 (MQ=255)
                                                                                     gTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGc                                                                >  2:49079/1‑149 (MQ=255)
                                                                                      ttGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGcc                                                               >  2:83599/1‑149 (MQ=255)
                                                                                         aaGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctg                                                             <  4:20190/148‑1 (MQ=255)
                                                                                           gCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgctg                                                          >  6:72041/1‑149 (MQ=255)
                                                                                               aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                      >  6:296955/1‑149 (MQ=255)
                                                                                               aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCAGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                      <  1:254124/149‑1 (MQ=255)
                                                                                                 gTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGcc                                                               <  2:321818/138‑1 (MQ=255)
                                                                                                 gTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGcc                                                               >  1:321818/1‑138 (MQ=255)
                                                                                                     aaGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAac                                                 >  3:268803/1‑148 (MQ=255)
                                                                                                         ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCCGCCTCGGCCGGGTTTTTCTTTGCCTGCGTTTTAACCGGCCCAGAAGGCCGCTCCTGCGGTCAGCCTGCGATAGCGATCCGCCCCATAGCCTGCTGCACCCCAACCACCt                                            >  1:143825/1‑149 (MQ=255)
                                                                                                           aGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTcc                                          <  5:72041/149‑1 (MQ=255)
                                                                                                               cAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCAGCTCCTGAGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCa                                       >  3:367012/1‑148 (MQ=255)
                                                                                                                     aTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCa                                <  1:83599/149‑1 (MQ=255)
                                                                                                                              aaGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAaccac                                              >  4:192625/1‑126 (MQ=255)
                                                                                                                              aaGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAaccac                                              <  3:192625/126‑1 (MQ=255)
                                                                                                                               aGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCgg                      <  4:268803/149‑1 (MQ=255)
                                                                                                                                    cATGAAAAACCCGGCCTGGGCCGGGTTTTTGTTTGCCTGGGTTTTAACCGTCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCGATAGGCTGCTGCTCCCCAACCTCCTCCCCAATCACCAGGATTGCCGGGCTa                  <  3:284912/148‑1 (MQ=255)
                                                                                                                                    aaagAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGGTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 >  5:107526/4‑149 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGc               >  1:158431/1‑148 (MQ=255)
                                                                                                                                        aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa             <  3:139282/149‑1 (MQ=255)
                                                                                                                                         aaaaCCCGGCGTCGGGCGTGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa             <  7:41630/148‑1 (MQ=255)
                                                                                                                                          aaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaac           >  5:157317/1‑149 (MQ=255)
                                                                                                                                          aaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTccc                                                     >  2:72407/1‑107 (MQ=255)
                                                                                                                                          aaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTccc                                                     <  1:72407/107‑1 (MQ=255)
                                                                                                                                               cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac      >  1:298539/1‑149 (MQ=255)
                                                                                                                                                   cTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATACCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGccc  <  2:116366/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CAGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2381254‑2381549

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: