Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 728,330 C→T V135M (GTG→ATG)  kdpA ← K(+) transporting P‑type ATPase subunit KdpA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913728,3300CT94.5% 162.8 / ‑6.2 55V135M (GTG→ATG) kdpAK(+) transporting P‑type ATPase subunit KdpA
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (21/31);  total (22/33)
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.00e+00

GGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCACCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTC  >  NC_000913/728088‑728573
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
gggTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGc                                                                                                                                                                                                                                             <  2:92225/251‑1 (MQ=255)
 ggtgaCAGCTGTTGCGCTCCGTCAACGGTAGTCACATCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCCCCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCTGCAGAAATAAAGTTTTGCTTCGTTAAGcc                                                                                                                                                                                                                                            <  2:94623/247‑1 (MQ=255)
                    cgtCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCGTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAg                                                                                                                                                                                                                         <  2:169530/249‑1 (MQ=255)
                         aCGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCaa                                                                                                                                                                                                                    <  1:548929/251‑1 (MQ=255)
                              aTTCACAGCCTGATAAGGCTGAAAGTTGTGCAGCTCACCTTGTTGAATAAAAAATAGTGCAATCAACATCGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGTTCTGGCGGGTAAACGCACGGATGATGGCAAAAATCTCCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCTTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTg                                                                                                                                                                                                                <  2:277017/250‑1 (MQ=255)
                                 cacaGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGtt                                                                                                                                                                                                             <  1:5120/250‑1 (MQ=255)
                                               gCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTTGTTTCACCGCTATAAGAt                                                                                                                                                                                              <  2:84679/251‑1 (MQ=255)
                                               gCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGGTTTGCATCGGTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCtcta                                                                                                                                                                                                  >  2:185920/1‑244 (MQ=255)
                                                         ttGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTg                                                                                                                                                                                     <  2:110816/250‑1 (MQ=255)
                                                            cAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTAt                                                                                                                                                                                 <  1:544390/251‑1 (MQ=255)
                                                               cgcACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGACATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTgg                                                                                                                                                                              <  2:77843/251‑1 (MQ=255)
                                                                       gttgAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGGGGGTTCACCGCTATAAGATTGCCCGTTGGTATTGGTGACaa                                                                                                                                                                        >  1:453846/1‑249 (MQ=255)
                                                                          gAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCt                                                                                                                                                                    <  2:224235/250‑1 (MQ=255)
                                                                            aTAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTTCCAGTTGGTATTGGTGACAAAGCTGAc                                                                                                                                                                 <  2:559990/251‑1 (MQ=255)
                                                                              aaaaaaCAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGAcgg                                                                                                                                                               >  2:443662/1‑251 (MQ=255)
                                                                                                caacGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGAAAGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCGATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGAt                                                                                                                                             <  2:20269/251‑1 (MQ=255)
                                                                                                    gCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGAtt                                                                                                                                                                                             <  1:126849/199‑1 (MQ=255)
                                                                                                    gCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGAtt                                                                                                                                                                                             >  2:126849/1‑199 (MQ=255)
                                                                                                         aGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGAc                                                                                                                                     >  2:251623/1‑250 (MQ=255)
                                                                                                                aGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCAGTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCCCCGCAAGCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAGATAGCTCAACGTGGTTTCGCCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTg                                                                                                                              <  2:377453/250‑1 (MQ=255)
                                                                                                                aGCACCCATAACGTGATGCGAAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGGAAACGCACGGATGCGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGAATGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACCACAGCCCTgg                                                                                                                             >  1:461300/1‑251 (MQ=255)
                                                                                                                   aCCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAATTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTTGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCaa                                                                                                                          <  1:366690/251‑1 (MQ=255)
                                                                                                                     ccATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAAct                                                                                                                        <  1:296385/251‑1 (MQ=255)
                                                                                                                       aTAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGGAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACCGCCCTGGCAActgc                                                                                                                      >  1:521533/1‑251 (MQ=255)
                                                                                                                               aTGCGTAGCAGATCGACCCAGGCCTTCCCGAGCGTGCTCATGCTCTGGCGGGGAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCCTCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAACGCCAGATCCCACGACAGCCCTGGCCACTGCTGGGGAtt                                                                                                              >  2:647637/1‑251 (MQ=255)
                                                                                                                                      gCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACTGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGc                                                                                                       <  1:580545/251‑1 (MQ=255)
                                                                                                                                         gATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGc                                                                                                       >  1:222964/1‑248 (MQ=255)
                                                                                                                                             gACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGGTGTGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGATAGAAAGTTTTGCATCGTTAAACCCGCCATCTGGCGGAACTAGCTCCACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTAATTAGCGCCCGATCGCACGGCAGCCCTGGCAACCGCTGTGTGTTAATCGGCAGAAAGt                                                                                                >  2:479461/1‑251 (MQ=255)
                                                                                                                                               cccAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGc                                                                                              <  2:578180/251‑1 (MQ=255)
                                                                                                                                                ccAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCt                                                                                             <  2:566009/251‑1 (MQ=255)
                                                                                                                                                 cAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTg                                                                                            >  1:10788/1‑251 (MQ=255)
                                                                                                                                                     cATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGGCCg                                                                                        >  2:325418/1‑251 (MQ=255)
                                                                                                                                                         cccGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGTCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGc                                                                                     <  2:523242/250‑1 (MQ=255)
                                                                                                                                                                  gCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGTCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATaaaaa                                                                           <  1:80180/251‑1 (MQ=255)
                                                                                                                                                                                  gTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGccc                                                           >  1:614976/1‑251 (MQ=255)
                                                                                                                                                                                     aacgcaacgATGAGGGCAAAACTCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCACCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGcccc                                                          >  2:262391/9‑249 (MQ=255)
                                                                                                                                                                                        gCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCAt                                                     >  2:266385/1‑251 (MQ=255)
                                                                                                                                                                                           cGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGtt                                                  >  1:448995/1‑251 (MQ=255)
                                                                                                                                                                                               tGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAgg                                              >  2:359428/1‑251 (MQ=255)
                                                                                                                                                                                                gAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGc                                             <  1:185920/251‑1 (MQ=255)
                                                                                                                                                                                                 aGGGCAAAAATCACCGCAAGCCGGCGGGCGGCAGAAAGAAAGTTTTGCATCTTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCTTTGGGTATTGGTGACAAAGCGGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAAGTTGTGGGGATTAAGCGGCAGATGGTGTGGACCGAGCAAAATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGcc                                            <  1:634377/251‑1 (MQ=255)
                                                                                                                                                                                                 aGGGCAAAAATCACCGCAACCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGCCCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGcc                                            <  1:325418/251‑1 (MQ=255)
                                                                                                                                                                                                  gggaaaaaaTCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCCCCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCg                                           <  1:251623/247‑1 (MQ=255)
                                                                                                                                                                                                  ggacAAAAATCACCGCACTCCCGTTGGCGGAAGAAAAAAAGTTTGGCACCGTTAAGCCCGCCATCGGGCAGAAAGAGGTAAACTGGGTTACACCGCTATAAGATGGCAAGGTGGTATTGGTGACAAAGCTGACGGCGGTATTGGGCGCCAGATCCCCCGACAGCCCGGGCGCCTGCTGTGGATTAAGCGGCAGATAGTGGTGACCGAGCAACAAAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCg                                           <  2:461300/248‑1 (MQ=255)
                                                                                                                                                                                                         aaTCACCGCAATCCCGCTGGCGGAAAAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGCCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTTACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATATTGCTGACCGAGCAACAAAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATgg                                    <  2:228720/251‑1 (MQ=255)
                                                                                                                                                                                                              ccGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGcaca                                >  2:316728/1‑250 (MQ=255)
                                                                                                                                                                                                                 cAATCCCGCTGGCGGAAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGTTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCATTTGGTTTTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGACCCCAGGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGTGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGa                            <  1:653423/251‑1 (MQ=255)
                                                                                                                                                                                                                  aaTCCCGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACaa                               >  1:511127/1‑247 (MQ=255)
                                                                                                                                                                                                                     cccGCTGGCGGCAGAAATAAAGTTTTGCATCTTTAAGCCCACAATCTGGCTGAAATAGCGAAACGTGTTTTCCCCGCTATAAGATTTCCATTTGTTATGGGTGACAAAGCGGACGGCGGTATTCAGCGCCAGACCCCACGACAGCCCTGGAATCTGCTGTGGTTTAAGCGGCAGATATTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATAtt                        <  2:192576/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccGCTGGCGGCAGAAAGAAAGTTTTGCTTCGTTAAGCCCGCCATCTGCCTGAAATAGCTCAACGTGGTTGCACCGCTATAAGATTGCCAGTTGGTTTTGGTGACAAAGCTGACGGCGGTATTCAGCCCCAGATCCCACGACAGCCCTGGAAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATAtt                        <  1:510125/251‑1 (MQ=255)
                                                                                                                                                                                                                     cccGCTGGCGGCAGAAAGAAAGTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATAtt                        >  1:4214/1‑251 (MQ=255)
                                                                                                                                                                                                                      ccGCTGGCGGCAGAAAGAAATTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGGCGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACAAAAAAAACAGGACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTg                       <  1:359428/251‑1 (MQ=255)
                                                                                                                                                                                                                                       aaggTTTTGCACCGTTAAGCCCGCCAACTGGCTGAATTAGCTCAACGTGGTTTCTCCGCTATAAGATTGCCATTTGGTTTTGGTGACAAAGCTGACGGCTGTATTCAGCCCCAGATCCCACGACAGCCCGGGAAACTGCGGTGGATTAGGCGGCAGATAGGGGTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCAc      <  2:615739/248‑1 (MQ=255)
                                                                                                                                                                                                                                       aaagTTTTGCATCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGAGACAAAGCTGACGGCGGTATTTAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAAAATAAAAAACAGCACCGCCAGCCCCAGCATGATCAGGCCGATATTGGCAAAAAGATATCGCTTACAGTTCAtctccc      >  2:654253/1‑249 (MQ=255)
                                                                                                                                                                                                                                              tGCATCGTTAAGCCAGCCATCTGGCTGAAATAGCGCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGACAGCTCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGATGACCGAGGAACATAAAAAACAGCACCGCCAGCCACAGCATGTTACGGCCGAGACTGGCACAAAGATATTGCTTCCAGGTCCTCTCAccgca  >  2:509237/1‑244 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                   
GGGTAACAGCTGTTGCGCTCCTTCAACGGTATTCACAGCCTGATAAGGCAGAAAGTTTTGCAGCGCACCTTGTTGAATAAAAAACAGTGCAATCAACAACGCCACAGGGACTAGCACCCATAACGTGATGCGTAGCAGATCGACCCAGGCATTCCCGAGCGTGCTCATGCTCTGGCGGGTAAACGCACGGATGAGGGCAAAAATCACCGCAATCCCGCTGGCGGCAGAAAGAAAGTTTTGCACCGTTAAGCCCGCCATCTGGCTGAAATAGCTCAACGTGGTTTCACCGCTATAAGATTGCCAGTTGGTATTGGTGACAAAGCTGACGGCGGTATTCAGCGCCAGATCCCACGACAGCCCTGGCAACTGCTGTGGATTAAGCGGCAGATAGTGCTGACCGAGCAACATAAAAAACAGCACCGCCAGCCCCAGCATGTTCAGGCCGAGAATGGCACAAAGATATTGCTTCCAGTTCATCTCACGGTC  >  NC_000913/728088‑728573

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑