Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,386,395 C→A 11.0% intergenic (+32/+55) 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,386,3950CA11.0% 582.7 / 49.9 191intergenic (+32/+55)oprF/cobAporin F/uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (85/85);  new base A (12/9);  total (97/94)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 191 total aligned reads displayed.

AGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2386247‑2386541
                                                                                                                                                    |                                                                                                                                                  
aGCCGTGTTGGCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCg                                                                                                                                                   >  1:5606830/1‑150 (MQ=255)
 gCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCgg                                                                                                                                                  <  2:544657/150‑1 (MQ=255)
  ccGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGc                                                                                                                                                 <  2:4987213/150‑1 (MQ=255)
  ccGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGc                                                                                                                                                 <  1:3316530/150‑1 (MQ=255)
   cGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc                                                                                                                                                <  2:5565281/150‑1 (MQ=255)
    gtgtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCg                                                                                                                                               <  1:1751278/150‑1 (MQ=255)
    gtgtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCg                                                                                                                                               <  1:305337/150‑1 (MQ=255)
     tgtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCgg                                                                                                                                              <  1:5309623/150‑1 (MQ=255)
      gtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCggg                                                                                                                                             <  2:828987/150‑1 (MQ=255)
         gACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGttt                                                                                                                                          <  2:2896619/150‑1 (MQ=255)
         gACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGttt                                                                                                                                          >  2:4649733/1‑150 (MQ=255)
          aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt                                                                                                                                         <  2:1026549/150‑1 (MQ=255)
          aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt                                                                                                                                         <  1:3545440/150‑1 (MQ=255)
          aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt                                                                                                                                         >  2:395972/1‑150 (MQ=255)
           cTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGACGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttttt                                                                                                                                        <  2:2858603/150‑1 (MQ=255)
            tCGGTAGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTc                                                                                                                                       >  2:2515876/1‑150 (MQ=255)
              ggTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCtt                                                                                                                                     >  1:4868315/1‑150 (MQ=255)
              ggTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGTTTTTcat                                                                                                                                     <  1:3893690/150‑3 (MQ=255)
                 tGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCAAGAAAAACCCGGCCGAGGCCGGGTTTTTcatgaa                                                                                                                                  <  2:3256790/150‑6 (MQ=255)
                   gCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt                                                                                                                                >  1:3863819/1‑150 (MQ=255)
                     tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTACTTGGGTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                              >  1:1633116/1‑150 (MQ=255)
                     tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                              <  2:1910388/150‑1 (MQ=255)
                        ggCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGtt                                                                                                                           >  1:1813499/1‑150 (MQ=255)
                           gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                        >  1:4926126/1‑150 (MQ=255)
                           gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                        >  1:2788959/1‑150 (MQ=255)
                           gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                        >  1:4162406/1‑150 (MQ=255)
                                 cGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGc                                                                                                                  <  2:1923307/150‑1 (MQ=255)
                                     cGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTg                                                                                                              <  1:3420508/150‑1 (MQ=255)
                                      ggTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGATTTTCTTTGCCTGGGTTTTAACCGGCCCTGa                                                                                                             <  2:2499699/150‑1 (MQ=255)
                                       gTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                            <  2:2652510/150‑1 (MQ=255)
                                          gCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTATCCGGCCCTGAAGGc                                                                                                         >  2:3126912/1‑150 (MQ=255)
                                            tGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAAGCGGCCCTGAAGGCCg                                                                                                       >  1:720239/1‑150 (MQ=255)
                                            tGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCg                                                                                                       >  2:388934/1‑150 (MQ=255)
                                                   gCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGc                                                                                                >  2:1433272/1‑150 (MQ=255)
                                                     cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCgg                                                                                              >  2:3822066/1‑150 (MQ=255)
                                                     cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGt                                                                                                                                            >  1:1408270/1‑104 (MQ=255)
                                                        cGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCa                                                                                           <  2:4134410/150‑1 (MQ=255)
                                                         gaagaaGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAATACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg                                                                                          >  2:3888670/1‑150 (MQ=255)
                                                         gaagaaGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAg                                                                                          >  1:3593728/1‑150 (MQ=255)
                                                              agGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgc                                                                                     <  2:407455/150‑1 (MQ=255)
                                                                gCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcga                                                                                   >  2:364512/1‑150 (MQ=255)
                                                                gCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGCGGCCGGGTTTTTCTTTGCCTGGGTTTTAAACGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcga                                                                                   >  1:4737772/1‑150 (MQ=255)
                                                                 ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgat                                                                                  <  1:4314177/150‑1 (MQ=255)
                                                                  cGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgatt                                                                                 <  2:2790502/150‑1 (MQ=255)
                                                                   gTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgattg                                                                                <  1:4093335/150‑1 (MQ=255)
                                                                       tGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgattgcgat                                                                            >  2:4733756/1‑150 (MQ=255)
                                                                        gTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGGTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCAGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATc                                                                           >  2:2253588/1‑150 (MQ=255)
                                                                         tCAACCGTCGCGTTGAAGCTCAGGTAGAATCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGTCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATcc                                                                          >  1:907168/1‑150 (MQ=255)
                                                                           aaCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGc                                                                        >  2:1315635/1‑150 (MQ=255)
                                                                              cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcccc                                                                     >  2:4533231/1‑150 (MQ=255)
                                                                                tCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCAt                                                                   >  2:259043/1‑150 (MQ=255)
                                                                                   cgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGc                                                                <  1:3815229/150‑1 (MQ=255)
                                                                                      tGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctg                                                             <  1:5329015/150‑1 (MQ=255)
                                                                                       gAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgc                                                            <  2:1773133/150‑1 (MQ=255)
                                                                                       gAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgc                                                            >  2:4454705/1‑150 (MQ=255)
                                                                                            tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTc                                                       >  1:1164255/1‑150 (MQ=255)
                                                                                            tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCGTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTc                                                       >  1:5146015/1‑150 (MQ=255)
                                                                                            tCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCCTAGCCTGCTGCTc                                                       >  1:796844/1‑150 (MQ=255)
                                                                                              aGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGACCGGGTTTTTCTTTGCCTgg                                                                                                                              <  2:3095484/77‑1 (MQ=255)
                                                                                                gTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGc                                                                                         >  2:4273457/1‑112 (MQ=255)
                                                                                                gTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATACGCCCCATAGCCTGCTGCTCCCCa                                                   >  1:2103670/1‑150 (MQ=255)
                                                                                                  aGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTACATGAAAAACCCGGCCTCGGCCCGGGTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCTATCCGCCCCATAGCCTGCTGCTCCCCAac                                                 >  2:678980/1‑150 (MQ=255)
                                                                                                   gAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGc                                                                                                >  1:3769367/1‑102 (MQ=255)
                                                                                                   gAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAacc                                                >  2:112484/1‑150 (MQ=255)
                                                                                                              cAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAAt                                     <  2:2567589/150‑1 (MQ=255)
                                                                                                              cAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAAt                                     <  2:1813499/150‑1 (MQ=255)
                                                                                                                     ttGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAgg                              >  1:298514/1‑150 (MQ=255)
                                                                                                                      tGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGa                             >  1:819589/1‑150 (MQ=255)
                                                                                                                        gTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGAtt                           <  1:3109299/150‑1 (MQ=255)
                                                                                                                         tttGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTg                          <  1:4644066/150‑1 (MQ=255)
                                                                                                                            gAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                       <  2:907168/150‑1 (MQ=255)
                                                                                                                             aaGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCgg                      <  2:1633116/150‑1 (MQ=255)
                                                                                                                               gCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCCACCACCTCCCCAATCACCAGGATTTCCGGGc                    >  1:411/1‑150 (MQ=255)
                                                                                                                               gCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGc                    <  1:4351720/150‑1 (MQ=255)
                                                                                                                                 ttCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTa                  <  2:4162406/150‑1 (MQ=255)
                                                                                                                                  tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 >  1:4722265/1‑150 (MQ=255)
                                                                                                                                  tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 <  1:2649899/150‑1 (MQ=255)
                                                                                                                                  tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 <  1:2275944/150‑1 (MQ=255)
                                                                                                                                  tCATGAAAAACCCGGCCTCGGCCGGGTTTGTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCACCCACCTCCCCAATCACCAGGATTGCCGGGCTAc                 <  2:5597950/150‑1 (MQ=255)
                                                                                                                                   cATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACg                >  1:4809600/1‑150 (MQ=255)
                                                                                                                                    aTGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcc                                                                       <  2:4040651/94‑1 (MQ=255)
                                                                                                                                    aTGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcc                                                                       >  1:4040651/1‑94 (MQ=255)
                                                                                                                                     tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGca              >  2:5240035/1‑150 (MQ=255)
                                                                                                                                     agAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGca              >  1:4319272/2‑150 (MQ=255)
                                                                                                                                      gAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCCGCCTGCGATTGCGATCCGCCCCATAGCCTGCTCCTCCCCAACCACCTCCCCACTCACCAGGATTGCCGGGCTCCGcca             >  2:3009282/1‑148 (MQ=255)
                                                                                                                                      gAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTGTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGCTCCGCCCCATAGCCAGCTGCTCCCCAACCCCCTCCCCAATCACCAGGATTGCCGGGCTACGcaa             >  2:4644018/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            >  2:5255035/1‑150 (MQ=255)
                                                                                                                                       aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            <  1:5297109/150‑1 (MQ=255)
                                                                                                                                       aaaaaCCCGGCCTCGGCCGGGATTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa            <  2:3584173/150‑1 (MQ=255)
                                                                                                                                        aaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaac           <  2:2110976/150‑1 (MQ=255)
                                                                                                                                         aaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacc          >  1:2359817/1‑150 (MQ=255)
                                                                                                                                          aaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacca         >  1:2763617/1‑150 (MQ=255)
                                                                                                                                           aCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaa        <  2:4843357/150‑1 (MQ=255)
                                                                                                                                           aCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaa        <  2:4843403/150‑1 (MQ=255)
                                                                                                                                           aCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaa        <  2:2694523/150‑1 (MQ=255)
                                                                                                                                            cccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaa       <  2:151239/150‑1 (MQ=255)
                                                                                                                                            cccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaa       >  1:2504275/1‑150 (MQ=255)
                                                                                                                                            cccGCCCTCGGCCGGCTATTACTTTGCCAGCGATTTAACCGCCCCTCAAGGCCGCTCCTGCGCACAGCCACCGATTGCGATCCCCCCCATAGCCAGCTGCTCCCCAACCACCTCCCCAATCACCAGCATAGCCGCGCTACGcaaaccaaa       <  2:5302586/150‑1 (MQ=255)
                                                                                                                                             ccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac      <  1:2152729/150‑1 (MQ=255)
                                                                                                                                              cGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACg     >  1:2875591/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
AGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCC  >  NC_002947/2386247‑2386541

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑