Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,155,616 G→C 5.4% Q118H (CAG→CAC mdtB → multidrug efflux pump RND permease subunit MdtB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,155,6160GC 5.4% 254.7 / 4.5 74Q118H (CAG→CACmdtBmultidrug efflux pump RND permease subunit MdtB
Reads supporting (aligned to +/- strand):  ref base G (31/39);  new base C (2/2);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

ACCGCGCCGCTAGAACGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGATATGGTGGA  >  NC_000913/2155479‑2155765
                                                                                                                                         |                                                                                                                                                     
aCCGCGCCGCTAGAACGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAAc                                                                                                                                           >  1:860691/1‑150 (MQ=255)
       cgcTAGAACGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGa                                                                                                                                    <  2:1116642/150‑1 (MQ=255)
         cTAGAACGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGAcc                                                                                                                                  <  1:449770/150‑1 (MQ=255)
             aaCGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACt                                                                                                                              <  1:528320/150‑1 (MQ=255)
                 ccAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACttgtt                                                                                                                          >  1:313268/1‑150 (MQ=255)
                   aGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGc                                                                                                                        <  1:519908/150‑1 (MQ=255)
                    gTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGcc                                                                                                                       >  1:509694/1‑150 (MQ=255)
                     ttCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCg                                                                                                                      <  2:314664/150‑1 (MQ=255)
                     ttCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCg                                                                                                                      <  1:1400836/150‑1 (MQ=255)
                       cGGGCAGATTTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGAGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAg                                                                                                                    >  1:1074750/1‑150 (MQ=255)
                       cGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAg                                                                                                                    >  1:1079796/1‑150 (MQ=255)
                        gggCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGc                                                                                                                   <  2:785096/150‑1 (MQ=255)
                        gggCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGc                                                                                                                   <  2:780074/150‑1 (MQ=255)
                           cAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGAt                                                                                                                >  1:234312/1‑150 (MQ=255)
                             gATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCt                                                                                                              <  2:1012479/150‑1 (MQ=255)
                             gATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCt                                                                                                              <  1:1162084/150‑1 (MQ=255)
                               tGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGc                                                                                                            >  1:245584/1‑150 (MQ=255)
                               tGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGc                                                                                                            <  1:1203847/150‑1 (MQ=255)
                                gTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGcc                                                                                                           <  2:288404/150‑1 (MQ=255)
                                  cTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTa                                                                                                         <  2:365249/150‑1 (MQ=255)
                                   tGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTaa                                                                                                        <  1:827546/150‑1 (MQ=255)
                                         gAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACccgcc                                                                                                  <  1:808808/150‑1 (MQ=255)
                                           aaCAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCt                                                                                                              <  1:176692/136‑1 (MQ=255)
                                           aaCAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCt                                                                                                              >  2:176692/1‑136 (MQ=255)
                                              aGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGttt                                                                                             <  1:953915/150‑1 (MQ=255)
                                                aTGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTAc                                                                                           <  2:313268/150‑1 (MQ=255)
                                                 tGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACttg                                                                                                                            <  1:613617/116‑1 (MQ=255)
                                                 tGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACttg                                                                                                                            >  2:613617/1‑116 (MQ=255)
                                                 tGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACa                                                                                          >  1:556596/1‑150 (MQ=255)
                                                  gtcgtcGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGc                                                                                                                   >  2:1357294/1‑124 (MQ=255)
                                                  gtcgtcGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGc                                                                                                                   <  1:1357294/124‑1 (MQ=255)
                                                  gtcgtcGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAg                                                                                         >  2:899407/1‑150 (MQ=255)
                                                   tcgtcgCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGc                                                                                        <  1:799769/150‑1 (MQ=255)
                                                     gtcgCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCACGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCaa                                                                                      >  1:928409/1‑150 (MQ=255)
                                                     gtcgCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCACGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCaa                                                                                      >  1:933171/1‑150 (MQ=255)
                                                      tcgCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCaaa                                                                                     >  2:366358/1‑150 (MQ=255)
                                                         cAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAACTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTg                                                                                  >  2:497419/1‑150 (MQ=255)
                                                               tCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCg                                                                            >  1:681090/1‑150 (MQ=255)
                                                                 cggcggTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGc                                                                          <  1:616828/150‑1 (MQ=255)
                                                                  ggcggTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGCTTTACAGCAAAGTGAACCCGGCa                                                                         >  2:729233/1‑150 (MQ=255)
                                                                   gcggTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACTCGGCAg                                                                        >  1:1185558/1‑150 (MQ=255)
                                                                        gCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccg                                                                   <  2:1219511/150‑1 (MQ=255)
                                                                         cGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCTGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccgc                                                                  <  1:860029/150‑1 (MQ=255)
                                                                          gTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccgcc                                                                 <  1:971625/150‑1 (MQ=255)
                                                                          gTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccgcc                                                                 <  2:1074750/150‑1 (MQ=255)
                                                                          gTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccgcc                                                                 <  2:1079796/150‑1 (MQ=255)
                                                                           gcAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATccgccg                                                                <  2:97136/149‑1 (MQ=255)
                                                                             aGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGAt                                                              >  2:1215847/1‑150 (MQ=255)
                                                                                 aTCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCTATCATg                                                          <  2:556596/150‑1 (MQ=255)
                                                                                 aTCACTTTGCAGTTCCAGCTAACATTACCGCGCGATGTCGCCGAGCAGGAAGTGCAGTCCGCGATTAACGCTGCGACGATCTTGTTTTCTTGCTATGTGGATAACCTTGATGTTTACAGAAAATTTAACTTGGCAGATACGACGATCATg                                                          <  2:1328720/150‑1 (MQ=255)
                                                                                  tCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGa                                                         <  2:860691/150‑1 (MQ=255)
                                                                                   cACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGAc                                                        >  1:882723/1‑150 (MQ=255)
                                                                                         gCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGc                                                  <  1:817693/150‑1 (MQ=255)
                                                                                         gCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGc                                                  >  2:388830/1‑150 (MQ=255)
                                                                                            gTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGt                                               >  2:359907/1‑150 (MQ=255)
                                                                                              tCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCa                                             <  2:374470/150‑1 (MQ=255)
                                                                                                cAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCAcc                                           <  1:912599/150‑1 (MQ=255)
                                                                                                 aGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCt                                          >  2:820354/1‑150 (MQ=255)
                                                                                                  gCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTc                                         >  2:1396988/1‑150 (MQ=255)
                                                                                                   cTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTc                                         >  1:351614/1‑149 (MQ=255)
                                                                                                   cTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTc                                         <  2:351614/149‑1 (MQ=255)
                                                                                                        aTTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGc                                   >  1:419348/1‑150 (MQ=255)
                                                                                                            ccGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATg                               <  2:343221/150‑1 (MQ=255)
                                                                                                                  gATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATg                         <  1:1065035/150‑1 (MQ=255)
                                                                                                                     gTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACg                      >  1:821328/1‑150 (MQ=255)
                                                                                                                           gAGCAGGAAGTGCACGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTg                <  2:928409/150‑1 (MQ=255)
                                                                                                                           gAGCAGGAAGTGCACGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTg                <  2:933171/150‑1 (MQ=255)
                                                                                                                             gCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGa              >  1:419649/1‑150 (MQ=255)
                                                                                                                                ggAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGa           >  1:617771/1‑150 (MQ=255)
                                                                                                                                ggAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGAGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGa           >  2:1341560/1‑150 (MQ=255)
                                                                                                                                 gAAGTGCAGGCCGCGATTAACTCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCAGACTGTTTAAAGCAAAGTGAACCCGGAAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGACGATGACGCAAGTGGAAGat          >  2:494901/1‑150 (MQ=255)
                                                                                                                                 gAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCACCCGCCCTGCCGATGCCGCAAGTGGAAGat          >  2:490056/1‑150 (MQ=255)
                                                                                                                                 gAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGat          <  2:882723/150‑1 (MQ=255)
                                                                                                                                       cAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGATAtggtg    <  1:209864/150‑1 (MQ=255)
                                                                                                                                         ggCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGATATGGTGGa  >  2:73750/1‑150 (MQ=255)
                                                                                                                                         ggCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGATATGGTGGa  >  2:814475/1‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
ACCGCGCCGCTAGAACGCCAGTTCGGGCAGATGTCTGGCCTGAAACAGATGTCGTCGCAAAGTTCCGGCGGTGCGTCAGTTATCACTTTGCAGTTCCAGCTAACATTACCGCTCGATGTCGCCGAGCAGGAAGTGCAGGCCGCGATTAACGCTGCGACCAACTTGTTGCCGAGCGATCTGCCTAACCCGCCGGTTTACAGCAAAGTGAACCCGGCAGATCCGCCGATCATGACGCTCGCCGTCACCTCAACCGCCATGCCGATGACGCAAGTGGAAGATATGGTGGA  >  NC_000913/2155479‑2155765

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