Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,378,8020AC51.4% 97.8 / 132.4 73V48G (GTA→GGA) ampCpenicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (32/3);  new base C (0/37);  total (32/41)
Fisher's exact test for biased strand distribution p-value = 2.29e-17
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.03e-13

CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGT  >  NC_000913/4378655‑4378949
                                                                                                                                                         |                                                                                                                                                     
ccTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGGGTGACGGGCTGCTTTTTGGCGATGTCGGCATAGCCCCAGGTAAAGTAATAAGGGTTACCCTGATAAATTaccg                                                                                                                                                   >  1:252665/1‑151 (MQ=255)
     aGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccacc                                                                                                                                             >  2:117808/1‑151 (MQ=255)
      gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccg                                                                                                                                            >  2:1088808/1‑151 (MQ=255)
      gCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccg                                                                                                                                            >  1:841912/1‑151 (MQ=255)
        aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATATGGTTTACCCTGATAAATTaccgccaccgcc                                                                                                                                          >  2:264173/1‑151 (MQ=255)
        aaTAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTaccgccaccgcc                                                                                                                                          >  2:928761/1‑151 (MQ=255)
         aTAGCTTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTAccgccaccgcca                                                                                                                                         <  2:300949/151‑1 (MQ=255)
              gTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATAccc                                                                                                                                    >  1:1066581/1‑151 (MQ=255)
                  ccccaacccccGCAAGTAATTTTTTTGCTGACCAAACCTAACTAAAACACCTTTTGCTGGGTGAGGGGTTGTTTTTTGGCGATGCCCCCATACCCCCAGAAAAAGTAAAAGGTTTTCCCTTAATAATTTCCCGCCACCGCCATACCCGGGa                                                                                                                                <  2:198855/141‑1 (MQ=255)
                      caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCtt                                                                                                                            >  1:920494/1‑151 (MQ=255)
                      caAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCtt                                                                                                                            >  1:85465/1‑151 (MQ=255)
                       aaGCCCGCCAGAAATTGTTTTGCTGCCCGACCCAAACTAAAACACCGTTTGTTGTGGGAGGGGCTGTTTTTGGGCGATGCCCCCATCCCCCCAGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCt                                                                                                                             <  1:544131/149‑1 (MQ=255)
                       aaGCCCGCAAGAAAATTTTTTGCTGCCCAAACCAAACCAAAACACCTTTTGTTGTGTGAGGGCTTGTTTTTGGGGGATGCCCGCATACCCCAGGGAAAAGAAATAAGTTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCttt                                                                                                                           <  1:715428/151‑1 (MQ=255)
                       aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCttt                                                                                                                           >  1:629317/1‑151 (MQ=255)
                       aaGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCttt                                                                                                                           >  2:542127/1‑151 (MQ=255)
                          cccGCAAGTAATTTTTTTTCTACCCAACCCAACCTAAAAAAGCTTTTGCTGTGTGGCGGGCTGTTTTTGGGCTTGCCCCCCATCCCCCCAGGAAAAATAATAGGTTTTCCCTTAAAAATTTCCCCCCACCGCCATACCCGGGATCTTttgt                                                                                                                        <  2:84965/149‑1 (MQ=255)
                          cACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTttgt                                                                                                                        >  2:930031/1‑151 (MQ=255)
                            cGCCAGTAATTGTTTCGCGGACCGACCCGAGCAAAAACACTTTTTGCTGTGTGGGGGGCTCTTTTTTGGCGTCGCCCGCATCCCCCCGGGAAAAATAATAAGGTTTTCCTTAAAAATTTCCCCCCTCCGCCATACCCGGGATCTTttgtt                                                                                                                       <  2:815479/150‑1 (MQ=255)
                            cGCCAGTAAGTGTTTTGCGTCCCAACCCAAACCAAAACAATTTTGGCTGTGTGAGGGGCTTTTTTTGGGCGATGCCCGCATGGCCCCGGTAAAAGTAAAAGGTTTTCCCTTAATAAATCCCCCCCACCGCCATACCCGGGATCTTttgttg                                                                                                                      <  2:312428/151‑1 (MQ=255)
                              caagaaattttttgGCTGCCCAACCCAAACTAAAACAACGTTTGTTGTGTGCGGGGCTTTTTTTTGGCAATGCCCGCATCCCCCAAGGAAAAGAAATAGGTTTTCCCTTAAAAAATCCCCCCCACCGCCATACCCGGGATCTTttgttg                                                                                                                      <  2:795294/139‑1 (MQ=255)
                                aGTAAGTGTTTTGCTACCCAAACTTATCAAAAACAACGTTTGTTTTGTGGCGGGCTTTTTTTTGGCGATCTCCACATACCCCCGGAAAAATTAAAAGGTTTCCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGc                                                                                                                     <  1:217504/148‑1 (MQ=255)
                                  tAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                >  1:1099730/1‑151 (MQ=255)
                                   aatttttttCCTGCCCAAACCCAACCAAAACAGCTTTTGCTTTGTGGGGGGTTGTTTTTGGGCATTGTCGCCAGCCCCCGAGGAAAAAAAATAAGTTTTCCAATAAAAATTACCCCCCCCCGCCATACCCGGGATCTTTTGTTGCTCta                                                                                                                 <  2:982997/144‑1 (MQ=255)
                                   aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                >  1:781989/1‑150 (MQ=255)
                                   aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                >  1:219118/1‑150 (MQ=255)
                                   aaaTGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCtat                                                                                                                >  2:557747/1‑150 (MQ=255)
                                   aaaTGTTTTGCTGACCAACCCTACCTAAAACAACGTTTGTTGTGTGACGGGCTTTTTTTTGGCGATGCCCGCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGctc                                                                                                                   <  1:718096/147‑1 (MQ=255)
                                    aaTTTTTCGCAGACCGAACTTAACAAAAACAACGTTTGTGGTGCGGTGGGCTTTTTGTGGGCGTTTCCCGCCTGGCCACGGGTAATGAAAAAGGTTTTGCCCTAAAAATTCTCTCCCACCGCCATACCCGGGATCTTTTGTTGCTCtata                                                                                                               <  2:276984/150‑1 (MQ=255)
                                      tGTTTTGCGGACGAAACCTAACAAAAACAACTTTGGCTGGGGGGCGGGGTGTTTTTGGGGGATCGCCGAATTGCCCCAGATAAAGAAATAGGGTTTCCCCTAAAAATTTCCCGCCTCCGCCATACCCGGGATCTTTTGTTGCTCTATAAg                                                                                                             <  1:418930/150‑1 (MQ=255)
                                        ttttgcgccccGAACCAAACTCAACCAAGGTTTGCTGTGTGGCGGGCTGTTTTTTGGCGTTCCCCGCATCCCCCCGGGAAATTTAATAATGTTTCCCCAGAAAAATCCCCGCCGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCg                                                                                                           <  2:470761/141‑1 (MQ=255)
                                          ttGCTGCCCGACCCTAACTCAAACAACTTTTGCTGTGTGACGGGCTTTTTTTGGGCGATGTCCGCATACCCCCAGGAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCggg                                                                                                         <  1:72868/150‑1 (MQ=255)
                                           tGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGt                                                                                                       >  1:1061592/1‑151 (MQ=255)
                                              tGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTaa                                                                                                     >  1:975954/1‑150 (MQ=255)
                                                  cGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATtg                                                                                                  >  1:68430/1‑149 (MQ=255)
                                                         aaCTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGc                                                                                          >  1:52483/1‑150 (MQ=255)
                                                          accccaaacatgttttGCTGTGTGACGGGTTCTTTTTTGGCATTTTCCGCATAGCCCAAGTTAAAGAAATAAGTTTTCCCCTAAAAAATTACCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATg                                                                                           <  1:792859/136‑1 (MQ=255)
                                                            ccaaaaatacGTTTGCTGTTTGCCGGTCTGTTTTTGGGCAATGCCCCCATACCCCAGATAAAATAAATAGTTTTTCCCCTTTAAAATCCCGCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGcaca                                                                                       <  2:267051/143‑1 (MQ=255)
                                                              cAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAat                                                                                     >  2:554700/1‑151 (MQ=255)
                                                               aaaCAACGTTTGTTGTGTGAGGGGCTGTTTTTGGGCGATGCCCGCATCCCCCAGGAAAAAGTAATAAGTTTTCCCCTGAAAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGAt                                                                                            <  2:1066425/143‑1 (MQ=255)
                                                                     tttttcttatgtgAGTGGCTGTTTTTTGGCAACGTCCGCCCGGCCCCAGAAAAAGAAATAGGGTTAACCCTATAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                               <  2:358179/141‑1 (MQ=255)
                                                                       tttGCTGTGTGAGGGGCTGTTTTTGGGCGATGCCCGCATACCCCCAGGAAAAGAAAAAAGTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                               <  1:981734/148‑1 (MQ=255)
                                                                       tttGCTGTGGGCGGGGCTGTTTTTGGGCGATGCCCGCATACCCCAAGAAAAAGTAATAATTTTTCCCCTAATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGtt                                                                               <  2:428947/148‑1 (MQ=255)
                                                                         tGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGAttt                                                                          >  1:517363/1‑151 (MQ=255)
                                                                              tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                      >  2:68502/1‑150 (MQ=255)
                                                                              tgtgACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                      >  2:461109/1‑150 (MQ=255)
                                                                               gtgAGGGGTTGCTTTTTGGGGTTGCCCGCATAGCCCCAGTTAAAGAAAAAGGGTTTCCCCTAATAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGa                                                                     <  1:1040541/150‑1 (MQ=255)
                                                                                tggcGGGTTGTTTTTTGGCGATGCCCGCTTAGCCCTGGGAAAATAAAGAAGTTTTCCCTTAAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATttgttg                                                                      <  1:1045841/145‑1 (MQ=255)
                                                                                  ggggtctgctTTTTGGAGTTGCCCGAATGCCCACAGTTAAAATAATATGTTTTCCCTGTAAAAATTCCCGCAGCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAgg                                                                   <  1:1063881/144‑1 (MQ=255)
                                                                                    gggCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgc                                                                >  1:540852/1‑150 (MQ=255)
                                                                                      gctgctTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGCAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTCGTTGAGGGgcag                                                              >  1:121503/1‑150 (MQ=255)
                                                                                           tttttGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAtg                                                        >  1:1099215/1‑151 (MQ=255)
                                                                                           ttttggggtaTGCCCGAATAGCCCCAGGAAAAATAATAAGTTTTCCCCTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAAt                                                         <  1:436207/141‑1 (MQ=255)
                                                                                           ttttgGGCGATGCCCCCATACCCCCAGGAAAAGAAATAAGTTTTCCCCTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGgcagca                                                            <  1:736983/147‑1 (MQ=255)
                                                                                                ggggATGCCCGCATACCCCAAGAAAAAGTAATAAGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGa                                                    <  2:187120/147‑1 (MQ=255)
                                                                                                ggCGATGTCCGCATACCCCCAGAAAAAGTAAAAGGTTTTCCCCTAAAAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTgg                                                     <  1:555881/149‑1 (MQ=255)
                                                                                                         cccaaacccccAGGAAAAATAAAAGGGTTTCCCCTAAAAAATTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCg                                          <  2:641742/144‑1 (MQ=255)
                                                                                                            agaccccaAGGAAAAGAAATGAGTTTTCCCCTAAAAAATCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGt                                        <  2:685126/146‑1 (MQ=255)
                                                                                                             tAGCCCAAGGAAAAGTAATAAGTTTTCCCTTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCgg                                         <  1:1043403/148‑1 (MQ=255)
                                                                                                             aaccccaAGGAAAAGTAAAAAGTTTTCCCTTAAAAAATCCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaa                                      <  2:28212/148‑1 (MQ=255)
                                                                                                                 cccAGGAAAAGTAATAAGGTTTACCCTGATAAATTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                  <  1:166723/151‑1 (MQ=255)
                                                                                                                 cccAGAAAAAGTAATATGTTTTCCCCTAAAAATTCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGtaattaa                                  <  2:643437/151‑1 (MQ=255)
                                                                                                                    aGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                >  2:652729/1‑150 (MQ=255)
                                                                                                                     ggTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaata                                >  1:242791/1‑149 (MQ=255)
                                                                                                                     ggTAAAGTAAAAGGTTTTACCCGAATAAATCCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                               <  1:217891/150‑1 (MQ=255)
                                                                                                                      gaaaaaGTAAAAAGTTTTCCCTTGAAAAATACCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATtaataa                               <  2:121069/147‑1 (MQ=255)
                                                                                                                         aaaGTAATAATCTTTCCCTTGAAAATTTCCCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                           <  2:1103860/151‑1 (MQ=255)
                                                                                                                         aaaGAAAAAAGTTTTACCTTGAAAATTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcg                           <  1:854924/151‑1 (MQ=255)
                                                                                                                          aaGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                          >  2:94969/1‑151 (MQ=255)
                                                                                                                            gTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                          >  1:1112914/1‑149 (MQ=255)
                                                                                                                              aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGgcgc                          >  1:183962/1‑147 (MQ=255)
                                                                                                                              aataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCag                        >  1:1036386/1‑149 (MQ=255)
                                                                                                                               ataaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCagag                      >  2:915492/1‑150 (MQ=255)
                                                                                                                                 aaGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGcgt                   >  1:981210/1‑151 (MQ=255)
                                                                                                                                    ttttCCCTTGAAAAATTCCGTCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGc                     <  1:1021795/145‑1 (MQ=255)
                                                                                                                                                 aTAAATTCCCGCCCCCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt       <  1:143717/150‑1 (MQ=255)
                                                                                                                                                  tAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGATCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACAt       >  2:1076975/1‑149 (MQ=255)
                                                                                                                                                   aaaTTCCCCCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAg     <  2:343461/150‑1 (MQ=255)
                                                                                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAggg   >  2:369580/1‑150 (MQ=255)
                                                                                                                                                     aTTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGt  >  2:793198/1‑151 (MQ=255)
                                                                                                                                                         |                                                                                                                                                     
CCTCGAGCAATAGCGTCGCCACCAAGCACGCCAGTAAATGTTTTGCTGACCGAACCTAACTCAAACAACGTTTGCTGTGTGACGGGCTGCTTTTTGGCGATGTCCGCATAGCCCCAGGTAAAGTAATAAGGTTTACCCTGATAAATTACCGCCACCGCCATACCCGGGATCTTTTGTTGCTCTATAAGCGGGGTAATTGTGCGATGCACAATATCGTTGATTTGTTGAGGGGCAGCAAATGTGGAGCAAGAGGCGGTAATTAATAAGGCGCAGAGCGTCGTTTTGAACATAGGGT  >  NC_000913/4378655‑4378949

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