Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,859 T→C 6.7% 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 6.7% 348.7 / 12.0 120H29R (CAT→CGT) ‡ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base T (54/57);  new base C (4/4);  total (59/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.31e-01

Only 100 of 120 total aligned reads displayed.

CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCG  >  NC_000913/4378714‑4379006
                                                                                                                                                 |                                                                                                                                                   
cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                                                                                 <  2:554734/150‑1 (MQ=255)
  cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                                                                               >  2:743061/1‑150 (MQ=255)
    aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGTAGTCCGCATAGCACCAGGTAAAGTAATAAGGTTTACGCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAatt                                                                                                                                             <  1:977041/150‑2 (MQ=255)
     acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAatta                                                                                                                                            <  1:464227/150‑3 (MQ=255)
       aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGATTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                                                                          <  1:1021322/150‑1 (MQ=255)
        acGTTAGCTGTGGGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAACGGTTACCCTGATAAATTACCGCCACCGCCATCCCCGGGATCTTTAGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAattaccc                                                                                                                                         <  1:117195/150‑6 (MQ=255)
         cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  2:1629694/150‑1 (MQ=255)
         cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  1:778887/150‑1 (MQ=255)
         cGTTTGCTGTGTGACGGGCTGATTTTTGGCGATGTCCGCAAAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  2:108519/150‑1 (MQ=255)
         agTTTGATGTGTGACGTGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGTTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCAATGCCATACCCGGGATCTTTTGTTTCTCTATAAGCGGGGTAATTGTGCGATGCAGAATATCGtt                                                                                                                                        <  2:63316/149‑1 (MQ=255)
          gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                       <  2:179506/150‑1 (MQ=255)
          gTTTGCTGGGTGACGGGCTGCCTTTTGGCGCTGGCCGCATATCCCCAGGTAAAGTAATAAGGGTTACCCTGATAAATTACCGCCACCGCCAAACCCGGGATCCTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                       <  1:621624/150‑1 (MQ=255)
           tttGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGa                                                                                                                                      <  1:1730602/150‑1 (MQ=255)
             tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAtt                                                                                                                                    <  1:808651/150‑1 (MQ=255)
                 gtgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCTGGGTAATTGTGCGTTGCACAATATCGTTGATttgtt                                                                                                                                >  1:253359/1‑150 (MQ=255)
                 gtgtgACCGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                                                                                <  1:1456133/150‑1 (MQ=255)
                   gtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                              <  1:617484/150‑1 (MQ=255)
                   gtgACGGGCGGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                                                                              <  2:502336/150‑1 (MQ=255)
                    tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCACCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                                                                             <  1:1300502/150‑1 (MQ=255)
                     gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                            >  2:407001/1‑150 (MQ=255)
                     gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                            <  1:775526/150‑1 (MQ=255)
                       cGGGCGGCTTTGTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATCAGGTTTACCCTGATAATTTACCGCCACCGCCATACCCGGGAGCTTTTGGTGCTCTATACGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                                                                          <  2:1524136/150‑1 (MQ=255)
                          gctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                                                                       >  2:788314/1‑150 (MQ=255)
                                  ttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTAGCCTGATAAATTACCGCCACCGCCATACCCGGGGTCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAGGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                                                                               >  1:110613/1‑150 (MQ=255)
                                  ttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                                                                               <  2:1808555/150‑1 (MQ=255)
                                  ttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTGTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgtg                                                                                                               <  2:1421815/150‑1 (MQ=255)
                                   tGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                                                                              <  2:1835712/150‑1 (MQ=255)
                                   tGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                                                                              >  2:1365254/1‑150 (MQ=255)
                                    ggCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                                                                             <  2:1474979/150‑1 (MQ=255)
                                      cGATTTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTATATGTGCGATGCACAATATCGTTGATTTGTTGAGGGGGAGCAAATGTGGAGc                                                                                                           >  1:349790/1‑150 (MQ=255)
                                       gATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                                                                          <  2:1828919/150‑1 (MQ=255)
                                          gTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                                                                       >  2:1581973/1‑150 (MQ=255)
                                            ccGCATAGCCCCAGGTAAAGTAATAAGGTTTTACCTGATAAATTACCGCCACCGCCATCCCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                                                                                     <  1:131825/150‑1 (MQ=255)
                                               cATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                                                                                  <  1:1675886/150‑1 (MQ=255)
                                                   gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATCCCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                                                                              <  1:1654784/150‑1 (MQ=255)
                                                   gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                                                                              >  1:663341/1‑150 (MQ=255)
                                                   gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                                                                              >  1:327520/1‑150 (MQ=255)
                                                   gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACAGCCACCGCCCTAGCCTGGCTCGGTTGTTGCTCTATACGCGGTGTAATTGTCCGATGCACAATATCGGTGATTTGTTGAGGGGCAGCAAATGTGGATCAAGAGTCGGtaat                                                                                              >  1:87090/1‑150 (MQ=255)
                                                            aaaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                                                                                     <  2:681958/150‑1 (MQ=255)
                                                             aaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                                                                                    >  2:1970669/1‑150 (MQ=255)
                                                             aaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGCTAATTGTGCGATGCACAATATAGTTGAATTTTTGAGTGGCAACAAATGTGTAGCAAGAGGCGCTAATTAATAAgccc                                                                                    >  2:341420/1‑147 (MQ=255)
                                                              aGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                                                                                   >  1:1523260/1‑150 (MQ=255)
                                                               gTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCa                                                                                  >  1:1050142/1‑150 (MQ=255)
                                                                taataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                                                                                 >  2:1479850/1‑150 (MQ=255)
                                                                taataaCGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                                                                                 <  1:1845656/150‑1 (MQ=255)
                                                                 aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:1126629/1‑150 (MQ=255)
                                                                 aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:628330/1‑150 (MQ=255)
                                                                  ataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                                                                               <  2:1895383/150‑1 (MQ=255)
                                                                   taaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                                                                              >  1:1956102/1‑150 (MQ=255)
                                                                      ggTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTATTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                                                                           >  2:1734167/1‑150 (MQ=255)
                                                                            cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATACGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                     <  1:891287/150‑1 (MQ=255)
                                                                            cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                     <  1:6891/150‑1 (MQ=255)
                                                                              cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                   >  1:1095939/1‑150 (MQ=255)
                                                                               tGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc                                                                  >  1:557086/1‑150 (MQ=255)
                                                                                  tAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATa                                                               >  2:1287976/1‑150 (MQ=255)
                                                                                    aaTTACCGCCACCGCTATACCCGGGATCTTTTGTTGCTCTATAAGCTGTGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGTACATAgg                                                             <  1:1682969/150‑1 (MQ=255)
                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATTCACAATATCGTTGATTTGGTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGTCGCATAGCGTCGTTTTGAACATAggg                                                            >  1:833074/1‑150 (MQ=255)
                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGCACATAggg                                                            >  2:1311546/1‑150 (MQ=255)
                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            <  2:1994451/150‑1 (MQ=255)
                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            >  1:704249/1‑150 (MQ=255)
                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTATATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTTTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGTTAATTAATAAGGCGCAGAGCGTCGGTTTGAACATAggg                                                            >  1:1450644/1‑150 (MQ=255)
                                                                                         ccgccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg                                                        >  1:1212445/1‑150 (MQ=255)
                                                                                            ccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGtt                                                     >  1:1872465/1‑150 (MQ=255)
                                                                                            ccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGtt                                                     >  2:1624237/1‑150 (MQ=255)
                                                                                             caccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGttt                                                    >  1:1864562/1‑150 (MQ=255)
                                                                                             caccgccaTACCCGGGATCTTTTGTTGCCCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGttt                                                    <  2:1439111/150‑1 (MQ=255)
                                                                                              accgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGGCGTTTTGAACATAGGGTCTGGTTTc                                                   <  1:703513/150‑1 (MQ=255)
                                                                                               ccgccaTACGCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTTTGCGATGCACAATATCGTTGATTTGTTGAGGCGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  <  2:932888/150‑1 (MQ=255)
                                                                                               ccgccaTACCCGGGATCTTTTGTTGCTCTATTAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCATATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  <  2:744421/150‑1 (MQ=255)
                                                                                               ccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGTGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  >  1:929713/1‑150 (MQ=255)
                                                                                               ccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  <  1:1235371/150‑1 (MQ=255)
                                                                                               ccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  <  1:1268546/150‑1 (MQ=255)
                                                                                               cccccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGGAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTcc                                                  <  2:303585/147‑1 (MQ=255)
                                                                                                 cccaaccccTGGCTCTTTTCTTGCTCTAAAACCGGGGAAATTGTGCGCTGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGATCAAGCGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCAt                                                <  2:1327674/144‑1 (MQ=255)
                                                                                                  ccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGGAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATa                                               <  2:1976476/150‑1 (MQ=255)
                                                                                                   caTACCCGCGATCTTTTGTTGCTCTATAAGCGGGCTAATTGTGCGCTGCACAATATCGTTGATTTGTTGAGGTGCAGCAAATATGCAGCAAGAGGCGCTAATTAATAAGGCTCAGAGCGTCGTATTGAACATAGGGACTGGATTCCATAc                                              >  2:728907/1‑150 (MQ=255)
                                                                                                     tACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACaa                                            >  2:774848/1‑150 (MQ=255)
                                                                                                         cGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTAGTTGAGGTGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTAGTTTTGAACATAGGGTCTGGTTTCAATAAAAAACg                                        >  1:40018/1‑150 (MQ=255)
                                                                                                             aTCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGccc                                    >  1:1027619/1‑150 (MQ=255)
                                                                                                                    gttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAg                             >  2:291900/1‑150 (MQ=255)
                                                                                                                     ttgCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAgg                            >  2:1062602/1‑150 (MQ=255)
                                                                                                                       gCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCg                          >  2:102414/1‑150 (MQ=255)
                                                                                                                        ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                         >  2:707188/1‑150 (MQ=255)
                                                                                                                          ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGAGTCGTTTTGAACATAGGGTCTAGTTTCCATACAAAACGGCCCGCCAAAGGCGGGc                       >  2:1007493/1‑150 (MQ=255)
                                                                                                                           tataAGCGGGGGAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGcc                      <  2:313586/150‑1 (MQ=255)
                                                                                                                            ataCGCGGGGTAAATGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTATGGCGGGCCg                     <  1:1827508/150‑1 (MQ=255)
                                                                                                                            ataAGCGGGGTAATTTTGCGATGCACAATATGGTTGATTTGTTGAGGGGCAGCAAATGTGGTGCAAGAGGCGGTAATAAATAAGGCGCTGAGCGTGGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCTTTGGGGGGCCg                     <  2:1930782/150‑1 (MQ=255)
                                                                                                                            ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                     <  1:1022785/150‑1 (MQ=255)
                                                                                                                                gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 >  2:886356/1‑150 (MQ=255)
                                                                                                                                 caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                >  1:607227/8‑150 (MQ=255)
                                                                                                                                 cGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                <  1:788314/150‑1 (MQ=255)
                                                                                                                                  ggggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa               <  1:1833029/150‑1 (MQ=255)
                                                                                                                                  aacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGTCGGGCCGGATTTa               >  2:1236140/7‑150 (MQ=255)
                                                                                                                                  aacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAACAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTa               >  1:1019622/7‑150 (MQ=255)
                                                                                                                                    ggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACa             <  2:61883/150‑1 (MQ=255)
                                                                                                                                       aaTTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTg          <  2:817008/150‑1 (MQ=255)
                                                                                                                                         ttGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGc        >  1:69522/1‑150 (MQ=255)
                                                                                                                                         ttGTGCGATGAACAATATCGTTGATTTTTTGAGGGGCAGCAAATGTGGAGCAAGAGGGGGTAATTAATAAGGCGCAGAGCTTCGTTTTGAACATAGTGTCTGGTTTCCATACAATACGGCCCGCCATAGGCGGGCCGGATTTACATTGGc        <  2:1246046/150‑1 (MQ=255)
                                                                                                                                          tgtgCGATGCACAATATCGTTGAGTTGTGGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATAAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCg       <  2:1216447/150‑1 (MQ=255)
                                                                                                                                               gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg  >  1:836743/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCG  >  NC_000913/4378714‑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: