Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,859 T→C 8.6% H29R (CAT→CGT) ‡ ampC ← penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8590TC 8.6% 302.9 / 17.5 117H29R (CAT→CGT) ‡ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base T (43/63);  new base C (6/4);  total (49/68)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-02

Only 100 of 122 total aligned reads displayed.

CTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTT  >  NC_000913/4378710‑4379008
                                                                                                                                                     |                                                                                                                                                     
cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                                                                                       <  2:159330/150‑1 (MQ=255)
cTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                                                                                       <  2:96837/150‑1 (MQ=255)
  aaCTCAAACCACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCGCTATAAGCGGGGTAATTGTGCGATGc                                                                                                                                                     <  1:925556/150‑1 (MQ=255)
    cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                                                                                   >  1:1183157/1‑150 (MQ=255)
      cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATCCCCTGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGAGATTCACaa                                                                                                                                                 >  1:955220/1‑150 (MQ=255)
         acaacGTATGCTGGGTGACGGGCTGCATTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGCGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAatat                                                                                                                                              <  1:240165/150‑1 (MQ=255)
           aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                                                                            >  2:478740/1‑150 (MQ=255)
            acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCATGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCAATAAGCGGGGTAATTGTGCATCGCACAattaccc                                                                                                                                           <  1:351917/150‑6 (MQ=255)
                ttGCTGGGTGCCGGGCTGCGTTTTGGCGATGTCCGCATAGCCCCATGTCAAGTAATAAGGTGTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                                                                                       <  2:666013/150‑1 (MQ=255)
                   cTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttg                                                                                                                                    <  2:1678129/150‑1 (MQ=255)
                     gtgtgACTGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                                                                                  <  2:1175918/150‑1 (MQ=255)
                      tgtgACGGGCTGCTTTTTGTAGATGTCCGCATAGCCCCATGTAAAGTAATAACGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 <  2:1275417/150‑1 (MQ=255)
                      tgtgACGGGCTGCTTTTTGGAGATGTCCGCATAGCACCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                                                                                 <  2:342265/150‑1 (MQ=255)
                       gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                                >  2:362161/1‑150 (MQ=255)
                       gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                                <  2:1095751/150‑1 (MQ=255)
                         taCGGGCTGATTTTTGGCGATGTCCGCATAGCCCAAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCATTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                              <  2:1280196/149‑1 (MQ=255)
                         gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCTGGGTAATTGTGCGATGCAAAATATCGTTGATTTGTTGAgg                                                                                                                              >  2:872995/1‑150 (MQ=255)
                            gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                                                                           <  1:809750/150‑1 (MQ=255)
                               ctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                                                                                        <  1:36777/150‑1 (MQ=255)
                                 gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                                                                                      >  1:1654857/1‑150 (MQ=255)
                                 gctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaa                                                                                                                      >  1:720706/1‑150 (MQ=255)
                                  ctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTTATAAGGTTTACCCTGATAAATTACCGCCAGCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCaaa                                                                                                                     <  2:1445161/150‑1 (MQ=255)
                                      ttGGCGATGTCCGCATATCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                                                                                 <  2:757976/150‑1 (MQ=255)
                                       tGGCGATGTCCGCATAGCCCCTTGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                                                                                <  1:384032/150‑1 (MQ=255)
                                       tGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                                                                                <  1:1247691/150‑1 (MQ=255)
                                          cGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGc                                                                                                             >  1:1120558/1‑150 (MQ=255)
                                           gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCa                                                                                                            >  1:1583807/1‑150 (MQ=255)
                                           gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                                                                            >  1:1568181/1‑150 (MQ=255)
                                           gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                                                                            >  2:512196/1‑150 (MQ=255)
                                                 cGCATAGCCCCATGTAAATTAATAAGGTTGACACTGATACATTACCGCCACCGCCAAACCCGGGATCATTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGAACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                                                                                      <  2:45563/150‑1 (MQ=255)
                                                 cGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                                                                                      >  2:772524/1‑150 (MQ=255)
                                                 cGCATAGCACCAGGTAAAGTAATAATGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTTTAAGCGTGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGc                                                                                                      <  1:880502/150‑1 (MQ=255)
                                                    aTAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                                                                                   <  2:68010/150‑1 (MQ=255)
                                                    aTAGCCCAAGGTAAAGTATTAATGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGAGCTTTTGTTGCTCTATAAGCGGGGTAATTGTTCGATTCACAATATCGTTGTTTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                                                                                   <  1:1469501/150‑1 (MQ=255)
                                                       gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTTTGCGATGCACAAAATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGAGGtaat                                                                                                >  2:130091/1‑150 (MQ=255)
                                                        ccccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCCCAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGAGGtaatt                                                                                               >  2:81763/1‑150 (MQ=255)
                                                        ccccAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGCGGTAATTGTGCGATGCACAATATCGATGATTTGTAGAGGGGCAGCAAATGTGGAGCAAGAGGCTGtaatt                                                                                               >  2:323106/1‑150 (MQ=255)
                                                         cccAGGTAAAGTAATAAGCTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaatta                                                                                              <  2:203290/150‑1 (MQ=255)
                                                           catgTAAAGTAATAAGGTGTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                                                                            <  1:1482693/147‑1 (MQ=255)
                                                           cAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                                                                            <  2:1412688/150‑1 (MQ=255)
                                                               tAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                                                                        <  1:516522/150‑1 (MQ=255)
                                                                aaaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGTGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                                                                                       >  1:223614/1‑150 (MQ=255)
                                                                aaaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                                                                                       >  2:1297772/1‑150 (MQ=255)
                                                                aaaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCCTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                                                                                       <  1:1391456/150‑1 (MQ=255)
                                                                     aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                  >  1:836895/1‑150 (MQ=255)
                                                                      ataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                                                                                 <  1:4591/150‑1 (MQ=255)
                                                                            tttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATAACGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt                                                                           >  2:697049/1‑150 (MQ=255)
                                                                             ttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGATGATATGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGTCGCAGAGCGTCGtt                                                                          >  2:343292/1‑150 (MQ=255)
                                                                              tACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGttt                                                                         <  2:1676739/150‑1 (MQ=255)
                                                                                cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                       <  1:1406207/150‑1 (MQ=255)
                                                                                cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                       <  1:1354827/150‑1 (MQ=255)
                                                                                  cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                     >  1:426629/1‑150 (MQ=255)
                                                                                       aaaTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAg                                                                <  1:590024/150‑1 (MQ=255)
                                                                                         aTTCCCGCCCCTGCCATAGCGGGGATCTTTTGTTGCTCTATTCGCGGGGTATTTGGGCGATTCAGAATATCGTTTTTTGGTTGAGGGGCAGCTAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                              <  1:172219/150‑1 (MQ=255)
                                                                                         aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                              >  1:1444901/1‑150 (MQ=255)
                                                                                         aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                              >  1:993622/1‑150 (MQ=255)
                                                                                         aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                              >  1:1182589/1‑150 (MQ=255)
                                                                                            accgccaccgccATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt                                                           >  2:1328885/1‑150 (MQ=255)
                                                                                              cgccaccgccaTACCCGGGAGCTTTTGTTGCTCTCTAAGCGGGGTAATTGTGCGAGGCAAAATATCGTTGATTTGTTGAGGGGCAGCAAATGGGGAGCAAGAGGCGGTAATTAATAAGGAGCAGAGCGTCGTTTTGAACATAGGGTCTgg                                                         <  2:450083/150‑1 (MQ=255)
                                                                                               gccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATTGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGt                                                        <  1:1696339/150‑1 (MQ=255)
                                                                                                 caccgccaTACCCGTGATCTTTTGTTGCTCAATAAGCGGGGTAATTGTGCGATTCACAATATCGTTGATTTGTTGAGTGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGttt                                                      <  2:976681/150‑1 (MQ=255)
                                                                                                 caccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGttt                                                      <  2:983488/150‑1 (MQ=255)
                                                                                                   ccgcccTACCCGGGAGCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCACGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                    <  2:294218/150‑1 (MQ=255)
                                                                                                     gccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                  <  1:478740/150‑1 (MQ=255)
                                                                                                     gccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                  <  2:234240/150‑1 (MQ=255)
                                                                                                         tACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCAATACaa                                              >  1:394359/1‑150 (MQ=255)
                                                                                                           cccGGGCTCTTTGGTTGCTCTATCAGCGGGGCATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAAc                                           <  1:1067244/150‑1 (MQ=255)
                                                                                                               ggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGc                                        <  2:1120558/150‑1 (MQ=255)
                                                                                                                gATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGcc                                       >  1:1676980/1‑150 (MQ=255)
                                                                                                                 aTCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGccc                                      <  1:1471973/150‑1 (MQ=255)
                                                                                                                    ttttGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGcc                                   >  2:841556/1‑150 (MQ=255)
                                                                                                                        gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                               <  2:577116/150‑1 (MQ=255)
                                                                                                                             tctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCggg                          <  2:1178712/150‑1 (MQ=255)
                                                                                                                              ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                         >  2:392785/1‑150 (MQ=255)
                                                                                                                               tataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                        <  2:1583807/150‑1 (MQ=255)
                                                                                                                               aattaGCGGGGTAATTGTGCGTTGCACAATATGGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                        <  2:1568181/146‑1 (MQ=255)
                                                                                                                                ataAGCGTGGTATTTGTGCGATGCACAATATCGTTGATTTGTTGTGCGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCCCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                       <  1:81763/150‑1 (MQ=255)
                                                                                                                                ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCg                       <  2:507018/150‑1 (MQ=255)
                                                                                                                                  aaGCGGGGTAATTGCGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCgtt                     <  1:474325/150‑3 (MQ=255)
                                                                                                                                   cgCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAt                    <  1:1457532/149‑1 (MQ=255)
                                                                                                                                   cgCGGGGTAATTGTGCGATGCACAATATCGGTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAt                    <  2:1695462/149‑1 (MQ=255)
                                                                                                                                   aGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACAGCCCGCCATAGGCGGGCCGGAt                    >  2:1049764/1‑150 (MQ=255)
                                                                                                                                    gctgggTATTTGTGCGATGCACATTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCATGAGGCGGGAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   <  1:343292/147‑1 (MQ=255)
                                                                                                                                     cgctgTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                  <  2:199623/146‑1 (MQ=255)
                                                                                                                                     caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGAGGGCCGGAttt                  >  2:136821/8‑150 (MQ=255)
                                                                                                                                     cGGGGTAATTGTGCGATGCACAATATGGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                  <  2:1632147/150‑1 (MQ=255)
                                                                                                                                      aacgataTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa                 >  2:585892/7‑150 (MQ=255)
                                                                                                                                        cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   >  1:261989/5‑146 (MQ=255)
                                                                                                                                        cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   <  2:1226679/142‑1 (MQ=255)
                                                                                                                                        cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   <  2:261989/142‑1 (MQ=255)
                                                                                                                                        cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                   >  1:1226679/5‑146 (MQ=255)
                                                                                                                                        cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa               >  2:28118/5‑150 (MQ=255)
                                                                                                                                         gataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAt              >  2:5395/4‑150 (MQ=255)
                                                                                                                                         gataTTGTGCATCGCACAATATCGTTGATGTGGTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTTTCCATACAATACGTCCCGCCATAGGCGGGCAGGATTTACAt              >  2:230770/4‑150 (MQ=255)
                                                                                                                                             ttGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGGGGAGCACGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGc          <  2:955220/150‑1 (MQ=255)
                                                                                                                                                tgCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAAATAATAAGGCGCAGAGCGTCGTTTTGAACATAGGTTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACAATGGCCAt       >  2:1513552/1‑150 (MQ=255)
                                                                                                                                                 gCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATg      >  1:750292/1‑150 (MQ=255)
                                                                                                                                                   gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg    <  1:512196/150‑1 (MQ=255)
                                                                                                                                                     tGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGtt  <  2:1182589/150‑1 (MQ=255)
                                                                                                                                                     tGCACAATAACGATGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACAGCCCGCCATAGGCGGGCCGGATTTACATTCgcgttgcgtt  >  2:624957/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTT  >  NC_000913/4378710‑4379008

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