Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,589,125 T→C 9.4% F511L (TTT→CTT)  acrD → aminoglycoside/multidrug efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,589,1250TC 9.4% 255.2 / 15.6 96F511L (TTT→CTT) acrDaminoglycoside/multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base T (40/47);  new base C (6/3);  total (46/50)
Fisher's exact test for biased strand distribution p-value = 3.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.44e-01

GCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAC  >  NC_000913/2588976‑2589273
                                                                                                                                                     |                                                                                                                                                    
gcACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCt                                                                                                                                                      <  1:13105/150‑1 (MQ=255)
       cGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGc                                                                                                                                               <  2:16596/150‑1 (MQ=255)
        ggTGCGATCTATCGCCAGTTCTCAATTACCATTGTTGCGGCGATGGAGCTGTCAGTACTGGTAGCGATGATGCTCACTCCGGCTCTGTGGGCCACACTACTTAAGCCACTGAAAAAAGGGGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                                                                              <  2:1427276/150‑1 (MQ=255)
        ggTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGcc                                                                                                                                              <  2:1119708/150‑1 (MQ=255)
          tGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTg                                                                                                                                            <  1:978100/150‑1 (MQ=255)
          tGCCATCTATAGCCAGTTTTCTATTAGCATTGTTTCGGCGATGGTGCTGTCAGTACTGTTAGCGATGATCCTCACTCCTGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGGGTGCATCATTGGCAAAAAGTCTTTTTTGCCTg                                                                                                                                            <  1:120233/150‑1 (MQ=255)
          ggCCATCTTTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGccca                                                                                                                                            <  2:1548500/149‑3 (MQ=255)
              aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                        <  2:1425432/150‑1 (MQ=255)
              aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                        <  1:1190754/150‑1 (MQ=255)
              aTCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGttt                                                                                                                                        <  2:27361/150‑1 (MQ=255)
               tCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGTCTGGTTTa                                                                                                                                       >  1:1250049/1‑150 (MQ=255)
                cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                                                                                                                                      <  2:835247/150‑1 (MQ=255)
                cTATCGCCAGTTCTCTATTACCATTGTTGCGGCGAGGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTaa                                                                                                                                      <  1:1320430/150‑1 (MQ=255)
                 tATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAc                                                                                                                                     <  1:1585937/150‑1 (MQ=255)
                  aTCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTACTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAAcc                                                                                                                                    >  1:996451/1‑150 (MQ=255)
                    cGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAg                                                                                                                                  <  2:1667337/150‑1 (MQ=255)
                     gCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGa                                                                                                                                 >  2:1606253/1‑150 (MQ=255)
                       catttCTCTATTTCCATTGTTGCGGCGTTGGGGCTGTCAGTACTTGTATAGATTATCCTCACTCCGGCGCTGTGGGGCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAATGCTTTTTTTCCTGGTTTAACCAGATg                                                                                                                               <  2:130866/147‑1 (MQ=255)
                         gTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                             <  2:537435/150‑1 (MQ=255)
                         gTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGtt                                                                                                                             >  1:962508/1‑150 (MQ=255)
                            ctctATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTaa                                                                                                                          >  2:1556495/1‑150 (MQ=255)
                                        tgttgCGGCGATGGTGCTGTCAGTACTGTTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGa                                                                                                              <  1:922979/150‑1 (MQ=255)
                                          ttgCGGCGATGGTGCTGTCAGTACTGGGAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAAc                                                                                                            <  2:30086/150‑1 (MQ=255)
                                             cggcgATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCt                                                                                                         <  1:294956/150‑1 (MQ=255)
                                                       gCTGTCAGTACTTGTCGCGATGATCCTCACTCCGGCTCTGTGGGCCACACTACTTAAGCCACTGAAAAAAGGGGAGCATCATCGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACAGCAACGCCGAACGCTACGAAAAAgg                                                                                               <  2:315668/150‑1 (MQ=255)
                                                       gCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgg                                                                                               <  1:1637527/150‑1 (MQ=255)
                                                        cTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAggg                                                                                              >  1:360190/1‑150 (MQ=255)
                                                         tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                                                                             <  2:420797/150‑1 (MQ=255)
                                                         tGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAgggg                                                                                             <  2:414544/150‑1 (MQ=255)
                                                          gTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGt                                                                                            >  2:34631/1‑150 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAGAAAAGCCTTTTTGcccat                                                                                                                                           >  2:630203/1‑98 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAGAAAAGCCTTTTTGcccat                                                                                                                                           <  1:630203/101‑4 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTACCCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                                                                          >  1:1059658/1‑150 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                                                                          <  2:800286/150‑1 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                                                                          >  1:885854/1‑150 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccat                                                                                                                                           >  2:282379/1‑98 (MQ=255)
                                                            cAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAAGCCTTTTTGcccat                                                                                                                                           <  1:282379/101‑4 (MQ=255)
                                                            cAGTACTGGGAGCGATGATCCTCTCTCCGGCTCTGTGTGCCACACTACTTAAGCCAGTGAAAAAAGGGGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTGTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTgg                                                                                          <  2:1279903/150‑1 (MQ=255)
                                                             aGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGc                                                                                         >  2:600359/1‑150 (MQ=255)
                                                              gTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCg                                                                                        >  2:171040/1‑150 (MQ=255)
                                                                 cTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGaaa                                                                                     >  2:1646720/1‑150 (MQ=255)
                                                                    gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                                                                  >  1:801648/1‑150 (MQ=255)
                                                                    gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                                                                  >  1:807677/1‑150 (MQ=255)
                                                                    gTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAtt                                                                                  >  1:808542/1‑150 (MQ=255)
                                                                            gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCg                                                                          >  2:1447773/1‑150 (MQ=255)
                                                                            gatCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCg                                                                          >  2:1576950/1‑150 (MQ=255)
                                                                                cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                                                                      <  2:363766/150‑1 (MQ=255)
                                                                                cTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACAGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGc                                                                      >  2:191965/1‑150 (MQ=255)
                                                                                  cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCCACGCCGAACGCTTCGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCt                                                                    >  1:89688/1‑150 (MQ=255)
                                                                                  cACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTGTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGAGAATTTTCGCCAccccttttt                                                                    <  1:1277283/150‑8 (MQ=255)
                                                                                      ccGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTgcgc                                                                >  2:155645/1‑150 (MQ=255)
                                                                                         gCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTgg                                                             >  1:332113/1‑150 (MQ=255)
                                                                                          ctctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGa                                                            <  1:1556495/150‑1 (MQ=255)
                                                                                           tctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAATAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAt                                                           <  1:1717468/150‑1 (MQ=255)
                                                                                           gctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAt                                                           <  2:221855/149‑1 (MQ=255)
                                                                                            ctGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGAtt                                                          >  2:486490/1‑150 (MQ=255)
                                                                                              gtgtgCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATtgt                                                        <  1:1676283/150‑1 (MQ=255)
                                                                                                   ccACACTTCTTAAGCCACTGAAAAAAGGTGATCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttt                                                   <  1:694779/150‑1 (MQ=255)
                                                                                                   ccACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAttt                                                   <  1:10032/150‑1 (MQ=255)
                                                                                                    cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTa                                                  >  1:151211/1‑150 (MQ=255)
                                                                                                    cacaCTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTa                                                  >  2:229418/1‑150 (MQ=255)
                                                                                                          actTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAAATCTCCACCGTAGCCTGCGCTGGATTGTCATTTATGTCct                                            >  2:685902/1‑150 (MQ=255)
                                                                                                             tAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctgct                                         >  1:1405097/1‑150 (MQ=255)
                                                                                                              aaGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCctgctg                                        >  1:1103048/1‑150 (MQ=255)
                                                                                                               aGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCGGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCtgctgc                                       <  2:1141410/150‑1 (MQ=255)
                                                                                                                gCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTgctgct                                      <  1:1624223/150‑1 (MQ=255)
                                                                                                                 ccACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCtt                                     >  2:1431708/1‑150 (MQ=255)
                                                                                                                    cTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggc                                  >  2:48703/1‑150 (MQ=255)
                                                                                                                     tGAAAAAAGGTGAGCATCCTGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACAGCAACGCCGAACTCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGCGATTTATGTCCTGCTGCTTggcg                                 <  1:298013/150‑1 (MQ=255)
                                                                                                                     tGAAAAAAGGTGAGCATCATGTGCAAAAAGGCTTGTTTGCCTGGTTTAACCAGATGTTTCACCGCAACGCCGAACGCTACGAAAACGGGGGGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcg                                 <  1:1491628/150‑1 (MQ=255)
                                                                                                                      gAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTggcgg                                >  2:598429/1‑150 (MQ=255)
                                                                                                                        aaaaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCa                              >  1:1310671/1‑150 (MQ=255)
                                                                                                                          aaaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATg                            >  1:1651810/1‑150 (MQ=255)
                                                                                                                           aaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTATCCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATgg                           <  1:1429877/150‑1 (MQ=255)
                                                                                                                           aaGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATgg                           <  1:437699/150‑1 (MQ=255)
                                                                                                                              gTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGgtgt                        >  2:1236796/1‑150 (MQ=255)
                                                                                                                                 aGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcc                     <  2:968839/150‑1 (MQ=255)
                                                                                                                                   catcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctg                   <  1:1606253/150‑1 (MQ=255)
                                                                                                                                    atcatGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgt                  >  1:44667/1‑150 (MQ=255)
                                                                                                                                    aaaccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgt                  >  1:295169/7‑150 (MQ=255)
                                                                                                                                      catGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttc                >  1:851225/1‑150 (MQ=255)
                                                                                                                                       ccaggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc               >  1:557599/4‑150 (MQ=255)
                                                                                                                                       atGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc               <  1:928082/150‑1 (MQ=255)
                                                                                                                                       atGGGCAAAAAGGCTTGTTTGCCGGGTTTAACCAGATGTTTAACCGCACCGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGtgttcctgttcc               <  1:894042/150‑1 (MQ=255)
                                                                                                                                        tGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct              >  2:1175219/1‑150 (MQ=255)
                                                                                                                                        caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGTTGCTTGGCGGCATGGTgttcctgttcct              >  1:577605/3‑150 (MQ=255)
                                                                                                                                        caggCAAAAAAGCCTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTgttcctgttcct              >  1:1110006/3‑150 (MQ=255)
                                                                                                                                          ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc            >  1:1512913/1‑150 (MQ=255)
                                                                                                                                          ggCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGc            >  2:1560969/1‑150 (MQ=255)
                                                                                                                                            cAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGt          >  2:1628423/1‑150 (MQ=255)
                                                                                                                                             aaacaGGCTTTTTTGCCTGGTTTAACCAGATGTTTTACCGCAACGCCGAAAGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGtt         <  2:874343/150‑1 (MQ=255)
                                                                                                                                             aaaaaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGtt         <  2:874341/150‑1 (MQ=255)
                                                                                                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGTCGAAAAGTCTCCACCGTAGCCTGCGCTGGAATGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc      <  2:830491/150‑1 (MQ=255)
                                                                                                                                                aaGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc      <  2:332113/150‑1 (MQ=255)
                                                                                                                                                aaGGCTTTGTTGCCTGGTTTAACCAGATGTTTATCCGCATCGGCGTACTCTACGTAAAAGGGGTGGCGTAAATTCTCCACGGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGc      <  1:398131/150‑1 (MQ=255)
                                                                                                                                                 aGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGAGTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGcc     <  1:347859/150‑1 (MQ=255)
                                                                                                                                                    cTTTATTGCCTGGTTTAACCAGATGTTTAACAGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAc  >  1:1608614/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCACCACCGGTGCCATCTATCGCCAGTTCTCTATTACCATTGTTGCGGCGATGGTGCTGTCAGTACTGGTAGCGATGATCCTCACTCCGGCTCTGTGTGCCACACTACTTAAGCCACTGAAAAAAGGTGAGCATCATGGGCAAAAAGGCTTTTTTGCCTGGTTTAACCAGATGTTTAACCGCAACGCCGAACGCTACGAAAAAGGGGTGGCGAAAATTCTCCACCGTAGCCTGCGCTGGATTGTGATTTATGTCCTGCTGCTTGGCGGCATGGTGTTCCTGTTCCTGCGTTTGCCGAC  >  NC_000913/2588976‑2589273

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