Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,882,278 Δ1 bp 8.2% coding (893/1416 nt) oprN → multidrug RND transporter outer membrane protein OprN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,882,2780C. 8.2% 349.6 / 25.6 110coding (893/1416 nt)oprNmultidrug RND transporter outer membrane protein OprN
Reads supporting (aligned to +/- strand):  ref base C (55/46);  new base . (3/6);  total (58/52)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.89e-01

Only 100 of 113 total aligned reads displayed.

GCTCGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCG  >  NC_002947/3882129‑3882427
                                                                                                                                                     |                                                                                                                                                      
gCTCGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCAc                                                                                                                                                        <  1:276355/150‑1 (MQ=255)
  tCGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCg                                                                                                                                                      >  1:582791/1‑150 (MQ=255)
  tCGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCg                                                                                                                                                      >  2:1451403/1‑150 (MQ=255)
   cGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCcgccgc                                                                                                                                                     <  2:163152/150‑1 (MQ=255)
     gCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCa                                                                                                                                                   >  1:1596412/1‑150 (MQ=255)
      ccAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCaa                                                                                                                                                  >  1:537582/1‑150 (MQ=255)
       cAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACg                                                                                                                                                <  2:159802/150‑1 (MQ=255)
           gcCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCAGCAATGTCGTTCGCCAGCTGGCGGCGGCCACCGCCAACgaga                                                                                                                                             >  2:1398386/1‑147 (MQ=255)
             ccTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGc                                                                                                                                           <  1:1046632/150‑1 (MQ=255)
             ccTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGc                                                                                                                                           <  2:1064482/150‑1 (MQ=255)
                 aCGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTgg                                                                                                                                      >  2:1707704/1‑150 (MQ=255)
                 aCGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTgg                                                                                                                                      >  2:2125084/1‑150 (MQ=255)
                   gcgcTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCGGTGGcc                                                                                                                                    <  1:1797753/150‑1 (MQ=255)
                     gcTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                  <  1:1062447/150‑1 (MQ=255)
                      cTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGTGGCCACGCCAACGTTGGCGGTGgccgcc                                                                                                                                 <  1:2217322/150‑1 (MQ=255)
                            gcataCCTGGCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCCGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                           >  2:1552885/5‑150 (MQ=255)
                              aGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCt                                                                                                                         <  1:1451403/150‑1 (MQ=255)
                                aCCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGt                                                                                                                       <  1:1027391/150‑1 (MQ=255)
                                  cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                  <  2:1784526/137‑1 (MQ=255)
                                  cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                  >  1:1785228/1‑137 (MQ=255)
                                  cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                  >  1:1784526/1‑137 (MQ=255)
                                  cTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCAc                                                                                                                                  <  2:1785228/137‑1 (MQ=255)
                                   tGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTcc                                                                                                                    >  2:828770/1‑150 (MQ=255)
                                       gCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcg                                                                                                                >  2:1363317/1‑150 (MQ=255)
                                        cccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCcgcgc                                                                                                               <  2:139042/150‑1 (MQ=255)
                                         ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTgg                                                                                                                                            >  2:730619/1‑121 (MQ=255)
                                         ccAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTgg                                                                                                                                            <  1:730619/121‑1 (MQ=255)
                                          cAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCggtg                                                                                                                                      >  2:2357711/1‑121 (MQ=255)
                                          cAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCggtg                                                                                                                                      <  1:2357711/124‑4 (MQ=255)
                                           aaGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTc                                                                                                            >  1:934877/1‑150 (MQ=255)
                                           aaGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTc                                                                                                            >  1:1236491/1‑150 (MQ=255)
                                            aGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCa                                                                                                           <  2:622734/150‑1 (MQ=255)
                                              gCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGc                                                                                                         <  1:1289077/150‑1 (MQ=255)
                                                 cTGCCGTCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGACGCCAGCGGGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGGAAAAGACGACCAGATCCCGCGCGTCAGCCTc                                                                                                      >  1:532539/1‑150 (MQ=255)
                                                       gCCATCGCCAAGGCGTTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCCCAGCGCCGAGCGCCACCTGGCGGCGGCCACCGCCAACGTTCGCGTGGCCACCGCCGACGTGTTCCCCCGCGTCAGCCTCAGTGGc                                                                                                <  2:663649/150‑1 (MQ=255)
                                                          aTCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCggtgg                                                                                                                                     <  2:2339777/109‑5 (MQ=255)
                                                          aTCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCggtgg                                                                                                                                     >  1:2339777/1‑105 (MQ=255)
                                                             gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTc                                                                                          <  2:2274734/150‑1 (MQ=255)
                                                             gCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTc                                                                                          <  2:537582/150‑1 (MQ=255)
                                                              ccAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCg                                                                                         >  1:2387832/1‑150 (MQ=255)
                                                                    ctcacccGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                           >  2:688607/6‑110 (MQ=255)
                                                                    ctcacccGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGAc                                                                                                                           <  1:688607/105‑1 (MQ=255)
                                                                           ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGt                                                                                                                                        >  2:754928/1‑90 (MQ=255)
                                                                           ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGt                                                                                                                                        <  1:754928/90‑1 (MQ=255)
                                                                           ggtcggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                            >  1:1217822/1‑150 (MQ=255)
                                                                               ggtGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGccgcccgc                                                                        >  1:1633318/1‑150 (MQ=255)
                                                                                   aCCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCt                                                                    >  1:1570703/1‑150 (MQ=255)
                                                                                    cccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTc                                                                   <  1:528869/150‑1 (MQ=255)
                                                                                     ccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                   >  1:795145/1‑117 (MQ=255)
                                                                                     ccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                   <  2:795145/117‑1 (MQ=255)
                                                                                     ccGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACGCCAACGTTGGCggtg                                                                                                                                      >  2:456761/1‑78 (MQ=255)
                                                                                       gggTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCa                                                                <  1:842081/150‑1 (MQ=255)
                                                                                        ggTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                            >  1:226020/1‑137 (MQ=255)
                                                                                        ggTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACcgccgc                                                                            <  2:226020/137‑1 (MQ=255)
                                                                                          tGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGAt                                                             >  1:178156/1‑150 (MQ=255)
                                                                                          tGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGAt                                                             >  1:177754/1‑150 (MQ=255)
                                                                                            aGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCg                                                           <  2:1633318/150‑1 (MQ=255)
                                                                                              ctgctgCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                   <  1:2046222/108‑1 (MQ=255)
                                                                                              ctgctgCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGt                                                                                                   >  2:2046222/1‑108 (MQ=255)
                                                                                              ctactgcACCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCAAGCGGCGGCCACCGCCAACGTTGGCGAGGCCACCGCCGAACTGTTCCCGCGCGTCAGCATCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGCTCgcc                                                         >  2:1700457/4‑148 (MQ=255)
                                                                                                  tgcGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCGAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCtt                                                     <  2:1294226/150‑1 (MQ=255)
                                                                                                    cgccgccgcCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcg                                                   <  2:1766847/150‑1 (MQ=255)
                                                                                                     gccgccgccCGGACATCCGCAGCGCCGAGCGACAGCTGGCGGCGCCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcgg                                                  >  1:2226599/1‑150 (MQ=255)
                                                                                                      ccgccgccCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTcggc                                                 <  2:1570703/150‑1 (MQ=255)
                                                                                                            cccGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGGTGGCGTGGCCACCGCCGACCTGTCCCGGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCac                                           <  1:163213/150‑2 (MQ=255)
                                                                                                            cccGGACATCCGCAGCGCCGAGCGCCAACTGGCGGCGGCCACCGCCAACGTTGGCGGGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCATATCGCCTCTTCGGCGGcccc                                           >  2:163213/1‑148 (MQ=255)
                                                                                                             ccGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAAc                                          >  2:1355518/1‑150 (MQ=255)
                                                                                                              cGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACg                                         >  2:595116/1‑150 (MQ=255)
                                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAc                                                                                  >  1:949270/1‑104 (MQ=255)
                                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTAc                                                                                  <  2:949270/104‑1 (MQ=255)
                                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                    >  2:1412958/1‑150 (MQ=255)
                                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                    >  2:490541/1‑150 (MQ=255)
                                                                                                                   aTCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTg                                     >  1:2035947/1‑149 (MQ=255)
                                                                                                                    tCCGCATCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                   <  1:646110/150‑1 (MQ=255)
                                                                                                                    tCCGCAGCGCCGGGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTggg                                   <  1:1700457/150‑1 (MQ=255)
                                                                                                                     ccGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTCCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGc                                  >  2:1918936/1‑150 (MQ=255)
                                                                                                                     ccGCAGCGCCGCCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                    >  2:1248954/1‑147 (MQ=255)
                                                                                                                     ccGCAGCGCCGCCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTgg                                    <  1:1248954/147‑1 (MQ=255)
                                                                                                                       gCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCAc                                >  1:258590/1‑150 (MQ=255)
                                                                                                                         aGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTCCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTg                              <  1:1918936/150‑1 (MQ=255)
                                                                                                                          gcgcCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGa                                                                                                                            >  2:2292883/1‑55 (MQ=255)
                                                                                                                          gcgcCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGa                                                                                                                            <  1:2292883/55‑1 (MQ=255)
                                                                                                                              cGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGc                                                                                       <  2:963358/88‑1 (MQ=255)
                                                                                                                              cGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGc                                                                                       >  1:963358/1‑88 (MQ=255)
                                                                                                                                aGCGGCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGa                       <  2:1946227/150‑1 (MQ=255)
                                                                                                                                       gCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTAc                <  1:2125084/150‑1 (MQ=255)
                                                                                                                                         tGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGc                                                         <  1:1820704/107‑1 (MQ=255)
                                                                                                                                         tGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGc                                                         >  2:1820704/1‑107 (MQ=255)
                                                                                                                                           gcggcggcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCACGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTgg            >  1:1120042/1‑150 (MQ=255)
                                                                                                                                                 gcCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGcc      <  2:901591/150‑1 (MQ=255)
                                                                                                                                                  ccACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTg                                     >  1:1514445/1‑118 (MQ=255)
                                                                                                                                                  ccACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCt     >  1:1586190/1‑150 (MQ=255)
                                                                                                                                                  ccACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCt     >  1:1810110/1‑150 (MQ=255)
                                                                                                                                                  ccACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCt     <  1:556005/150‑1 (MQ=255)
                                                                                                                                                  ccACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCt     >  1:453193/1‑150 (MQ=255)
                                                                                                                                                   cACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCtt    >  1:311319/1‑150 (MQ=255)
                                                                                                                                                    aCCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCgg                                                          >  2:1757592/1‑95 (MQ=255)
                                                                                                                                                    aCCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTc   <  1:1970832/150‑1 (MQ=255)
                                                                                                                                                    aCCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTc   >  2:86253/1‑150 (MQ=255)
                                                                                                                                                     ccGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCg  >  2:294079/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GCTCGGCCAGCGCCCTGACGCGCTCAGCGTAGACCTGTCGCCCAAGGCCCTGCCGGCCATCGCCAAGGCGCTGCCGGTCGGTGACCCGGGTGAGCTGCTGCGCCGCCGCCCGGACATCCGCAGCGCCGAGCGCCAGCTGGCGGCGGCCACCGCCAACGTTGGCGTGGCCACCGCCGACCTGTTCCCGCGCGTCAGCCTCAGTGGCTTCCTCGGCTTTACCGCCGCCCGCGGCTCGCAGATCGGCTCTTCGGCGGCCAACGCCTGGGCACTGGGCCCGAGCATTACCTGGGCGGCCTTCG  >  NC_002947/3882129‑3882427

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: