Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,673,540 T→G 6.1% intergenic (‑39/‑373) 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,5400TG 6.1% 258.4 / 7.6 82intergenic (‑39/‑373)mdtJ/tqsAmultidrug/spermidine efflux pump membrane subunit MdtJ/autoinducer 2 exporter
Reads supporting (aligned to +/- strand):  ref base T (41/36);  new base G (3/2);  total (44/38)
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.04e-01

CATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATA  >  NC_000913/1673391‑1673687
                                                                                                                                                     |                                                                                                                                                   
cATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCtt                                                                                                                                                     >  1:163937/1‑150 (MQ=255)
cATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCtt                                                                                                                                                     <  2:1468578/150‑1 (MQ=255)
   tAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAt                                                                                                                                                  <  1:563252/150‑1 (MQ=255)
    aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAtt                                                                                                                                                 >  2:762820/1‑150 (MQ=255)
    aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAtt                                                                                                                                                 >  2:937094/1‑150 (MQ=255)
    aaaaTAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAtt                                                                                                                                                 >  1:852012/1‑150 (MQ=255)
     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:236179/150‑1 (MQ=255)
                                      cccATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCcgcg                                                                                                               <  1:1094641/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)
                                                 gACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAt                                                                                                                                                  <  2:1926235/104‑1 (MQ=255)
                                                 gACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAAt                                                                                                                                                  >  1:1926235/1‑104 (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)
                                                                                                                                                     |                                                                                                                                                   
CATTAAAATAAAGCCGCCATTTCCCTCACTGACGCTCGCCCATTTCATTGACAGCGTACCGGTAATTTCTGTAGCAATAGCCAGACCTAATAAAATCCAATAAATATACATTGTCCTTCTCCTGCAAGAGAATTATTTTAATTTTCGCTTAATTCAGCGAAACCAAAGTAAAGTGGTCAAGCCCCGCGCATAGCGCAAGCTTTGGCAGAAGAAAGGACTAAAGCGCGTTGCGCCAGTGCTGCTCACCTACGAGCAAGATAGTGGATGAGGTACGAAGATGGGTAAATGTAGAAAATA  >  NC_000913/1673391‑1673687

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: