Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,545 C→A 6.4% intergenic (‑44/‑368) mdtJ ← / → tqsA multidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,673,5450CA 6.4% 261.6 / 5.6 78intergenic (‑44/‑368)mdtJ/tqsAmultidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base C (38/35);  new base A (3/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.66e-01

AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTC  >  NC_000913/1673396‑1673692
                                                                                                                                                     |                                                                                                                                                   
aaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCTGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTc                                                                                                                                                     <  1:1434638/150‑1 (MQ=255)
aaaTAAAGCCGCCATTTCACTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCGTTCTCCTGCAAGAGAATTATTTTAATTTTCGATTAATTc                                                                                                                                                     >  1:1744437/1‑150 (MQ=255)
 aaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCa                                                                                                                                                    >  1:1307799/1‑150 (MQ=255)
   tAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGc                                                                                                                                                  >  2:1863125/1‑150 (MQ=255)
   tAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGc                                                                                                                                                  <  2:582554/150‑1 (MQ=255)
     aaGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGa                                                                                                                                                <  2:1679696/150‑1 (MQ=255)
       gacgccATTTCCCTCACTGACGATCGCCGATTTCATTGACAGCGTACCGGTAATTTATGAAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGaaa                                                                                                                                              >  2:1082738/3‑150 (MQ=255)
        ccgccATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAAc                                                                                                                                             >  1:1552272/1‑150 (MQ=255)
              tttCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagt                                                                                                                                       <  1:430258/150‑1 (MQ=255)
              tttCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTGAATTAAGCGaaaattaaaa                                                                                                                                       <  2:1930720/150‑8 (MQ=255)
                tCCCTCACTGAAGCTCTCCCATTTCATTGACAGCTTAACGGTAATTTCTGTATCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCaaagtaa                                                                                                                                     <  1:1403309/150‑1 (MQ=255)
                    tCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTg                                                                                                                                 <  1:246072/150‑1 (MQ=255)
                     cACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTgg                                                                                                                                <  1:1346069/150‑1 (MQ=255)
                     cACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTgg                                                                                                                                <  2:1082631/150‑1 (MQ=255)
                       cTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTCCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTc                                                                                                                              <  2:258430/150‑1 (MQ=255)
                         gACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCaa                                                                                                                            >  2:1612406/1‑150 (MQ=255)
                          aCGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAg                                                                                                                           <  1:527339/150‑1 (MQ=255)
                               cGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcg                                                                                                                      >  2:634143/1‑150 (MQ=255)
                               cGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcg                                                                                                                      <  2:1290025/150‑1 (MQ=255)
                                gCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgc                                                                                                                     >  1:914963/1‑150 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                    <  1:1094641/150‑1 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                    <  1:236179/150‑1 (MQ=255)
                                 cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                                    >  1:1253485/1‑150 (MQ=255)
                                     tttCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATa                                                                                                                >  1:615719/1‑150 (MQ=255)
                                      ttCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAg                                                                                                               <  1:835039/150‑1 (MQ=255)
                                       tCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAgc                                                                                                              >  2:42603/1‑150 (MQ=255)
                                                  gTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTgg                                                                                                   <  2:939000/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATAGCCACACCAAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGc                                                                                                  <  1:1157752/150‑1 (MQ=255)
                                                   tACCGGTAATTTCTGTAGCAATACCCAGACCTATTAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAAACAAAGTAAAGTGGTCAAGCCCCTCGCATAGCGCAAGCTTTGGc                                                                                                  <  2:175383/150‑1 (MQ=255)
                                                    aCCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCa                                                                                                 <  2:1307799/150‑1 (MQ=255)
                                                     ccGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCag                                                                                                >  1:1097007/1‑150 (MQ=255)
                                                       ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAATGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaa                                                                                              >  2:1664873/1‑150 (MQ=255)
                                                       ggTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaa                                                                                              <  2:529789/150‑1 (MQ=255)
                                                         tAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCagaaga                                                                                            >  1:577699/1‑150 (MQ=255)
                                                             ttCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAgg                                                                                        <  2:577699/150‑1 (MQ=255)
                                                                tGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACt                                                                                     >  2:883766/1‑150 (MQ=255)
                                                                agTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTAGCTTAATTCAGAGAAACCAAAGTAAATTGGTAAAGCCCCGCGCATAGCGCAAGCTATGCCAGAAGAAAGGAca                                                                                     >  2:153903/2‑149 (MQ=255)
                                                                 gTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTa                                                                                    >  2:1308452/1‑150 (MQ=255)
                                                                 gTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTa                                                                                    >  1:944818/1‑150 (MQ=255)
                                                                  tAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTaa                                                                                   >  2:1996756/1‑150 (MQ=255)
                                                                    gCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAg                                                                                 <  1:814469/150‑1 (MQ=255)
                                                                          gCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGtt                                                                           >  1:311582/1‑150 (MQ=255)
                                                                               aCCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGcc                                                                      <  2:87736/150‑1 (MQ=255)
                                                                                 cTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAg                                                                    <  1:634143/150‑1 (MQ=255)
                                                                                  taataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGt                                                                   >  1:866737/1‑150 (MQ=255)
                                                                                    ataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtgc                                                                 >  2:1711404/1‑150 (MQ=255)
                                                                                    ataaAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtgc                                                                 <  2:163937/150‑1 (MQ=255)
                                                                                       aaaTCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGtgctgc                                                              >  2:1147097/1‑150 (MQ=255)
                                                                                            cAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCAcc                                                         <  2:1253485/150‑1 (MQ=255)
                                                                                            cAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCAcc                                                         <  2:1403256/150‑1 (MQ=255)
                                                                                                 aataTACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGa                                                    >  1:1097346/1‑150 (MQ=255)
                                                                                                   tataCATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGc                                                  >  1:572338/1‑150 (MQ=255)
                                                                                                     taCATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCaa                                                >  2:786811/1‑150 (MQ=255)
                                                                                                     taCATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCaa                                                <  1:42603/150‑1 (MQ=255)
                                                                                                       cATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGa                                              >  1:1007040/1‑150 (MQ=255)
                                                                                                            tCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTg                                         >  2:1015503/1‑150 (MQ=255)
                                                                                                              cTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGa                                       <  1:762820/150‑1 (MQ=255)
                                                                                                                tctcCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATg                                     >  1:1830963/1‑150 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                                >  1:1929465/1‑150 (MQ=255)
                                                                                                                     tGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTa                                >  2:694787/1‑150 (MQ=255)
                                                                                                                      gCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTAc                               >  1:1528895/1‑150 (MQ=255)
                                                                                                                       cAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACg                              <  2:1097346/150‑1 (MQ=255)
                                                                                                                                  ttttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTaa                   >  2:1764101/1‑150 (MQ=255)
                                                                                                                                   tttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTaaa                  <  1:1711404/150‑1 (MQ=255)
                                                                                                                                   tttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTaaa                  <  2:311582/150‑1 (MQ=255)
                                                                                                                                    ttAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAAt                 >  1:1251456/1‑150 (MQ=255)
                                                                                                                                     ttggtttCGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATg                >  2:804847/5‑150 (MQ=255)
                                                                                                                                     tAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCTTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATg                <  2:565502/150‑1 (MQ=255)
                                                                                                                                       aTTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTa              >  2:464642/1‑150 (MQ=255)
                                                                                                                                          ttCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaa           >  2:1402042/1‑150 (MQ=255)
                                                                                                                                           tCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaaa          <  2:615719/150‑1 (MQ=255)
                                                                                                                                            cGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaaaa         <  1:1433243/150‑1 (MQ=255)
                                                                                                                                            cGCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGaaaa         >  1:115962/1‑150 (MQ=255)
                                                                                                                                             gCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTg                                         >  2:1033354/1‑117 (MQ=255)
                                                                                                                                             gCTGAATTAAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTg                                         <  1:1033354/117‑1 (MQ=255)
                                                                                                                                              cTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATa       >  1:871486/1‑150 (MQ=255)
                                                                                                                                                 aaTTCAGCGAAACCATAGTAAAGTGTTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCATTTCTGCTCACCTACGAGCAATATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACg    <  1:1578335/150‑1 (MQ=255)
                                                                                                                                                 aaTTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACg    <  1:1147097/150‑1 (MQ=255)
                                                                                                                                                  aTTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGt   >  1:75416/1‑150 (MQ=255)
                                                                                                                                                   ttCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTAGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGAAGGGTAAATGAAGAAAATAACgac  >  1:1598630/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATAACGTC  >  NC_000913/1673396‑1673692

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: