Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,912,748 T→G 5.1% intergenic (‑172/+20) htpX ← / ← prc putative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,912,7480TG 5.1% 391.7 / 8.1 131intergenic (‑172/+20)htpX/prcputative endopeptidase/carboxy‑terminal protease for penicillin‑binding protein 3
Reads supporting (aligned to +/- strand):  ref base T (60/63);  new base G (4/3);  total (65/66)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 136 total aligned reads displayed.

AACAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912599‑1912896
                                                                                                                                                      |                                                                                                                                                    
  cAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAAttt                                                                                                                                                    >  1:1389168/1‑150 (MQ=255)
    aaaCGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCt                                                                                                                                                  <  2:613962/150‑1 (MQ=255)
        gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTtgtg                                                                                                                                              >  1:1003061/1‑150 (MQ=255)
        gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTtgtg                                                                                                                                              <  1:1052486/150‑1 (MQ=255)
        gCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTtgtg                                                                                                                                              <  2:1156248/150‑1 (MQ=255)
         cAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATATACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGc                                                                                                                                             <  2:726647/150‑1 (MQ=255)
              tGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAAGTTCTTGTGCCtgat                                                                                                                                        >  1:1444730/1‑150 (MQ=255)
               gCGATACCCACATCGTATGGTTTACGCGACTATTTTCAAGTCTGGATAGTGCTTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCtgatt                                                                                                                                       <  1:159704/150‑1 (MQ=255)
               gCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCtgatt                                                                                                                                       <  2:508721/150‑1 (MQ=255)
                cGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCATGCACAATTTCTTGTGCCtgattg                                                                                                                                      >  2:605453/1‑150 (MQ=255)
                 gATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTATACATTTTGTAGCATCTGATTTACGGCATCTTGCCGCTGTT‑AAAAAATCAGGCACAAGAAATTGTGCCtgatttt                                                                                                                                     <  1:81481/150‑3 (MQ=255)
                     cccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATAtt                                                                                                                                 >  1:1063889/1‑150 (MQ=255)
                     cccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATAtt                                                                                                                                 >  2:469735/1‑150 (MQ=255)
                     cccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCCCTGTT‑AAAAAATCAGGCACAATTTCTAGTGCCTGATTGATAtt                                                                                                                                 >  2:1615230/1‑150 (MQ=255)
                      ccACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTa                                                                                                                                >  1:549877/1‑150 (MQ=255)
                       cacaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTAc                                                                                                                               <  1:243024/150‑1 (MQ=255)
                        acaTCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACt                                                                                                                              <  1:1052725/150‑1 (MQ=255)
                            cGTATGGGTTACGCGACTATTTTCAATTCTGGCTAGTGCGTAAGTCAGCAGAAAGACACAACTTTAATTTTTGTAGCATCTGATTTACGGCATGTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                          <  2:315155/150‑1 (MQ=255)
                            cGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                          <  1:910816/150‑1 (MQ=255)
                            cGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGAc                                                                                                                          <  1:882834/150‑1 (MQ=255)
                             gTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACg                                                                                                                         >  2:914049/1‑150 (MQ=255)
                                      aCGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGgt                                                                                                                <  2:794246/150‑1 (MQ=255)
                                            cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTcc                                                                                                          >  1:1518447/1‑150 (MQ=255)
                                            cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTcc                                                                                                          >  2:843209/1‑150 (MQ=255)
                                            cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTcc                                                                                                          >  1:529479/1‑150 (MQ=255)
                                            cTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGA‑TTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCATGATTGATATTACTTGACGGGAGCGGGTTGTTCcg                                                                                                         >  2:338484/1‑150 (MQ=255)
                                               ttttCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACGTTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCcgcg                                                                                                       >  1:1535266/1‑150 (MQ=255)
                                                   cAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGGGAGGGTTGTTCCGCGGGTc                                                                                                   >  1:164816/1‑150 (MQ=255)
                                                   cAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTc                                                                                                   >  1:1344164/1‑150 (MQ=255)
                                                     aGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTg                                                                                                 >  1:1208257/1‑150 (MQ=255)
                                                        cTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCCCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCt                                                                                              <  1:1320161/150‑1 (MQ=255)
                                                          ggATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCttt                                                                                            >  1:1003252/1‑150 (MQ=255)
                                                            aTAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCttttt                                                                                          >  2:785195/1‑150 (MQ=255)
                                                              aGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCa                                                                                        >  1:304955/1‑150 (MQ=255)
                                                               gTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCACGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCaa                                                                                       <  1:21910/150‑1 (MQ=255)
                                                                tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      >  2:1900479/1‑150 (MQ=255)
                                                                tGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAg                                                                                      <  1:162080/150‑1 (MQ=255)
                                                                   gTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGCCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCtt                                                                                   <  1:1019776/150‑1 (MQ=255)
                                                                      aGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGc                                                                                <  1:666566/150‑1 (MQ=255)
                                                                      aGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGAACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGAGGGTTGTTCCGCGGGACTGGCTTTTTCAAGCTTCGc                                                                                >  2:799226/1‑150 (MQ=255)
                                                                       gTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTc                                                                                         <  2:516069/140‑1 (MQ=255)
                                                                       gTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGcc                                                                               <  2:1074904/150‑1 (MQ=255)
                                                                       gTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGcc                                                                               <  1:1277385/150‑1 (MQ=255)
                                                                       gTCACCAGAAAGACACAACTTTACATTATGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATAATACTTGACGGCAGCGGGATGTTCCGCGGGTCTGCCTTTTTCAAGCTTCGcc                                                                               >  2:1212104/1‑150 (MQ=255)
                                                                        tCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCa                                                                              <  2:544492/150‑1 (MQ=255)
                                                                            cAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATc                                                                          >  2:1545667/1‑150 (MQ=255)
                                                                                aaGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGt                                                                      >  1:978110/1‑150 (MQ=255)
                                                                                 aGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTg                                                                     <  1:1095456/150‑1 (MQ=255)
                                                                                  gACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGATCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGc                                                                    >  1:1687201/1‑150 (MQ=255)
                                                                                  gACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGAGc                                                                    >  2:1680906/1‑150 (MQ=255)
                                                                                  cacacacCTATACATTTTGTACCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAACTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTGGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGc                                                                    <  2:74982/149‑1 (MQ=255)
                                                                                    cacaACTTTACTTTTTTTATCATCTGATTTACGGCTTCTTGTCGCTGTT‑AAAACATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGGTGTTCCGCGGGTCTTGCTTTTTCAAGCTTCGCCAGATCTAGTGCgt                                                                  <  1:1093784/150‑2 (MQ=255)
                                                                                    cacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  >  1:1314909/1‑150 (MQ=255)
                                                                                    cacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  >  2:1460612/1‑150 (MQ=255)
                                                                                    cacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAAATTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  >  1:1801801/1‑150 (MQ=255)
                                                                                    cacaAATTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  <  1:1277243/150‑1 (MQ=255)
                                                                                    aacaACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGa                                                                  <  1:629503/149‑1 (MQ=255)
                                                                                      caACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGata                                                                <  1:1619227/150‑1 (MQ=255)
                                                                                      caACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGata                                                                >  2:1689859/1‑150 (MQ=255)
                                                                                        aCTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATAtt                                                              <  2:965571/150‑1 (MQ=255)
                                                                                         cTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTc                                                             <  2:794743/150‑1 (MQ=255)
                                                                                          tttACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCa                                                            >  2:1640118/1‑150 (MQ=255)
                                                                                            tACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCAcc                                                          <  2:1389168/150‑1 (MQ=255)
                                                                                             aCATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCg                                                         <  2:1067248/150‑1 (MQ=255)
                                                                                                ttttGTAGCATCTGATTTACGGCTTCTTGTCGCTGTT‑AAAAATTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGTGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  1:1026763/150‑1 (MQ=255)
                                                                                                ttttGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGtct                                                      <  2:1063889/150‑1 (MQ=255)
                                                                                                    gTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGTCACAATTTATTGTGCCTGATTGATATTAATAGACAGGAGCGGGTTGTTACGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATc                                                  >  1:1146511/1‑150 (MQ=255)
                                                                                                     tAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATcc                                                 <  2:71860/150‑1 (MQ=255)
                                                                                                      aGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                                >  1:795229/1‑150 (MQ=255)
                                                                                                      aGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCa                                                <  2:1003252/150‑1 (MQ=255)
                                                                                                         aTCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAt                                             >  1:1568256/1‑150 (MQ=255)
                                                                                                         aTCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGAt                                             <  2:196038/150‑1 (MQ=255)
                                                                                                            tGATTTACGGCATCTTGTCGCTGTTAAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                           <  1:1401116/150‑1 (MQ=255)
                                                                                                            tGATTTACGGCATCTTGTCGCTGTTAAAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATaa                                           <  1:1401109/150‑1 (MQ=255)
                                                                                                             gATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAgg                                         >  1:602235/1‑150 (MQ=255)
                                                                                                             gATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAgg                                         >  2:393087/1‑150 (MQ=255)
                                                                                                                   cGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCgg                                   >  1:1880620/1‑150 (MQ=255)
                                                                                                                    ggCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTTGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGc                                  <  1:914049/150‑1 (MQ=255)
                                                                                                                     gCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCt                                 >  1:1846092/1‑150 (MQ=255)
                                                                                                                       aTCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCGTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTTTTCCGCGGCTCTGGCTTTTTCAAGCTTCGCCAGAGCGCGTGCGATATTCACCGTCTCCTCCAGCTAAGGATCCGGCTcc                               >  1:234013/1‑150 (MQ=255)
                                                                                                                          ttGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTgg                            <  2:1649967/150‑1 (MQ=255)
                                                                                                                           tGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGt                           <  2:1003061/150‑1 (MQ=255)
                                                                                                                           tGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGt                           >  1:435715/1‑150 (MQ=255)
                                                                                                                              cGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGACAGATCGAGTGCGATATTCACCGTCTCATCCAGATAATGATCCGCCTCCTGGTAAt                        >  2:349217/1‑150 (MQ=255)
                                                                                                                               gCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAAtg                       >  1:748139/1‑149 (MQ=255)
                                                                                                                                   tt‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTc                   <  2:1437747/150‑1 (MQ=255)
                                                                                                                                    t‑AAAAACTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                  <  1:1672187/150‑1 (MQ=255)
                                                                                                                                    t‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCg                  <  1:518458/150‑1 (MQ=255)
                                                                                                                                      aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCgg                 <  2:800478/150‑1 (MQ=255)
                                                                                                                                      aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGAGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCgg                 >  1:1855394/1‑150 (MQ=255)
                                                                                                                                      aaaaaaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGCCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCCCCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCgg                 >  1:220517/1‑150 (MQ=255)
                                                                                                                                       ataaaTCTGGCACTATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGt                <  1:902980/148‑1 (MQ=255)
                                                                                                                                          aaTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGa             <  2:230992/150‑1 (MQ=255)
                                                                                                                                           aTCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            <  1:436423/150‑1 (MQ=255)
                                                                                                                                           aTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            >  1:120636/1‑150 (MQ=255)
                                                                                                                                           aTCAGGCACAAGAAATTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGat            >  2:41955/1‑150 (MQ=255)
                                                                                                                                            tCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatc           <  2:1444730/150‑1 (MQ=255)
                                                                                                                                            tCAGGCACAATTTCTTGTGCCGGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCGTCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatc           <  2:1438259/150‑1 (MQ=255)
                                                                                                                                             cAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGatca          <  1:831155/150‑1 (MQ=255)
                                                                                                                                                     aTTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCTAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGttt  >  1:222588/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                    
AACAAAACGCAATATGCGATACCCACATCGTATGGGTTACGCGACTATTTTCAAGTCTGGATAGTGCGTAAGTCACCAGAAAGACACAACTTTACATTTTGTAGCATCTGATTTACGGCATCTTGTCGCTGTT‑AAAAAATCAGGCACAATTTCTTGTGCCTGATTGATATTACTTGACGGGAGCGGGTTGTTCCGCGGGTCTGGCTTTTTCAAGCTTCGCCAGATCGAGTGCGATATTCACCGTCTCATCCAGATAAGGATCCGGCTCCTGGTAATCTTTCGGTAGATCATCCAGTTT  >  NC_000913/1912599‑1912896

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