Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,481,8040GA40.8% 109.2 / 31.1 71P125S (CCA→TCA) emrYputative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (36/6);  new base A (0/29);  total (36/35)
Fisher's exact test for biased strand distribution p-value = 2.37e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCT  >  NC_000913/2481654‑2481954
                                                                                                                                                        |                                                                                                                                                       
atatCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTgg                                                                                                                                                         >  2:163810/1‑151 (MQ=255)
  atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGaa                                                                                                                                                       <  2:1171454/151‑1 (MQ=255)
          ccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGc                                                                                                                                              <  2:319896/151‑1 (MQ=255)
            aGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGccc                                                                                                                                            >  1:941786/1‑151 (MQ=255)
            aGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGccc                                                                                                                                            >  2:832703/1‑151 (MQ=255)
            aGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGccc                                                                                                                                            >  1:136932/1‑151 (MQ=255)
              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCg                                                                                                                                          >  1:988845/1‑151 (MQ=255)
                aaaGTTATCACAAATATACCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAACCTCTGTGACAGTGGAATTAACGGCCCCGcc                                                                                                                                        <  1:636959/151‑1 (MQ=255)
                 aaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGcc                                                                                                                                        >  2:1255657/1‑150 (MQ=255)
                 aaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCa                                                                                                                                       <  1:1214219/151‑1 (MQ=255)
                    ttATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCAtt                                                                                                                                     >  2:856316/1‑150 (MQ=255)
                             atatAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGa                                                                                                                            >  1:13545/1‑150 (MQ=255)
                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAc                                                                                                                          >  1:207956/1‑151 (MQ=255)
                              tataACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAAc                                                                                                                          >  2:801564/1‑151 (MQ=255)
                                             tCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAatat                                                                                                           >  1:983567/1‑151 (MQ=255)
                                               ggCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATc                                                                                                          >  1:549545/1‑150 (MQ=255)
                                                 ccaaaatatGGGAGCAATAACCAGGTTAATTACCAATAATGCAAAACAAATTTTTTTTTTTTCGGGGGGATATTTCTTAAAAAAAAACCTCTGTGAAAGAGAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAat                                                                                                             <  2:824702/140‑1 (MQ=255)
                                                  ccACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                                       >  1:1268893/1‑150 (MQ=255)
                                                     caTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACAt                                                                                                   >  1:276499/1‑151 (MQ=255)
                                                      atatGGGAGAGATAATCCCGGTCATTGACAATAATGCAAGAGAATATGTTTTTTTTTCGGGGGAATAATTCTTAAAAAATAACCTCTGTGAAAGAAAAATCAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                                                     <  2:413114/148‑1 (MQ=255)
                                                       tatCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATc                                                                                                  >  1:304552/1‑150 (MQ=255)
                                                        atCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCa                                                                                                 >  1:1297456/1‑150 (MQ=255)
                                                            gAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAGTAATAAACTCTGTGACAGTGGAATTAACGGCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGAACATCCAGAtt                                                                                            >  2:80605/1‑151 (MQ=255)
                                                              gCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                          >  1:331004/1‑151 (MQ=255)
                                                              gCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                          >  1:1197640/1‑151 (MQ=255)
                                                                 aTAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                                       >  2:218722/1‑151 (MQ=255)
                                                                   aaTCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGCCTCTAAAAAATATCAGCACATCAAGATTGGTCGa                                                                                      >  2:788428/1‑150 (MQ=255)
                                                                    accgcGGTTATTAACCATAAGCCAAGAACAATTTTTTTTTTTTCGGGCGAATAATCCCTAAAAAAAAAACTCTGAGACGGTGAAATTAGCGCCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                                                    <  1:479767/147‑1 (MQ=255)
                                                                       aCGGTCATTGACCAAAATGCCAGAGCAATTTTTTTTTTTTCGGGGGGATAATTCTTAAATAAAAAACTCTGTGAAAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATa                                                                                    <  1:879253/148‑1 (MQ=255)
                                                                              ttGACCATAATGCAAGACAAATTTTTTTTTTTTCGGGCGGATAATTCTTAAAAAAAAACCCTGGTAAAAGTGAAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                                                             <  1:575067/148‑1 (MQ=255)
                                                                               tcaCCAAAATCCAAGAGAAATTTTTTTTTTTTCTGGGGAATATTTCCTAAATAAAAAACTCTGTAACGGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                                                          <  1:19823/148‑1 (MQ=255)
                                                                               tGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                         >  1:392259/1‑151 (MQ=255)
                                                                                gACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                         >  1:362433/1‑150 (MQ=255)
                                                                                 aCCATAATGCCAGAGAAAATGTTCTTTTTTCGGGCGGATATTTCCTAAAAAAAAAACTCGGTGAAAGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAAGCTaca                                                                             <  1:440596/145‑1 (MQ=255)
                                                                                  ccATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTaa                                                                       <  1:150283/150‑1 (MQ=255)
                                                                                  caaTAATCCAAGAGAAATTTTTTTTTTTCCTGGCGAAAATTTCTTAAAAAAAAACCTCTGTAAAAGTAAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                         <  2:867057/146‑1 (MQ=255)
                                                                                   cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                     >  1:971270/1‑151 (MQ=255)
                                                                                      aaTGCCAGAGAAATTTTTTTTTTTTCGGGCGAATAATTCTTTAATAAAAACCTCTGGGACAGGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAg                                                                  <  1:673635/151‑1 (MQ=255)
                                                                                       aTGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGGGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                                    >  2:846629/1‑148 (MQ=255)
                                                                                       aTGCCAAAGAAAATTTTTTTTTTTCGGGCGGATAATTCCTTAAAAAAAAACTCTGAAACAGTGAAATAAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGa                                                                 <  2:1074476/151‑1 (MQ=255)
                                                                                          ccAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGa                                                                 >  2:402161/1‑148 (MQ=255)
                                                                                            agagaatttGTTTTTTTTTCTGGCGATTATTTCCTTAATAAAACCCTCTGTGAAAGTGAATTTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                                                               <  2:1228098/140‑1 (MQ=255)
                                                                                             gataaaaTGTTTTTTTTTCGGGCGTATATTCCATTAATAAGACATTCTTTAAAGGTGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                                <  1:485389/142‑1 (MQ=255)
                                                                                             gagCAAGTTTTCTTTTTTTTGGCGGATATGTTCTAAATAAAAACCATAGAGCGAATGAAATCGCCGGCCCCGCCATTAACCCCTGAACGACTCGAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGc                                                             <  2:724644/149‑1 (MQ=255)
                                                                                                 aattttttttttttCGGGGGGATATTTCATAAAAAAAAACCTCTGTGAAAATAAAATAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaa                                                         <  2:760235/141‑1 (MQ=255)
                                                                                                 aattttttttttttCGGGCGAATATTTCTTTAAAAAAAACCTCTGTGAAGGGGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                       <  1:263482/143‑1 (MQ=255)
                                                                                                 aattttttttttttCGGGCGAATATTTCCAAAAAAAAAAACTTTGAGAAAGGGAAATAAGCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                           <  2:986516/139‑1 (MQ=255)
                                                                                                 aacgttttttttttCGGGAGAATATTTCCTAAATAATAAACTTTGGGCCATTGGAATAAAGGTCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                       <  1:1087466/143‑1 (MQ=255)
                                                                                                 aaaTTTTTTTTTTTCGGGCGAATAATTCTTTAAAAAAAACCTCTGTGACAGTGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaa                                                         <  1:584103/149‑1 (MQ=255)
                                                                                                 aaaTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaa                                                        >  2:1218445/1‑150 (MQ=255)
                                                                                                      ttCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAg                                                   >  2:612126/1‑150 (MQ=255)
                                                                                                            tttCTGGGGGATATTTCTTAAAAAAAAACCTCGGTGACAGTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTg                                             <  1:232478/150‑1 (MQ=255)
                                                                                                            tttCTGGCGAATAATTCCTAAAAAAAAACCTCTGTGAAAGTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTg                                             <  2:1106564/150‑1 (MQ=255)
                                                                                                               cTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaa                                           >  2:666503/1‑149 (MQ=255)
                                                                                                               cTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaa                                           >  2:219255/1‑149 (MQ=255)
                                                                                                                           tCCTTAAAAAAAAACCTCGGTAACAGTGAAATAAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaa                                      <  1:408290/143‑1 (MQ=255)
                                                                                                                              ttAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                           >  1:120766/1‑150 (MQ=255)
                                                                                                                              ttAAAAAAAATCTCGGTGAAAGAGAAATTAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaataat                               <  2:1006607/146‑1 (MQ=255)
                                                                                                                                aaaaataaACTCTGGGACAGGGAAATTAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCtt                            <  1:1123531/144‑1 (MQ=255)
                                                                                                                                  taataaACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                      >  2:212860/1‑151 (MQ=255)
                                                                                                                                   aaaaaaCTCTGTGATAGTAAAATTAAGGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                      <  2:647456/147‑1 (MQ=255)
                                                                                                                                   aaaaaaCTCTGGGACAGGGAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                      <  1:626758/147‑1 (MQ=255)
                                                                                                                                    aaaaacctctGTGAAAGTAAAATTACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                      <  2:1262446/144‑1 (MQ=255)
                                                                                                                                      aaaCTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt                   >  1:1387272/1‑150 (MQ=255)
                                                                                                                                      aaaCTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCt                   >  1:503646/1‑150 (MQ=255)
                                                                                                                                      aaaCTCTGTGACAGTAAAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTa                  <  1:1039477/151‑1 (MQ=255)
                                                                                                                                        aCTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAtt                >  2:843671/1‑151 (MQ=255)
                                                                                                                                          ctctGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTc               >  1:1231101/1‑151 (MQ=255)
                                                                                                                                             tctggaagagggAAATAACGCGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAt                 <  2:535652/138‑1 (MQ=255)
                                                                                                                                               gtgAGAGTGAATTTAGCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt            <  1:913694/150‑1 (MQ=255)
                                                                                                                                               gtgAAAGGGAAATTAACGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCt              <  2:13546/148‑1 (MQ=255)
                                                                                                                                                  aaaGTGAAATTACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt            <  2:1197612/145‑1 (MQ=255)
                                                                                                                                                      tGGAATTAACGGCCCCGCCATTAACCCCTGAACGTCTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGCGACTGAAAGTACAAATAATCTTAATTCGCTTATTGTTTGTGCCa      >  2:476524/1‑149 (MQ=255)
                                                                                                                                                        gAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCt  >  2:1024974/1‑151 (MQ=255)
                                                                                                                                                        gAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCt  >  2:103685/1‑151 (MQ=255)
                                                                                                                                                        |                                                                                                                                                       
ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAACGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACCT  >  NC_000913/2481654‑2481954

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