Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,213,020 T→A 6.5% M15K (ATG→AAG)  mexC → multidrug RND transporter membrane fusion protein MexC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,213,0200TA 6.5% 241.0 / 7.0 77M15K (ATG→AAG) mexCmultidrug RND transporter membrane fusion protein MexC
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (2/3);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGC  >  NC_002947/3212873‑3213167
                                                                                                                                                   |                                                                                                                                                   
cTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATgg                                                                                                                                                   <  2:1484691/150‑1 (MQ=255)
     ggggTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATc                                                                                                                                              <  2:629222/150‑1 (MQ=255)
     ggggTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATc                                                                                                                                              <  1:307360/150‑1 (MQ=255)
               agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATgg                                                                                                                                                   >  1:80554/1‑135 (MQ=255)
               agagTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATgg                                                                                                                                                   <  2:80554/135‑1 (MQ=255)
                 agTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATCGGCCAAGGCGATCGCCATGGccaat                                                                                                                                  <  2:1409660/150‑4 (MQ=255)
                   tCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGccaatgc                                                                                                                                <  2:416719/150‑6 (MQ=255)
                     ttttATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGAGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGtgtg                                                                                                                              >  1:1169782/1‑150 (MQ=255)
                     ntttATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTACGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGggtg                                                                                                                              >  2:1101841/2‑150 (MQ=255)
                        tATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGc                                                                                                                           <  1:1500828/150‑1 (MQ=255)
                         aTGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGcc                                                                                                                          >  2:523853/1‑150 (MQ=255)
                           gATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATgg                                                                                                                                                   >  1:122691/1‑123 (MQ=255)
                           gATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATgg                                                                                                                                                   <  2:122691/123‑1 (MQ=255)
                            aTTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgc                                                                                                                       >  1:715255/1‑150 (MQ=255)
                             ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATg                                                                                                                                                    >  1:893031/1‑120 (MQ=255)
                             ttttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATg                                                                                                                                                    <  2:893031/120‑1 (MQ=255)
                               ttAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGNTNGGCGCANTGGCCATGGCGANCGCCTTGGCCGGGTGTGAGCCCgcngc                                                                                                                    <  1:1101841/150‑1 (MQ=255)
                                tAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcagca                                                                                                                   >  2:1632851/1‑150 (MQ=255)
                                    gagaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                               >  1:1449926/1‑150 (MQ=255)
                                    gagaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCACATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagc                                                                                                               <  2:1541933/150‑1 (MQ=255)
                                     agaCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                              >  2:473675/1‑150 (MQ=255)
                                       aCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAgg                                                                                                            >  2:968471/1‑150 (MQ=255)
                                       aCAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAgg                                                                                                            >  1:160030/1‑150 (MQ=255)
                                        cAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGa                                                                                                           <  1:1273614/150‑1 (MQ=255)
                                         aaTATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGAc                                                                                                          >  2:175233/1‑150 (MQ=255)
                                           tatTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAGGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTGGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACgg                                                                                                        >  1:53903/1‑150 (MQ=255)
                                           tatTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACgg                                                                                                        >  1:661834/1‑150 (MQ=255)
                                              tGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGc                                                                                                     <  2:197384/150‑1 (MQ=255)
                                               gACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCa                                                                                                    <  2:715255/150‑1 (MQ=255)
                                                    aTGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCaa                                                                                               <  2:942431/150‑1 (MQ=255)
                                                          gTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAg                                                                                         >  1:1344560/1‑150 (MQ=255)
                                                           tttGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGc                                                                                        >  1:1057487/1‑150 (MQ=255)
                                                              gAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                              <  1:22494/125‑1 (MQ=255)
                                                              gAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCagcagagca                                                                                                              >  2:22494/1‑125 (MQ=255)
                                                               aaaaGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTAcc                                                                                    <  1:968471/150‑1 (MQ=255)
                                                                   gTAGCATCCGCCCCGTTTTCAATCAAGGTACGGGTGTTTGATGTGCAGTACATGTTCGATAAGGGTGGGCGCATGGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGt                                                                                <  1:772571/150‑1 (MQ=255)
                                                                    tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                      >  2:65114/1‑95 (MQ=255)
                                                                    tAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCAAGGCGATCGCCATGGc                                                                                                                                      <  1:65114/95‑1 (MQ=255)
                                                                            gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                       >  2:1057420/1‑150 (MQ=255)
                                                                            gCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGt                                                                       <  1:417513/150‑1 (MQ=255)
                                                                              cccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCAAAAACCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCAAATGCGAGCTACCCGGTGGAAATCGTGa                                                                     >  1:516028/1‑150 (MQ=255)
                                                                               ccGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTggc                                                                    >  1:277645/1‑148 (MQ=255)
                                                                                   tttCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTg                                                                <  2:845215/150‑1 (MQ=255)
                                                                                      caatcaaGGTTCGGGTGTGTGGTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATTGCCTTGGCCGGgtg                                                                                                                                >  1:1106052/1‑83 (MQ=255)
                                                                                      caatcaaGGTTCGGGTGTGTGGTGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATTGCCTTGGCCGGgtg                                                                                                                                <  2:1106052/83‑1 (MQ=255)
                                                                                       aatcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGgcgc                                                            >  2:563170/1‑150 (MQ=255)
                                                                                        atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcag                                                                                                                     >  2:565458/1‑92 (MQ=255)
                                                                                        atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcag                                                                                                                     <  1:565458/92‑1 (MQ=255)
                                                                                        atcaaGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGcc                                                           >  1:421250/1‑150 (MQ=255)
                                                                                            aGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAg                                                       <  1:473675/150‑1 (MQ=255)
                                                                                                tCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAgg                                                   <  2:1057487/150‑1 (MQ=255)
                                                                                                    gtgtgtGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                      >  1:776830/1‑95 (MQ=255)
                                                                                                    gtgtgtGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTg                                                                                                      <  2:776830/95‑1 (MQ=255)
                                                                                                        gtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcag                                                                                                                     <  1:380928/76‑1 (MQ=255)
                                                                                                        gtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCgcag                                                                                                                     >  2:380928/1‑76 (MQ=255)
                                                                                                        gtgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTg                                           >  1:9275/1‑150 (MQ=255)
                                                                                                         tgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGc                                                      <  1:845031/138‑1 (MQ=255)
                                                                                                         tgATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGc                                                      >  2:845031/1‑138 (MQ=255)
                                                                                                                 aaGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCa                                                     >  2:930361/1‑131 (MQ=255)
                                                                                                                 aaGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCa                                                     <  1:930361/131‑1 (MQ=255)
                                                                                                                  aGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                          >  1:1084104/1‑93 (MQ=255)
                                                                                                                  aGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                          <  2:1084104/93‑1 (MQ=255)
                                                                                                                      aGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgcc                             <  2:277645/150‑1 (MQ=255)
                                                                                                                       ggTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccg                            <  2:661834/150‑1 (MQ=255)
                                                                                                                             aTAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTgccgggccgg                      <  2:713023/150‑1 (MQ=255)
                                                                                                                                aGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTc                   <  2:160030/150‑1 (MQ=255)
                                                                                                                                 gggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAg                                                       >  1:1163316/1‑113 (MQ=255)
                                                                                                                                 gggTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAg                                                       <  2:1163316/113‑1 (MQ=255)
                                                                                                                                       gcgcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                          <  1:1355518/72‑1 (MQ=255)
                                                                                                                                       gcgcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGa                                                                                          >  2:1355518/1‑72 (MQ=255)
                                                                                                                                       gcgcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                     >  1:900767/1‑125 (MQ=255)
                                                                                                                                       gcgcATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCt                                     <  2:900767/125‑1 (MQ=255)
                                                                                                                                          cATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc         <  1:563170/150‑1 (MQ=255)
                                                                                                                                          acccggCCAAGGCGATCGCCATGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGACGGTTATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgc         >  1:282450/5‑150 (MQ=255)
                                                                                                                                           aTTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATgcg        >  2:1066305/1‑150 (MQ=255)
                                                                                                                                               gCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTc    >  2:119597/1‑150 (MQ=255)
                                                                                                                                                 cATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGc  >  1:40715/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGAAGGGGTCATTAAGAGTCTTTTATGATTTTAATGAGACAATATTGACTCATGTGGGTTTGAAAAGTAGCATCCGCCCCGTTTTCAATCAAGGTTCGGGTGTGTGATGGGCAAGACAGGTTCGATAAGGGTGGGCGCATTGGCCATGGCGATCGCCTTGGCCGGGTGTGAGCCCGCAGCAGAGCAGGACGGTGCAGCCAATGCGAGCTACCCGGTGGAAATCGTGACATTGGCGCCCGAGCAGGTCGCCTTGGCGGCTGAACTGCCGGGCCGGGTCGAGCCGATGCGCGTCGC  >  NC_002947/3212873‑3213167

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