Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,378,891 C→A 17.0% A19S (GCT→TCT)  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,8910CA17.0% 137.1 / 4.7 54A19S (GCT→TCT) ampCpenicillin‑binding protein, beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base C (32/12);  new base A (2/8);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 3.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG  >  NC_000913/4378656‑4379124
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTACTAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCCTACCCGGGATCTTTTGTTGTCCTATAAGCGGGGTAATTGTGCGATGAACAAAATGTTTTATTTGTTGAGGGGCCGCaaa                                                                                                                                                                                                                                                 >  1:486529/1‑239 (MQ=255)
         aTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAATTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTCTAAGCGGGGTAATTGTGCGATGTACATAATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                                                                                                                                                                                                        >  2:421821/1‑239 (MQ=255)
               tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTTGGCGGTGCACAATATCGTTGTTTTTTTTAGGGGCAGAAAATGTGGGGAAAAAggg                                                                                                                                                                                                                                  >  1:95819/1‑238 (MQ=255)
                        aGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCGGAAAATGTGGAGAAAGAGGGGGtaattaa                                                                                                                                                                                                                        >  2:908188/1‑239 (MQ=255)
                             gCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                                                                                                                                                                                                                   >  2:778614/1‑239 (MQ=255)
                                         tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGAAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCCCTCTAAACGGGGTAATTGTGCGATCCACAATTTCGTTGATTTGTTGCGGGGCCGCAAATGTGGAGCAAGAGGGGGGAATTTAAAAAGGCGAGAGcg                                                                                                                                                                                                          >  1:22080/1‑236 (MQ=255)
                                         tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTATTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGTTTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGGGGGAATtaaaaaagggaaaggggggg                                                                                                                                                                                                       >  1:881724/1‑225 (MQ=255)
                                                                acaacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGAACTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGACCAATATCGTTGATTTGTTGAGGGGCGGCAAATTGGGAGCAAGAGGCTGGAATTTAAAAAGGTCCGAGCCCCGTTTTGacacaagagcttgggtt                                                                                                                                                                                 >  1:16754/1‑222 (MQ=255)
                                                                        tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                                                                                                                                                                                                        >  2:536942/1‑191 (MQ=255)
                                                                                    ggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGCCGTTTTGAACAAAGGGCCTGGTTTCCTTACAAACCGGACCGccta                                                                                                                                                           >  2:268970/1‑237 (MQ=255)
                                                                                       tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCGGAGCGTCGTTTTGAACATAGGGACTGGTTTCCATACAAAACGGCCCGCCATAgg                                                                                                                                                        >  2:96403/1‑239 (MQ=255)
                                                                                       tgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGGTTCCATACAAAACGGCCCGCCATAgg                                                                                                                                                        >  2:829904/1‑239 (MQ=255)
                                                                                         ctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCCTAGCCg                                                                                                                                                     >  2:653885/1‑239 (MQ=255)
                                                                                              ttgggcgatcaccATAGCCCCCGATAAAAGAAAATTGTTTTCCCTAATAAAAATCCCCGCCCCCGCCTTCCCCGGGATCTTTTGTTGCTCTAAAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCTTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGc                                                                                                                                                 <  2:881724/227‑1 (MQ=255)
                                                                                                                  cAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                                                                                                                                                                                                             >  2:368717/1‑160 (MQ=255)
                                                                                                                    ggTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCCTAGGCCGGCCGGATTTTCATTgggggatgggtt                                                                                                                        >  2:638491/1‑229 (MQ=255)
                                                                                                                    ggTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGGTTTTAACATAGGGCCTGGTTTCCATACAAAAAGGCCCGCCCTAggcggggcggggtttcactggggagacgtgt                                                                                                                        >  1:330996/1‑218 (MQ=255)
                                                                                                                         aGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTGTCCATACAAAACAGCCCGCCCTAGGCGGGCCGGGTTTCACTTGGCCAGGCGTTAAATTg                                                                                                                  >  2:928860/1‑239 (MQ=255)
                                                                                                                          gTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGGTTCCATACAAAACGGCCCGCCCTAGGCGGGGCCGGTTTTCAATTgtgcagtggttaaattt                                                                                                                   >  2:373951/1‑223 (MQ=255)
                                                                                                                              taaGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAAGCGCAGGGCGGCCTTTTGGACATAGGGGCTGGGTTCCATACAAAACGGCCCCCCCAAAGCGGGCCCGGATTTCAATTgcgaggcgagaaatggaaaaa                                                                                                               >  1:511652/1‑219 (MQ=255)
                                                                                                                                     tACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGGCGGGGTTTCAATGTCGATGTGATAGAAGGGAACAAGACCaa                                                                                                      >  1:970221/1‑239 (MQ=255)
                                                                                                                                     tACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAAAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTTCATTTGCGGAGCGTTGGATTGTTAAGAccccca                                                                                                      >  2:678793/1‑237 (MQ=255)
                                                                                                                                     tACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAAATAGGGTCTGGTTTCCCTACAAAACGGCCCGCCATAGGCGGGCCGGGATTTCCATTGGGGGTGGTTTAGAtttacagcacccaa                                                                                                       >  2:926589/1‑227 (MQ=255)
                                                                                                                                     tACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGACGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCCGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGGTTTTAATTTGCGATTCgatcatttttaaacaaccaaa                                                                                                       >  1:288296/1‑219 (MQ=255)
                                                                                                                                         cttatttattgccccccaccccccaTCCCCGGGTGTTTTTTTTGCTTTATAAGCGGGGAAATTGGGGGATGCACAATATCGTTGATTGGGTGAGGGGCAGCAAATGGGGAGCAAGAGGTGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTTTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATc                                                                                                    <  1:542961/217‑1 (MQ=255)
                                                                                                                                            aT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTAAATTTGCGGGTGGTTAGATTGGaaacaaccccaacaacgg                                                                                               >  2:287068/1‑223 (MQ=255)
                                                                                                                                             t‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTTGCGATGCGGTAGATTGTAAAGAGACCAAACAGCGTg                                                                                              >  2:205425/1‑239 (MQ=255)
                                                                                                                                                aatctccccccccggccaTACCCGGGTTTTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGa                                                                                             <  2:330996/225‑1 (MQ=255)
                                                                                                                                                  ttccccccccccgcctTCCCCGGGTTTTTTTGTTGCTTAATAAGCGGGGTAATTGGGCGAAGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAGGTGGAGCAAGAGGCGTTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATGGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACa                                                                                          <  1:205425/229‑1 (MQ=255)
                                                                                                                                                       ccaccccgccaTCCCCGGGTTTTTTTTTTGCTTTAAAAGGGGGTTAATTGGGGGATGCAAAATGTTGTTTATTTGTTGGGGGGAAGCAAAGGGGGATCAGGGGGCGGTAATTAAAAAGGCGCAGAGGGTCTTTTTGAACATGGGGTCGGTTTCCAATACAAAAGGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGt                                                                                     <  1:638491/234‑1 (MQ=255)
                                                                                                                                                             cggcaTCCCCGGGATTTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACATTATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                               <  2:486529/236‑1 (MQ=255)
                                                                                                                                                             cgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt                                                                                                                                                                                             <  2:430276/135‑1 (MQ=255)
                                                                                                                                                             cgccaTACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAAAAt                                                                                                                                                                                             >  1:430276/1‑135 (MQ=255)
                                                                                                                                                             ccctttccccGGGTTTTTTTGTTGCTCTAAAGGGGGGGTAATTGTGGGATGCAAAATATGGTGGATTGGTTGAGGGGGAGAAAATGTGGAGAAAGGGGGGGTAATTAATAAGGCCCAGGGCGTGGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGAt                                                                               <  1:829904/232‑1 (MQ=255)
                                                                                                                                                                       gggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                                                                                                                                                                                                                   >  2:4935/1‑88 (MQ=255)
                                                                                                                                                                       gggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAgg                                                                                                                                                                                                                                   <  1:4935/88‑1 (MQ=255)
                                                                                                                                                                       gggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa                                                                                                                                                                                                 <  2:916224/121‑1 (MQ=255)
                                                                                                                                                                       gggATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGa                                                                                                                                                                                                 >  1:916224/1‑121 (MQ=255)
                                                                                                                                                                       gggATCTTTGGTTGCTCTATAGGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGAAAATGGGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGa                                                                    <  2:970221/239‑1 (MQ=255)
                                                                                                                                                                                       ttataAGCGGGGTAATTGTGCGATGCAAAATATCGTTGATTTGTTGGGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGa                                                             <  2:759755/229‑1 (MQ=255)
                                                                                                                                                                                                 ggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGGGCAGAGCGTCGTTTTGAAAATAGGGGCTGGGTTCCATAAAAAAAGGCCCGCCCTAGGCGGGCCGGGGTTTCATTTGCGGGGCGTTTGATTGTAACTCACCCCATACGCGGGACAACTGTCAGGGAAGCAACCCAAACCGAAAAAAAccaacttccgtctct                                          >  1:759755/1‑226 (MQ=255)
                                                                                                                                                                                                 ggTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGCGATGCGTTTGATTTGAACCGACCCAATCAACGGGACAACTGTCAGGAAAGAAGCCAAACCCTGAAAAAACCACTTGGCCGCCgg                                         >  2:887233/1‑239 (MQ=255)
                                                                                                                                                                                                  gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGCGCTGCGTTAGAATTTAAAGACACCAATCAGCGGGACAAACCTCAGGGTAGCAACCAGACCGTTAAAAAACCACTTTCCCCCAggg                                        >  2:979529/1‑238 (MQ=255)
                                                                                                                                                                                                  gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTGGCGATGCGGTAGATTGTAACGACACCAATCAGCGGGACAACTGTCAGGATAGCAGCCAGACCGTAGaaaatgtctcttatacacaa                                          >  2:218256/1‑222 (MQ=255)
                                                                                                                                                                                                  gTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGAACGGTTAGATTGTAAAGAAACCAAACAGCGTGGCAAATGGCCGGGTAACAGCCAAACCCTAGAAAACCCATTTGCCCCCAggg                                        >  2:475811/1‑238 (MQ=255)
                                                                                                                                                                                                     aTTGTGGGACGCAAAATCTTGTTGATTGTTTGGGGGGAAGAAAATGGGGAGAAAGGGGCGGTAATTAATAAGGCGAAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACg                                     <  2:288296/239‑1 (MQ=255)
                                                                                                                                                                                                                cacaATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCCGGTTTTCAATTGGGGTGCGTTTGATTTTTACCAACCCCAACCAGCGGACAAACTGTCGGGTAAGAACCCAACCCTACAAAAACCAATTTCCCGCaagcacagggaattttca                          >  1:24730/1‑222 (MQ=255)
                                                                                                                                                                                                                             tttttggagggggCAACAAAGGGGGAGCAAGGGGCGGTAATTAATAAGGCCCAGGGCGTCGTTTTGAACATGGGGTCTGGTTTCCATAAAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATcgcg             <  2:88309/230‑1 (MQ=255)
                                                                                                                                                                                                                              tttGTTGAGGGGCAGAAAATGGGGAGCAAGAGGGAGTAATTAATAAGGCGCAGGGCGTCGTTTTGAAAATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgt            <  2:24730/239‑1 (MQ=255)
                                                                                                                                                                                                                                 tttgGGGGGAAGAAAAGGGGGAGAAAGGGGGGGGAATAAAGAAGGGGCAGAGCGTGGTTTGGAAAATGGGGTTTGGTTTCCATAAAAAACGGGCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt         <  1:926589/238‑1 (MQ=255)
                                                                                                                                                                                                                                 ggtgAGGGGAAGAAAATGGGGAGCAAGGGGCGGAAATTAAAAAGGCGAAGGGCGTCGTTTTGAACAAAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCgtggt         <  2:95819/237‑1 (MQ=255)
                                                                                                                                                                                                                                  ttgAGGGGCAGAAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGa                                                                                                                      <  2:138649/133‑1 (MQ=255)
                                                                                                                                                                                                                                      ggggCAAAAAAGGGGGAGAAAGAGGTGGTAATTAATAGGGGCCAGAGCGTCTTTTTGAACAGAGGTTTTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCAt     <  2:598813/238‑1 (MQ=255)
                                                                                                                                                                                                                                       gggCAGCAAATGGGGAGCAAGGGTCGGTAATTAAAAAGGCGCAGAGCGTGGTTTTGAACATAGGTTGTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATa    <  2:932181/238‑1 (MQ=255)
                                                                                                                                                                                                                                         gcagcaAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGGCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGGTTTACATTTGCGATGCGTTAGATTGTAAAGACACCAATCAGCGTGAAAACTGGCAGGATAGCAGCCAGACCCTAGAAAACCCCATTTGCCGCAGGGACGTGGGATTTCAAATCTTGCaagcccgggggccaacg  >  2:948592/1‑222 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                  
CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGAT‑‑AAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCTGGTTTCCATACAAAACGGCCCGCCATAGGCGGGCCGGATTTACATTGGCGATGCGTTAGATTGTAACGACACCAATCAGCGTGACAACTGTCAGGATAGCAGCCAGACCGTAGAAAACCCATTTGCCCGCAGGTACGTGGATTTTCAGATCGTGCATCGCGTGGTGCATACG  >  NC_000913/4378656‑4379124

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