Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,857 G→A 6.7% R30C (CGC→TGC)  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,8570GA 6.7% 357.1 / 10.3 119R30C (CGC→TGC) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base G (54/57);  new base A (4/4);  total (58/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 = 7.05e-01

Only 100 of 121 total aligned reads displayed.

CCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCG  >  NC_000913/4378709‑4379006
                                                                                                                                                    |                                                                                                                                                     
ccTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGa                                                                                                                                                      >  1:681958/1‑150 (MQ=255)
     cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                                                                                 <  2:554734/150‑1 (MQ=255)
       cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACaa                                                                                                                                               >  2:743061/1‑150 (MQ=255)
         aacaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGTAGTCCGCATAGCACCAGGTAAAGTAATAAGGTTTACGCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAatt                                                                                                                                             <  1:977041/150‑2 (MQ=255)
          acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcatcgcac                                                                                                                                            <  2:681900/150‑7 (MQ=255)
          acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAatta                                                                                                                                            <  1:464227/150‑3 (MQ=255)
            aacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGATTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                                                                          <  1:1021322/150‑1 (MQ=255)
             acGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAattaccc                                                                                                                                         <  1:714064/150‑6 (MQ=255)
             acGTTAGCTGTGGGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAACGGTTACCCTGATAAATTACCGCCACCGCCATCCCCGGGATCTTTAGTTGCTCTATAAGCGGGGTAATTGTGCATCGCACAattaccc                                                                                                                                         <  1:117195/150‑6 (MQ=255)
              cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  1:778887/150‑1 (MQ=255)
              cGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  2:1629694/150‑1 (MQ=255)
              cGTTTGCTGTGTGACGGGCTGATTTTTGGCGATGTCCGCAAAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                                                                                        <  2:108519/150‑1 (MQ=255)
              agTTTGATGTGTGACGTGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGTTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCAATGCCATACCCGGGATCTTTTGTTTCTCTATAAGCGGGGTAATTGTGCGATGCAGAATATCGtt                                                                                                                                        <  2:63316/149‑1 (MQ=255)
               gTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTg                                                                                                                                       <  2:766373/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)
                         tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCACCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAg                                                                                                                             <  1:1300502/150‑1 (MQ=255)
                          gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                            >  2:407001/1‑150 (MQ=255)
                          gACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                                                                            <  1:775526/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)
                                         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:327520/1‑150 (MQ=255)
                                                        gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaat                                                                                              >  1:663341/1‑150 (MQ=255)
                                                        gCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACAGCCACCGCCCTAGCCTGGCTCGGTTGTTGCTCTATACGCGGTGTAATTGTCCGATGCACAATATCGGTGATTTGTTGAGGGGCAGCAAATGTGGATCAAGAGTCGGtaat                                                                                              >  1:87090/1‑150 (MQ=255)
                                                                 aaaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                                                                                     >  1:633642/1‑150 (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)
                                                                    gTAATAAGGTTTACCCTGATAAATTAACGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgaa                                                                                  >  2:100071/1‑148 (MQ=255)
                                                                     taataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                                                                                 >  2:1479850/1‑150 (MQ=255)
                                                                     taataaCGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                                                                                 <  1:1845656/150‑1 (MQ=255)
                                                                      aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:628330/1‑150 (MQ=255)
                                                                      aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:1887712/1‑150 (MQ=255)
                                                                      aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:1126629/1‑150 (MQ=255)
                                                                      aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                >  2:1134482/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)
                                                                                 cccTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTg                                                                     <  1:6891/150‑1 (MQ=255)
                                                                                   cTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGaa                                                                   >  1:1095939/1‑150 (MQ=255)
                                                                                    tGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAc                                                                  >  1:557086/1‑150 (MQ=255)
                                                                                     gATAAATAACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACa                                                                 <  2:410827/150‑1 (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                                                            >  1:704249/1‑150 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            <  2:2002143/150‑1 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg                                                            <  2:1994451/150‑1 (MQ=255)
                                                                                          aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTATATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTTTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGTTAATTAATAAGGCGCAGAGCGTCGGTTTGAACATAggg                                                            >  1:1450644/1‑150 (MQ=255)
                                                                                              ccgccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTg                                                        >  1:1212445/1‑150 (MQ=255)
                                                                                                gccaccgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGt                                                      >  1:836297/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)
                                                                                                    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)
                                                                                                                  aTCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGccc                                    >  1:1027619/1‑150 (MQ=255)
                                                                                                                          ttgCTCTACAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATTCAAAACGGCCCGCCATAgg                            >  1:437831/1‑150 (MQ=255)
                                                                                                                            gCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCg                          >  2:102414/1‑150 (MQ=255)
                                                                                                                             ctctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCgg                         >  2:707188/1‑150 (MQ=255)
                                                                                                                               ctATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                       >  2:1321493/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)
                                                                                                                                 ataAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCg                     <  1:1022785/150‑1 (MQ=255)
                                                                                                                                     gCGGGGTAATTGTGCGATGCACAATATCGTTGTTTTGTTGAGGGGCAGCAAATGTGGCGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 <  2:1477478/150‑1 (MQ=255)
                                                                                                                                     gCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAtt                 >  2:886356/1‑150 (MQ=255)
                                                                                                                                      caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGTCGGTAATTAATAAGTCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                >  2:1709029/8‑150 (MQ=255)
                                                                                                                                      caacgataTTGTGCATCGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGAttt                >  1:607227/8‑150 (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)
                                                                                                                                               tgtgCGATGCACAATATCGTTGAGTTGTGGCGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATAAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCg       <  2:1216447/150‑1 (MQ=255)
                                                                                                                                                    gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCg  >  1:836743/1‑150 (MQ=255)
                                                                                                                                                    gATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCAAAGGCGGGCCGGATTTACATTGGCGATGCg  >  1:1879751/1‑150 (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: