Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC39.7% 111.5 / 14.7 59V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (35/0);  new base C (0/23);  total (35/24)
Fisher's exact test for biased strand distribution p-value = 1.14e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378653‑4378951
                                                                                                                                                       |                                                                                                                                                     
ccccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑a                                                                                                                                                       <  1:502632/150‑1 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accg                                                                                                                                                    >  2:46661/1‑151 (MQ=255)
  ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accg                                                                                                                                                    <  2:1217552/151‑1 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTATCCTGATAAATT‑accgcc                                                                                                                                                  >  2:1089994/1‑150 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accgcc                                                                                                                                                  >  2:817042/1‑150 (MQ=255)
     cGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accgc                                                                                                                                                   <  1:586273/149‑1 (MQ=255)
       aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accgccacc                                                                                                                                               >  1:1030014/1‑151 (MQ=255)
          aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑accgccaccgcc                                                                                                                                            >  2:88114/1‑151 (MQ=255)
                 tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCg                                                                                                                                     >  1:692532/1‑151 (MQ=255)
                 tCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCg                                                                                                                                     >  1:962855/1‑151 (MQ=255)
                   gCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCggg                                                                                                                                   >  1:1316273/1‑151 (MQ=255)
                        caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCtt                                                                                                                              >  1:327574/1‑151 (MQ=255)
                            cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                          >  1:784307/1‑151 (MQ=255)
                             aCGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTttgtt                                                                                                                         >  1:696139/1‑151 (MQ=255)
                                ccAGAAATTGTTTTGCTGCCCAACCCAACCAAAAACACGGTTTGGTGTGTGGGGGTTTGTTTTTTGGCGTTGCCCGCACACCCCGAGAAAAAGTAAA‑AAGTTTTCCCCTAATAATTC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                       <  2:505483/150‑1 (MQ=255)
                                      atttttttGCTGACCAACCCAAACTAAAACACCGTTTGTTGTGGGAGGGGTTGTTTTTTGGGGATGCCCGAAAACCCCCGGATAAAGAAAT‑AAGTTTTCCCCTAATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATaa                                                                                                                <  1:828995/147‑1 (MQ=255)
                                        tGTTTTGCTAACCGACCCAAACTAAAACAACGTTTGCTGTGTGACGGGTTGTTTTTTGGGGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTCCCTTGATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                               <  1:622012/150‑1 (MQ=255)
                                           tttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACAGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                           >  2:3701/1‑151 (MQ=255)
                                            ttccTGCCCAAACCTACCACAAACAACGTTTGTTGTGTGACGGCCTTTTTTTTGGCGATGCCCCCATCCCCCAAGTTAAAGTAAT‑AAGGTTTCCCCTAATAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                           <  1:645867/147‑1 (MQ=255)
                                            ttGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCgggg                                                                                                          >  2:724344/1‑151 (MQ=255)
                                              gCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                         >  1:1091927/1‑150 (MQ=255)
                                               cggcccGAACCTACCCCAAACACCTTTTGCGTTGTGACGGTTTTTTTTTTGGCGATGCCCGCTCCCCCCCGGGTAAAGTAAA‑AGGGTTTCCCTTAAAAAATT‑GCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                         <  2:88905/145‑1 (MQ=255)
                                                tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAAt                                                                                                      >  1:222430/1‑151 (MQ=255)
                                                     gAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtgtg                                                                                                  >  2:1210388/1‑150 (MQ=255)
                                                      accaaaccgcAAAGAAAGTTTGCTGGGAGGGGGGTTTCTTTGTGGTGATGCCCACATCGCCCGAGGTAAATAATT‑TTGTTTTCCCCTGAAAAATT‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCg                                                                                                <  2:154454/142‑1 (MQ=255)
                                                             cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                          >  2:983655/1‑150 (MQ=255)
                                                             cTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                          >  1:501997/1‑150 (MQ=255)
                                                                aaaCACCGTTGGCTGTGTGAGGTGTTTTTTTTGGGCGAGTCCCGCACAGCCCCAGAAAAAGTAAT‑AAGTTTTCCCTTAATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                         <  1:141668/148‑1 (MQ=255)
                                                                aaaCAACTTTTGCTGTGTGCGGGGCTGTTTTTGGGCGATGCCCGCATACCCCCGGAAAAAGTAAA‑AAGTTTTCCCCTAAAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                       <  1:257651/150‑1 (MQ=255)
                                                                aaaCAACTTTTGCTGTGTGAGGGGCTGTTTTTGGGGGATGTCCGAATCCCCCCAGGAAAAGAAAT‑AGGTTTTCCCCTAATAAATC‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcac                                                                                          <  1:718517/147‑1 (MQ=255)
                                                                aaaCAACGTTTGCTGTGGGGCGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGAAAT‑AAGTTTTCCCCTGATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                       <  1:929636/150‑1 (MQ=255)
                                                                   caacGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCg                                                                                   >  1:830926/1‑151 (MQ=255)
                                                                                gtgAGGGGCTGTTTTTTGGCGTTGCCCGCATACCCCCAGGAAAAGTAAT‑AAGTTTTCCCCTAATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgtt                                                                         <  1:331218/148‑1 (MQ=255)
                                                                                 tgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTGAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                     >  2:611802/1‑151 (MQ=255)
                                                                                   aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                   >  2:1272289/1‑151 (MQ=255)
                                                                                   aCGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                   >  2:422710/1‑151 (MQ=255)
                                                                                     ggttttttttttGGCGTTGCCCCCACAGCCCGGGAAAAAGAAAT‑AGTTTTTCCCTTAATAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                                 <  2:157569/144‑1 (MQ=255)
                                                                                     gggCTTTTTTTGGGCGATGCCCGCAAACCCCCAGTAAAATAAAT‑AAGTTTTCCCATAAAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgca                                                                 <  1:187716/151‑1 (MQ=255)
                                                                                     gggCTTTTTTTGGGCGATGCCCGCAAACCCCCAGGAAAAGTAAT‑AAGTTTTCCCTTAATAATTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgggg                                                                   <  1:1273651/149‑1 (MQ=255)
                                                                                     gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                  >  1:567970/1‑150 (MQ=255)
                                                                                        ctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcag                                                                >  2:242057/1‑149 (MQ=255)
                                                                                          tttttttGGTTATGCCCGCAT‑CCCCCGAGTAAAATAAA‑TAGGTTTTCCCTTAAAAATTCCCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagc                                                               <  2:479454/146‑1 (MQ=255)
                                                                                             ttttGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                          >  2:791942/1‑150 (MQ=255)
                                                                                              tttGGCGGTGTCCGCAGGGCCCAGGGTAAA‑TAATGAAGTTTTACCGGGAAAATTC‑CCCCCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtgt                                                         <  1:914380/150‑1 (MQ=255)
                                                                                                tGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                       >  2:365628/1‑150 (MQ=255)
                                                                                                    gATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCa                                                   >  2:311898/1‑150 (MQ=255)
                                                                                                         ccGCATACCCCCAGGAAAATTATA‑AGGGTTTCCCCTAAAAAATT‑CCCCCAACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAaga                                                <  2:1284152/148‑1 (MQ=255)
                                                                                                              tAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGta                                         >  1:152248/1‑150 (MQ=255)
                                                                                                              tAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                          >  2:186683/1‑149 (MQ=255)
                                                                                                                   ccAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaat                                   >  1:1207505/1‑151 (MQ=255)
                                                                                                                        tAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAgg                               >  1:380541/1‑150 (MQ=255)
                                                                                                                        tAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgc                              >  1:824263/1‑151 (MQ=255)
                                                                                                                         aaaaTAAT‑AAGTTTTACCTTAAAAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                             <  1:687384/151‑1 (MQ=255)
                                                                                                                              aat‑aaGTTTTCCCCTGAAAATTT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCaga                         <  1:724007/150‑1 (MQ=255)
                                                                                                                                    gttttcccttaaaaaaTT‑CCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcg                      <  1:144692/135‑1 (MQ=255)
                                                                                                                                    gttttCCCCTAATAAATT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtc                    <  1:369548/148‑1 (MQ=255)
                                                                                                                                    gttttCCCCTAATAAATT‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                          <  1:960180/142‑1 (MQ=255)
                                                                                                                                      tttCCCCTAAAAAATC‑CCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgtcg                   <  1:88075/149‑1 (MQ=255)
                                                                                                                                        tACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGtttt               >  1:1271973/1‑151 (MQ=255)
                                                                                                                                           ccTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGAAATTAATAAGGCGCAGAGCGTCGtttt               >  2:1129972/1‑148 (MQ=255)
                                                                                                                                                 aaaTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg      >  2:1245327/1‑151 (MQ=255)
                                                                                                                                                 aaaTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg      >  1:587537/1‑151 (MQ=255)
                                                                                                                                                 aaaTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAgg      >  1:975434/1‑151 (MQ=255)
                                                                                                                                                  aaTT‑CCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATa        <  2:300036/148‑1 (MQ=255)
                                                                                                                                                   aTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg     >  2:703793/1‑150 (MQ=255)
                                                                                                                                                   aTT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg     >  1:752345/1‑150 (MQ=255)
                                                                                                                                                     t‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCt  >  1:91401/1‑151 (MQ=255)
                                                                                                                                                       |                                                                                                                                                     
CCCCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAAT‑AAGGTTTACCCTGATAAATT‑ACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGTCT  >  NC_000913/4378653‑4378951

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