Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,481,804	0	G	A	41.9%	80.1 / 27.2	62	P125S (CCA→TCA)	emrY	putative multidrug efflux system
Reads supporting (aligned to +/- strand):  ref base G (26/10);  new base A (0/26);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.30e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACC  >  NC_000913/2481654‑2481953                                                                                                                                                       |                                                                                                                                                       atatCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTgg                                                                                                                                                        <  2:39384/151‑1 (MQ=255)   atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                                                       <  1:469865/150‑1 (MQ=255)   atCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGa                                                                                                                                                       <  1:451264/150‑1 (MQ=255)     ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAGACTCTGTGACAGTGGAAtt                                                                                                                                                    >  1:301155/1‑151 (MQ=255)     ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAAtt                                                                                                                                                    <  1:387368/151‑1 (MQ=255)               cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCg                                                                                                                                         >  1:861781/1‑151 (MQ=255)               cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCg                                                                                                                                         >  1:217293/1‑151 (MQ=255)                  aaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCa                                                                                                                                      <  1:864892/151‑1 (MQ=255)                          aaaaaTAAACCCCCCAAATATCGCCCCACATATGGGAGCGATAATCACGGTCATGGCCCATAATCCAAGAGAAATTGTTTTTTTTTCGGGCGGAAAATTCTTAAAAAATAACCTCTGTGACAGTGGAATAAA‑CGGCCCCGCCATTAAcccc                                                                                                                              <  1:629950/149‑1 (MQ=255)                           caAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCt                                                                                                                             >  2:65923/1‑151 (MQ=255)                                                   ccACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                                      >  1:989481/1‑150 (MQ=255)                                                       atatggggcagtaaaTCAGGGTTTTTGACCAAATCGCAAGAGAAATTTTTTTTTTTTCGGGCGGATAATTCCTAAAAAAAAAACTTGGTGACAGAAAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGc                                                                                                      <  2:140417/133‑1 (MQ=255)                                                       atagCGGAGAAATAATCCCGTTCATTAACCAAAAGCCCAGAGAAAATGTTTTTTTTTCTGGGGAATAATTCCTAAAAAAAAACCTCTTTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcaca                                                                                                   <  1:748318/149‑1 (MQ=255)                                                        tatGGGAGCGATAATCCCGGTCATAGACCATAATGCCAGAAAAAATGTTTTTTTTTCTGGGGAATAATCCCTAAAAAAAAAACTCTGTGAAAGAAAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGca                                                                                                     <  2:696521/146‑1 (MQ=255)                                                        tatCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATc                                                                                                 >  1:382463/1‑150 (MQ=255)                                                               gCGATAATCACGGTCATGGACAATAATGCCAGAAAAAATGTTTTTTTTTCGGGCGGATAATTCCTAAAAAAAAAACTCTGTGACAGTGGAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                         <  1:336324/151‑1 (MQ=255)                                                                 aaTAACCCCGATCATTGACCAAAAGGCCAGACAAATTGTTCTTTTTTCGGGCGAATAATCCTTAAAAAAAAAACACTGGGACAGAAAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                         <  2:125787/148‑1 (MQ=255)                                                                    aaacACGGTCATTGACCATAATGCCAGAACAAATGTTCTTTTTTCGGGCGGATAATTCCTTAAAAAAAAACTCGGTGACAGGGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                                      <  1:417925/146‑1 (MQ=255)                                                                    aaTCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                    >  1:256687/1‑151 (MQ=255)                                                                     aTCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                    >  2:498775/1‑150 (MQ=255)                                                                       cACGGTCATTGACCATAATGCAAGAGAAAATTTTTTTTTTTCTGGCGAATAATTCCTAAAAAATAAACTCTGTAAAAGTGAAATTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAgg                                                                                 <  2:470032/151‑1 (MQ=255)                                                                             cATTGACCATAATGCCAGAGAAAATGTTTTTTTTTCTGGCGGATATTTCCTTAAAAAAAACCTCTGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                                                            <  1:739252/150‑1 (MQ=255)                                                                                tGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                        >  1:719369/1‑151 (MQ=255)                                                                                tGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                                                         >  2:135571/1‑150 (MQ=255)                                                                                  aCCATAATGCCAGAGAAAATTTTCTTTTTTCGGGCGGATAATTCCTTAAAAAAAAACTCTGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                        <  1:619014/149‑1 (MQ=255)                                                                                   ccATAATGCCAGAGAAAATGTTTTTTTTTCGGGCGAATATTTCCTTAAAAAAAAACTCTGTGACAGTGAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                        <  1:736956/148‑1 (MQ=255)                                                                                   ccAAAATGCCAGACAAATTTTTTTTTTTTCTGGCGAATAATCCCTTAAAAAAAACCTCTGTGACAGT‑GAATTAACCGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAt                                                                         <  1:149788/147‑1 (MQ=255)                                                                                    cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                    >  2:914166/1‑151 (MQ=255)                                                                                     aTAATGCCAAAGAAATTGTTTTTTTTTCTGGCGGATAATTCCTTAAAAATAATCTGTGAGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                                                     <  1:632727/149‑1 (MQ=255)                                                                                      tAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGa                                                                   >  2:440630/1‑150 (MQ=255)                                                                                            cAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGc                                                            >  2:643446/1‑151 (MQ=255)                                                                                            cAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGc                                                            >  2:900709/1‑151 (MQ=255)                                                                                             agagCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                                                           >  2:958783/1‑151 (MQ=255)                                                                                             agagAAATTTTTTTTTTTTCGGGAGGATAATTCTTAAAAAAAAAACTCGGTGACAGTAAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                               <  2:721482/147‑1 (MQ=255)                                                                                               agCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                          >  1:855840/1‑150 (MQ=255)                                                                                                 aaattttttttttttCTGGCGAATATTTCCTAAAAAAAAAACTCTGTGAAAGGGAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaa                                                        <  1:534553/141‑1 (MQ=255)                                                                                                  aaaTTTTTTTTTTTCGGGCGGATATTTCCTAAAAAAAAAACTCTGTGACAGTAAAATAAG‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGa                                                          <  2:588532/147‑1 (MQ=255)                                                                                                   aaTGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaa                                                      >  1:658734/1‑150 (MQ=255)                                                                                                     tttttttttttCTGGCGAATATTTCCTAAAAAATAAACCCTGTGAAAGAGAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGaaaaaaa                                                    <  2:749958/145‑1 (MQ=255)                                                                                                       tttttttttTGGGCGAATAATTCCTTAAAAAAAACCTCTGTGACAGGGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAg                                                  <  2:603300/147‑1 (MQ=255)                                                                                                         cTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTg                                                >  1:310625/1‑150 (MQ=255)                                                                                                          ttttttCTGGCGGATAATTCCTAAAAAAAAAACTCGGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGa                                               <  2:452709/150‑1 (MQ=255)                                                                                                             tttCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTg                                            >  1:891286/1‑150 (MQ=255)                                                                                                               tCTGGGGGATATTCCCTTAAAAATAAACTCTGTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaa                                          <  1:283367/150‑1 (MQ=255)                                                                                                                cTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGaa                                          >  1:639660/1‑149 (MQ=255)                                                                                                                  ggCGGTTAATTCTTTAAAAAAAACCTCTGGGACAGAGAAATTAC‑CGGCTCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa                                      <  1:710515/151‑1 (MQ=255)                                                                                                                   gCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTa                                      >  1:702943/1‑150 (MQ=255)                                                                                                                       aTAATTCCTAAAAAAAAACCTCTGAAAAAGGGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTaaa                                    <  2:980448/148‑1 (MQ=255)                                                                                                                         aaTTCCTAAAAAAAAACCTCTGGAACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaat                                 <  1:793418/149‑1 (MQ=255)                                                                                                                           ttCCTAAATAAAAACCTCTGTAACAGGGAATTTAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAaataa                               <  1:43411/149‑1 (MQ=255)                                                                                                                               ttAATAATAAACTGTGTGAAAGTGGAATTAA‑CGGCGCCGCCATTAAATAATGAACGACTCTAAAAAATATAAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                          >  2:89965/1‑150 (MQ=255)                                                                                                                               ttAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                          >  1:571531/1‑150 (MQ=255)                                                                                                                               ttAAAAAAAAACTCTGTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                          <  1:855413/150‑1 (MQ=255)                                                                                                                                 aaaaaaaacctctGTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAAtt                       <  1:490989/142‑1 (MQ=255)                                                                                                                                  ataataaACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGATACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                      >  1:976846/1‑151 (MQ=255)                                                                                                                                  aaaaaaaaCTCTGTGACAGTGGAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                      <  1:627253/146‑1 (MQ=255)                                                                                                                                    aaaaacctctGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                     <  1:442856/144‑1 (MQ=255)                                                                                                                                        aaCTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAt                >  1:333851/1‑151 (MQ=255)                                                                                                                                          ctctGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTa                 <  2:374790/148‑1 (MQ=255)                                                                                                                                          ctctGGGACAGTAAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTAtt               <  2:65440/150‑1 (MQ=255)                                                                                                                                              gtgAAAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt           <  1:138286/150‑1 (MQ=255)                                                                                                                                              gggACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCttt           <  1:517549/148‑1 (MQ=255)                                                                                                                                                  cAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGcc      >  2:677089/1‑151 (MQ=255)                                                                                                                                                      ggAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAcc  >  2:75404/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       ATATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAACTCTGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAACC  >  NC_000913/2481654‑2481953

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: ‑