Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,381,402 T→G 7.3% 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,4020TG 7.3% 338.3 / 9.7 115intergenic (+31/+56)oprF/cobAporin F/uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (46/60);  new base G (6/3);  total (52/63)
Fisher's exact test for biased strand distribution p-value = 2.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 115 total aligned reads displayed.

TTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2381262‑2381549
                                                                                                                                            |                                                                                                                                                      
ttGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGGCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCggg                                                                                                                                                <  6:28560/149‑1 (MQ=255)
  gACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGt                                                                                                                                               <  8:2157/148‑1 (MQ=255)
  gACTCGGTTGGCTACGGAGAAACCCGTCCGGTAGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTACAAGCCCAGGCCAAGTAAATTGTTAGAAGCTTCATGAAAAAACCGGACGCGGACCGGtt                                                                                                                                              >  7:318835/1‑149 (MQ=255)
   aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                             <  6:98706/149‑1 (MQ=255)
   aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                             >  6:129797/1‑149 (MQ=255)
   aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                             >  8:196723/1‑149 (MQ=255)
   aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                             <  5:129797/149‑1 (MQ=255)
     tCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttttt                                                                                                                                           >  2:73345/1‑149 (MQ=255)
     tCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGCCCGTACTGTCAACCGACCCGTTGAACCTCAGGAAGAACCCCAGGCCAAGTAATTGGTTAGAAGCTTCATTAAAAACCCGGCCTAGGCCGGGtttt                                                                                                                                            >  3:276901/1‑148 (MQ=255)
       ggTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTc                                                                                                                                          <  1:364163/148‑1 (MQ=255)
       ggTTGGCTACGGCGAAACCCGCCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttttt                                                                                                                                           <  1:137686/147‑1 (MQ=255)
        gTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCATGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGTTTTTcat                                                                                                                                        <  5:8811/149‑3 (MQ=255)
         ttGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCttt                                                                                                                                       >  7:86714/1‑149 (MQ=255)
         ttGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCttt                                                                                                                                       <  4:231456/149‑1 (MQ=255)
            gCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCg                                                                                                                                                  >  5:2308/1‑135 (MQ=255)
             cTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt                                                                                                                                   >  3:238603/1‑149 (MQ=255)
             cTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt                                                                                                                                   >  2:332763/1‑149 (MQ=255)
             cTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                               >  8:114661/1‑137 (MQ=255)
             cTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                               <  7:114661/137‑1 (MQ=255)
              tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTg                                                                                                                                  >  7:113022/1‑149 (MQ=255)
              tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTg                                                                                                                                  >  5:354094/1‑149 (MQ=255)
                cGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                                 <  5:261440/148‑1 (MQ=255)
                   cGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGttt                                                                                                                             >  6:181776/1‑149 (MQ=255)
                     aaaCCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGtttt                                                                                                                            >  3:18898/1‑148 (MQ=255)
                     aaaCCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGGTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGtttt                                                                                                                            >  8:200551/1‑148 (MQ=255)
                      aaaccGTCCGGTTGCTGACATCGCCACCGATGAATGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCGGGGTTTTTCTTTGCCTGGGTTTTaa                                                                                                                          <  8:318835/146‑1 (MQ=255)
                         ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg                                                                                                                       <  3:24012/149‑1 (MQ=255)
                         ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCg                                                                                                                       >  7:335564/1‑149 (MQ=255)
                            tCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGcc                                                                                                                    <  8:335564/149‑1 (MQ=255)
                                gTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTg                                                                                                                 >  1:60291/1‑148 (MQ=255)
                                gTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGa                                                                                                               <  7:77240/149‑1 (MQ=255)
                                 ttGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                              <  8:113022/149‑1 (MQ=255)
                                     tGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGcc                                                                                                          >  6:319700/1‑149 (MQ=255)
                                        aaaCGCCACCGAAGAAGGCCGTGCTGTCAACCTTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCTTCAGGGCCggt                                                                                                       <  4:87271/148‑3 (MQ=255)
                                         aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg                                                                                                                                                 <  5:158290/107‑1 (MQ=255)
                                         aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg                                                                                                                                                 >  6:158290/1‑107 (MQ=255)
                                         aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGtt                                                                                                                              <  2:11633/126‑1 (MQ=255)
                                         aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGtt                                                                                                                              >  1:11633/1‑126 (MQ=255)
                                         aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc                                                                                                      <  2:268160/149‑1 (MQ=255)
                                           cGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCt                                                                                                    >  2:350561/1‑149 (MQ=255)
                                            gCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                              >  5:72462/1‑138 (MQ=255)
                                            gCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                              <  6:72462/138‑1 (MQ=255)
                                              cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg                                                                                                 <  1:38113/149‑1 (MQ=255)
                                               aCCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCg                                                                                                         >  1:216008/1‑140 (MQ=255)
                                               aCCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCg                                                                                                         <  2:216008/140‑1 (MQ=255)
                                               aCCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg‑‑g                                                                                              <  2:60291/149‑1 (MQ=255)
                                                 cGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCgg                                                                                                                      <  6:283915/126‑1 (MQ=255)
                                                  gaagaaGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCa                                                                                           <  4:238603/149‑1 (MQ=255)
                                                          ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcga                                                                                   <  5:183379/149‑1 (MQ=255)
                                                          ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcga                                                                                   >  2:106240/1‑149 (MQ=255)
                                                             tGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgattg                                                                                <  1:194870/149‑1 (MQ=255)
                                                             tGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgattg                                                                                <  5:369903/149‑1 (MQ=255)
                                                             tGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGCCCGGGTTTTTCTTTGCCGTGGTTTTAACAGGACCTGAAGGCAGCTACTGCG‑‑GTCAGCCtgcgaatg                                                                                >  8:142730/1‑149 (MQ=255)
                                                                tGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgattgcga                                                                             >  6:44240/1‑149 (MQ=255)
                                                                  tCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgattgcgat                                                                            <  5:319700/148‑1 (MQ=255)
                                                                    aaCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCg                                                                         >  6:185975/1‑149 (MQ=255)
                                                                       cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGccc                                                                      >  4:182114/1‑149 (MQ=255)
                                                                       cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGccc                                                                      <  1:332763/149‑1 (MQ=255)
                                                                       cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGccc                                                                      <  6:354094/149‑1 (MQ=255)
                                                                                    cTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCtgctgc                                                         >  6:241740/1‑149 (MQ=255)
                                                                                    cTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTAATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGTGTTTTAAACGGCCCTGAAGTCCGCTCCTCCG‑‑GTCAGCCTGCTATTGCGATCCGCACCAGAGCCtgctgc                                                         >  2:290634/1‑149 (MQ=255)
                                                                                       aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc                                                                                                      >  8:219394/1‑103 (MQ=255)
                                                                                       aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTc                                                                                                      <  7:219394/103‑1 (MQ=255)
                                                                                       aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgattgcg                                                                              <  8:197469/125‑1 (MQ=255)
                                                                                       aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                      <  1:106240/149‑1 (MQ=255)
                                                                                       aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                      <  5:241740/149‑1 (MQ=255)
                                                                                            gAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGTGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCaa                                                  <  7:142730/148‑1 (MQ=255)
                                                                                            gAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAac                                                 <  6:143902/149‑1 (MQ=255)
                                                                                            gAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAac                                                 <  8:310193/149‑1 (MQ=255)
                                                                                              aGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAacca                                               >  2:225946/1‑149 (MQ=255)
                                                                                                 ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGc                                                                                                                     <  5:318846/79‑1 (MQ=255)
                                                                                                 ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGc                                                                                                                     >  6:318846/1‑79 (MQ=255)
                                                                                                  cAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTc                                           <  8:80571/149‑1 (MQ=255)
                                                                                                     tctaaGTTATTGGTTTGATGCTTCATGAAAAACCCGGCCTCTGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCGCTGATGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTccg                                         <  2:295956/145‑2 (MQ=255)
                                                                                                        aaGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCaa                                      >  4:263894/1‑148 (MQ=255)
                                                                                                         aGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATc                                    <  1:350561/149‑1 (MQ=255)
                                                                                                         aGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCATTGCCTGGGTTTTAAACGGCACTGAAGGCCGCTACTGCG‑‑GTCAGGCAGCGATTGCGATACGCCCCATACACTGCTGAACCCCAACCACCTCACCacat                                    >  1:295956/1‑146 (MQ=255)
                                                                                                             aTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgatt                                                                                 <  8:264125/100‑1 (MQ=255)
                                                                                                             aTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtgcgatt                                                                                 >  7:264125/1‑100 (MQ=255)
                                                                                                                     gAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGcc                        >  5:242105/1‑149 (MQ=255)
                                                                                                                      aaGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                       >  7:210067/1‑149 (MQ=255)
                                                                                                                           tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGcccc                                                                     >  6:215436/1‑98 (MQ=255)
                                                                                                                           tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGcccc                                                                     <  5:215436/98‑1 (MQ=255)
                                                                                                                            aaagAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGTTAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGc                    >  1:98790/4‑147 (MQ=255)
                                                                                                                            aaagAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 >  7:89469/4‑149 (MQ=255)
                                                                                                                             aTGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 <  1:225946/148‑1 (MQ=255)
                                                                                                                              tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtg                                                                                      >  5:255396/1‑78 (MQ=255)
                                                                                                                              tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCtg                                                                                      <  6:255396/78‑1 (MQ=255)
                                                                                                                              tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                       <  8:17597/141‑1 (MQ=255)
                                                                                                                              tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                       >  7:17597/1‑141 (MQ=255)
                                                                                                                               gAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGAGTAAACATGCGTTTGCGAACCGCCCCATAGCATGCTGATCCCCAACCACATCCCAAATCACAAGGTTTGCCGGGATACg                >  1:111871/1‑149 (MQ=255)
                                                                                                                                aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGca              >  4:168015/1‑148 (MQ=255)
                                                                                                                                aaaaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGCTTGCCGGGCTACGca              >  6:256770/1‑148 (MQ=255)
                                                                                                                                 aaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            <  8:305478/149‑1 (MQ=255)
                                                                                                                                 aaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            >  7:124310/1‑149 (MQ=255)
                                                                                                                                  aaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            <  6:345006/148‑1 (MQ=255)
                                                                                                                                       cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac      <  3:263894/149‑1 (MQ=255)
                                                                                                                                       cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac      <  5:44240/149‑1 (MQ=255)
                                                                                                                                       cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacc          <  1:361354/145‑1 (MQ=255)
                                                                                                                                       cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCAGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCGCCCCAACCACCTCCCCAAACACCACGAATGCCCGGCTACGcaaaccaag       >  2:361354/1‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                                      
TTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCG‑‑GTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2381262‑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: