Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,858 A→T 9.0% H29Q (CAT→CAA) ‡ 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,8580AT 9.0% 302.6 / 15.9 119H29Q (CAT→CAA) ‡ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (43/65);  new base T (6/5);  total (49/70)
Fisher's exact test for biased strand distribution p-value = 3.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.87e-01

Only 100 of 120 total aligned reads displayed.

CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCG  >  NC_000913/4378709‑4379006
                                                                                                                                                     |                                                                                                                                                    
ccTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                                                                                      >  2:1696339/1‑150 (MQ=255)
 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)
            aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAattacc                                                                                                                                          <  2:504421/150‑5 (MQ=255)
               ttttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTAGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                       <  2:887738/149‑1 (MQ=255)
                 ttGCTGGGTGCCGGGCTGCGTTTTGGCGATGTCCGCATAGCCCCATGTCAAGTAATAAGGTGTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAt                                                                                                                                     <  2:666013/150‑1 (MQ=255)
                  tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                                                                                    >  2:590024/1‑150 (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:1095751/150‑1 (MQ=255)
                        gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                              >  2:362161/1‑150 (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)
                                      tttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgt                                                                                                                >  1:798342/1‑150 (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:1561242/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)
                                                     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)
                                                              ggTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCATGAGGCGGTAATtaataa                                                                                        >  2:1097986/1‑150 (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)
                                                                        taaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                                                                              >  2:286666/1‑150 (MQ=255)
                                                                         aaGGTTTAACCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGATGCTCAATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                                                                             <  1:881137/150‑1 (MQ=255)
                                                                             tttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATAACGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcgt                                                                         >  2:697049/1‑150 (MQ=255)
                                                                              ttACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGATGATATGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGTCGCAGAGCGTCGtt                                                                        >  2:343292/1‑150 (MQ=255)
                                                                                 cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                     <  1:1406207/150‑1 (MQ=255)
                                                                                 cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                     <  1:1354827/150‑1 (MQ=255)
                                                                                   cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                   >  1:426629/1‑150 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            >  1:993622/1‑150 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            >  1:1182589/1‑150 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            >  1:1444901/1‑150 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATCCACAATATCGTTGATATGTTGAGGGGCAGCAAATCTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGACGTTTTGAACATaccg                                                            >  1:1711718/1‑147 (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                                                    >  1:414705/1‑150 (MQ=255)
                                                                                                      gccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                <  1:478740/150‑1 (MQ=255)
                                                                                                      gccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                <  1:488002/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)
                                                                                                                         gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                             <  2:712239/150‑1 (MQ=255)
                                                                                                                         gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                             <  2:577116/150‑1 (MQ=255)
                                                                                                                              tctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCggg                        <  2:1178712/150‑1 (MQ=255)
                                                                                                                               ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                       >  2:391983/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)
                                                                                                                                     gctgggTATTTGTGCGATGCACATTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCATGAGGCGGGAATTAATAAGGCGCAGAGCGTCGTGTTGAACATAGGGTCTGGTATCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 <  1:343292/147‑1 (MQ=255)
                                                                                                                                     gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 <  2:1654857/150‑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                 <  2:261989/142‑1 (MQ=255)
                                                                                                                                         cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 >  1:261989/5‑146 (MQ=255)
                                                                                                                                         cgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 <  2:1226679/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)
                                                                                                                                                   cGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGCTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGc   <  1:286666/150‑1 (MQ=255)
                                                                                                                                                    gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg  <  1:512196/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCG  >  NC_000913/4378709‑4379006

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: