Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,383,368 C→A 7.3% 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,383,3680CA 7.3% 629.0 / 24.8 186intergenic (+32/+55)oprF/cobAporin F/uroporphyrinogen‑III C‑methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (95/77);  new base A (6/8);  total (101/85)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 187 total aligned reads displayed.

CAGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCCG  >  NC_002947/2383219‑2383515
                                                                                                                                                     |                                                                                                                                                   
    cGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGGAAAACCCGGCCTCGGcc                                                                                                                                                 <  2:22170/150‑1 (MQ=255)
    cGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc                                                                                                                                                 <  2:1796442/150‑1 (MQ=255)
    cGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGcc                                                                                                                                                 <  1:3982798/150‑1 (MQ=255)
       gtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCggg                                                                                                                                              <  2:2773079/150‑1 (MQ=255)
       gtTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCggg                                                                                                                                              <  2:127347/150‑1 (MQ=255)
          gACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGttt                                                                                                                                           >  2:611090/1‑150 (MQ=255)
          gACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGttt                                                                                                                                           <  1:2708017/150‑1 (MQ=255)
           aCTCGGTTGGCTCCGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGtttt                                                                                                                                          <  2:3348860/150‑1 (MQ=255)
           aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGtttt                                                                                                                                          <  1:2112010/150‑1 (MQ=255)
           aCTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGGACTGGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGtttt                                                                                                                                          <  2:3348816/150‑1 (MQ=255)
             tCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTc                                                                                                                                        >  1:894155/1‑150 (MQ=255)
              cGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGTTTTTca                                                                                                                                       <  1:1609432/150‑2 (MQ=255)
               ggTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCGAGGCCGGGTTTTTcat                                                                                                                                      <  1:4311518/150‑3 (MQ=255)
                gtgggCTCCGGCGAAACCCGTCCGGTTGCTGACAACGCCCCCGCAGAAGGCCGTGCTGTCACCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCACGTACTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCttt                                                                                                                                     <  2:3780290/147‑1 (MQ=255)
                 ttGTCTACGGAGAAACCCGTCCGGTTGCTGACAAAGCCACCGAAGAAGGCCGTGCTGTAAACCGTCTCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTg                                                                                                                                    <  1:4033858/150‑1 (MQ=255)
                    gCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCt                                                                                                                                 >  2:2453963/1‑150 (MQ=255)
                     cTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTAGAAGCTCAGGTAGACGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTg                                                                                                                                <  1:73267/150‑1 (MQ=255)
                      tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                               >  1:3334131/1‑150 (MQ=255)
                      tACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTgg                                                                                                                               >  2:3510168/1‑150 (MQ=255)
                        cGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGt                                                                                                                             <  1:1204648/150‑1 (MQ=255)
                          gCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGttt                                                                                                                           <  1:3510168/150‑1 (MQ=255)
                            gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                         <  1:2811924/150‑1 (MQ=255)
                            gAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTa                                                                                                                         >  2:2870299/1‑150 (MQ=255)
                               aCCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAAcc                                                                                                                      >  1:4101565/1‑150 (MQ=255)
                                 ccGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCgg                                                                                                                    <  2:2141545/150‑1 (MQ=255)
                                   gTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGcc                                                                                                                  >  2:2521379/1‑150 (MQ=255)
                                    tCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGccc                                                                                                                 <  2:750/150‑1 (MQ=255)
                                    tCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAAACGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGccc                                                                                                                 >  1:4427542/1‑150 (MQ=255)
                                     ccGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCt                                                                                                                <  1:3152295/150‑1 (MQ=255)
                                        gTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGaa                                                                                                             <  1:3433478/150‑1 (MQ=255)
                                           gCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGc                                                                                                          <  2:4118482/150‑1 (MQ=255)
                                                 aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCTTCAGGGCCggtta                                                                                                    <  1:2400074/150‑5 (MQ=255)
                                                 aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTcc                                                                                                    >  1:4216361/1‑150 (MQ=255)
                                                 aaCGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTcc                                                                                                    >  2:1201695/1‑150 (MQ=255)
                                                   cGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTg                                                                                                  >  1:2549865/1‑150 (MQ=255)
                                                    gCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGc                                                                                                 >  1:3231294/1‑150 (MQ=255)
                                                     ccACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGATTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCg                                                                                                <  1:1338563/150‑1 (MQ=255)
                                                      cACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCgg                                                                                               >  2:1840059/1‑150 (MQ=255)
                                                         cGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCa                                                                                            >  1:870994/1‑150 (MQ=255)
                                                           aagaagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGc                                                                                          <  2:1889164/150‑1 (MQ=255)
                                                            agaagGCCGTGCTGTCAACCGTCGCGTTTAAGCTCACGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGcc                                                                                         <  1:512191/150‑1 (MQ=255)
                                                            agaagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGcc                                                                                         <  2:2580447/150‑1 (MQ=255)
                                                             gaagGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCt                                                                                        >  2:1300115/1‑150 (MQ=255)
                                                                ggCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcg                                                                                     >  2:4506580/1‑150 (MQ=255)
                                                                  ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgat                                                                                   <  1:2028844/150‑1 (MQ=255)
                                                                  ccGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgat                                                                                   >  2:2842940/1‑150 (MQ=255)
                                                                       cTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCtgcgattgcga                                                                              <  2:2310246/150‑1 (MQ=255)
                                                                          tCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATcc                                                                           >  2:1683404/1‑150 (MQ=255)
                                                                            aaCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGc                                                                         >  2:4012707/1‑150 (MQ=255)
                                                                            aaCCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGc                                                                         >  2:3752023/1‑150 (MQ=255)
                                                                               cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcccc                                                                      >  1:2343870/1‑150 (MQ=255)
                                                                               cGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGcccc                                                                      >  1:3280163/1‑150 (MQ=255)
                                                                                gTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCGTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATAGCGATCCGCCCCa                                                                     >  1:3043283/1‑150 (MQ=255)
                                                                                  cgcgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATa                                                                   <  1:99011/150‑1 (MQ=255)
                                                                                   gcgTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTGGTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCCGCCTGCGATTGCGATCCGCCCCATAg                                                                  >  2:1137465/1‑150 (MQ=255)
                                                                                    cgTTGAAGCTCAGGTAGAAGCCCAGGCCATGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGc                                                                 <  1:4150384/150‑1 (MQ=255)
                                                                                         aaGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgct                                                            <  2:3098864/150‑1 (MQ=255)
                                                                                          aGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCctgctg                                                           <  1:1102575/150‑1 (MQ=255)
                                                                                              cAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                       >  2:1364940/1‑150 (MQ=255)
                                                                                              cAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTcc                                                       >  1:3160532/1‑150 (MQ=255)
                                                                                                   aGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAac                                                  <  2:2373592/150‑1 (MQ=255)
                                                                                                         ccAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTc                                            >  1:1170109/1‑150 (MQ=255)
                                                                                                          cAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTcc                                           >  1:855929/1‑150 (MQ=255)
                                                                                                          cAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTcc                                           >  1:3137320/1‑150 (MQ=255)
                                                                                                          cAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCCGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTcc                                           >  1:3864724/1‑150 (MQ=255)
                                                                                                           aGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTccc                                          >  2:3674579/1‑150 (MQ=255)
                                                                                                              ccAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCaa                                       <  1:2329113/150‑1 (MQ=255)
                                                                                                              ccAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCaa                                       <  1:609293/150‑1 (MQ=255)
                                                                                                                aaGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATc                                     <  1:647318/150‑1 (MQ=255)
                                                                                                                    aaTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCa                                 <  2:1326678/150‑1 (MQ=255)
                                                                                                                    aaTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCa                                 <  1:4065922/150‑1 (MQ=255)
                                                                                                                    aaTTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCa                                 <  1:2211307/150‑1 (MQ=255)
                                                                                                                      ttGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAgg                               >  1:4135808/1‑150 (MQ=255)
                                                                                                                       tGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGa                              >  1:87810/1‑150 (MQ=255)
                                                                                                                         gTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGAtt                            >  1:2438771/1‑150 (MQ=255)
                                                                                                                          tttGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTg                           >  1:442151/1‑150 (MQ=255)
                                                                                                                             gAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                        >  1:2950117/1‑150 (MQ=255)
                                                                                                                             gAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCg                        >  1:3097373/1‑150 (MQ=255)
                                                                                                                                   tCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTAc                  <  1:1554572/150‑1 (MQ=255)
                                                                                                                                      tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGca               >  1:2517023/1‑150 (MQ=255)
                                                                                                                                      tGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGca               >  1:2073933/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa              >  2:3367295/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa              >  2:2423035/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa              >  2:1261055/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCACCCCAATCACCAGGATTGCCGGGCTACGcaa              >  1:1937997/1‑150 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCTCGGCCGGGTTTTTCTTGGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATCGCCTGCTGCTCCCCAACCACCTCCCCACTCACCAGGAGTGCCGGGCGCCGcac              >  2:3367502/1‑149 (MQ=255)
                                                                                                                                       gAAAAACCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcac              >  2:2972452/1‑149 (MQ=255)
                                                                                                                                        aaaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaa             >  2:1801494/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaac            >  2:4201207/1‑150 (MQ=255)
                                                                                                                                         aaaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaac            <  2:894081/150‑1 (MQ=255)
                                                                                                                                          aaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacc           >  1:4182769/1‑150 (MQ=255)
                                                                                                                                          aaaCCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaa              >  2:3281754/1‑147 (MQ=255)
                                                                                                                                          aaaCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaacc           >  2:2173281/1‑150 (MQ=255)
                                                                                                                                            aCCCGGCCGAGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaa         >  1:981512/1‑150 (MQ=255)
                                                                                                                                              ccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac       <  1:3899391/150‑1 (MQ=255)
                                                                                                                                              ccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac       >  2:4081713/1‑150 (MQ=255)
                                                                                                                                              ccGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGcaaaccaaac       >  2:1807821/1‑150 (MQ=255)
                                                                                                                                                ggCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGc     <  1:617496/150‑1 (MQ=255)
                                                                                                                                                ggCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGc     <  1:617848/150‑1 (MQ=255)
                                                                                                                                                   cTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCcg  >  1:488423/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAGCCGTGTTGACTCGGTTGGCTACGGCGAAACCCGTCCGGTTGCTGACAACGCCACCGAAGAAGGCCGTGCTGTCAACCGTCGCGTTGAAGCTCAGGTAGAAGCCCAGGCCAAGTAATTGGTTTGAAGCTTCATGAAAAACCCGGCCTCGGCCGGGTTTTTCTTTGCCTGGGTTTTAACCGGCCCTGAAGGCCGCTCCTGCGGTCAGCCTGCGATTGCGATCCGCCCCATAGCCTGCTGCTCCCCAACCACCTCCCCAATCACCAGGATTGCCGGGCTACGCAAACCAAACGCCCG  >  NC_002947/2383219‑2383515

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: