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

TATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAC  >  NC_000913/2481655‑2481952                                                                                                                                                         |                                                                                                                                                       tatCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTgg                                                                                                                                                        >  2:150214/1‑150 (MQ=255)    ccAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAAtt                                                                                                                                                    >  1:463980/1‑151 (MQ=255)     cAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTa                                                                                                                                                   <  2:1123096/151‑1 (MQ=255)       aCCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑c                                                                                                                                                >  1:242211/1‑151 (MQ=255)        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑Cg                                                                                                                                               <  1:1251920/151‑1 (MQ=255)        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑Cg                                                                                                                                               <  1:244159/151‑1 (MQ=255)        ccccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑Cg                                                                                                                                               <  1:988039/151‑1 (MQ=255)          ccAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGc                                                                                                                                             <  1:676291/151‑1 (MQ=255)              cTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGcccc                                                                                                                                          >  2:132792/1‑150 (MQ=255)                aaaGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCg                                                                                                                                         <  2:486957/149‑1 (MQ=255)                      aTCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAAc                                                                                                                                 >  2:814032/1‑151 (MQ=255)                                 aaCCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGAc                                                                                                                      >  2:866776/1‑151 (MQ=255)                                   cccccaaaTACCGGCCCACATTTCGGAGCGATAACCACGGTCATTGACCATAATGCCAGAGAAATTGTTTTTTTTTCGGGCGGATAATTCTTTAATAAAAAA‑CT‑C‑TGTGACAGGGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGAct                                                                                                                     <  1:629143/144‑1 (MQ=255)                                    cGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGActct                                                                                                                   >  1:87348/1‑151 (MQ=255)                                         aaaagcgcccaaaatatCGGGGAGATAACCACGGTCATTGCCAATAACCCAAAACAAAATGTTTTTTTTTCTGGCGAATAATTCTTAAAAAAAAAA‑CT‑C‑TGGAACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaa                                                                                                                <  1:77787/136‑1 (MQ=255)                                         aaTATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaaa                                                                                                               >  2:753466/1‑150 (MQ=255)                                         aaTATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaaa                                                                                                               >  2:246570/1‑150 (MQ=255)                                         aaTACCGCCCCACATAGCGGAGCAATAATCACGTTCATTGACCAAAATCCAAGAAAAAATTTTTTTTTTTCGGGGGGATAATTCTTAAAAAAAAAC‑CT‑C‑TGTGAAAGTAAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTaaaaa                                                                                                              <  1:1122870/151‑1 (MQ=255)                                                   cacaTATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcac                                                                                                    >  2:884398/1‑151 (MQ=255)                                                      atatCGGAGCTTTAATCGCGTTTATTGACAATAATGCCAGAAAAAATTTTTTTTTTTCTGGGGGATATTTCCTTAATAAAAAA‑CT‑C‑GGTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGcaca                                                                                                   <  2:319186/149‑1 (MQ=255)                                                       tatCGGGGCAATAATCACGGTCATTGCCCAAAAGGCCAGAAAAAATTTTTTTTTTTCGGGCGAATATTTCTTTAAAAATAAA‑CT‑T‑TGGGAAAGTGAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATc                                                                                                 <  2:826445/150‑1 (MQ=255)                                                           gggcggaaaaTCTCGGTTATTGACCATTATCACAAATCAAGTGTTTTTTTTTCAGGCGGAAAATCCCTTAAAAATAAA‑TT‑T‑TGGGCAAGGGTAATCAG‑GGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGa                                                                                             <  2:276248/142‑1 (MQ=255)                                                           ggaggaaaaaTCACGGTAATTAACCATAATGCCAGAGCAAATTTTCTTTTTGCGGGCGAATAACCCCTAAAAAAAAAA‑CT‑C‑TGAAAAAGAGAAATAAA‑CGGCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGAt                                                                                            <  1:522902/143‑1 (MQ=255)                                                           ggAGCGATAACCACGGTCATGGACATAAAGGCCAGAGCAATTTTTTTTTTTTCGGGCGTATAATTCTTAAAAAAAAAA‑CT‑C‑GGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCaa                                                                                               <  2:985374/148‑1 (MQ=255)                                                                  tAATCCCGGTCATTACCAATAATGCCAGAGAAATTTTTTTTTTTTCGGGCGAATAATTCCTAAAAAAAAAC‑CT‑C‑TGTGACAGTGAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTgg                                                                                         <  2:518973/147‑1 (MQ=255)                                                                  tAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCg                                                                                      >  1:280702/1‑150 (MQ=255)                                                                   aaTCGCGTTCTTTGACAATAAGCCCAGAAAAAATTTTCTTTTTTCGGGCGTATATTCCCTAAAAAAAAAC‑CC‑C‑TGTGACAGTAAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                    <  1:726159/151‑1 (MQ=255)                                                                     tCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGAt                                                                                    >  2:968590/1‑149 (MQ=255)                                                                       aCGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGc                                                                                >  1:653708/1‑151 (MQ=255)                                                                         ggTCATTGACAATAATGCAAGAGCAAATTTTTTTTTTTCTGGCGGATAATTCTTAAAAAAAAAC‑CT‑C‑TGTGACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCt                                                                               <  2:64088/150‑1 (MQ=255)                                                                             aTTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTacac                                                                           >  2:975579/1‑150 (MQ=255)                                                                              ttGACCAAAAGGCCAGAGAAATTTTTTTTTTTTCGGGCGAATATTTCCTAAAAAATAAC‑CT‑C‑TGTGACAGTGAAATAAC‑GGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTaca                                                                            <  2:575081/148‑1 (MQ=255)                                                                                gACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACAtt                                                                        >  2:1220728/1‑150 (MQ=255)                                                                                   cATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATg                                                                    >  1:535092/1‑151 (MQ=255)                                                                                   cATAATGCCAGAGAAAATTTTTTTTTTTTGGGGGAATAATTCTTTAATAAAAAA‑CT‑C‑TGGGAAAGTGAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTa                                                                       <  2:629003/148‑1 (MQ=255)                                                                                      aacccgagagAAAATTTTTTTTTTTCGGGGGAAAATTTCCTAAAAAATAAC‑TT‑C‑GGTGCCAGTGAATTTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAAt                                                                     <  2:1053859/141‑1 (MQ=255)                                                                                        ttcaagagAAATTTTTTTTTTTTCGGGCGGTTATTTCCTAAATAAAAAA‑CC‑C‑GCGAAAAGTGAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGaa                                                                  <  2:1196711/144‑1 (MQ=255)                                                                                          ccAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGAc                                                               >  2:1137265/1‑149 (MQ=255)                                                                                          ccAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATCAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACa                                                              <  2:937792/150‑1 (MQ=255)                                                                                            agagAAATTGTTTTTTTTTCTGGCGAATAATTCTTAAAAAATAAT‑CT‑C‑TGTGACAGGGAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCg                                                           <  2:1215824/151‑1 (MQ=255)                                                                                              agaaaaTTTTTTTTTTTCTGGCGGATAATTCCTTAAAAATAAC‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑aa                                                        <  2:559509/148‑1 (MQ=255)                                                                                              agCAATTTTTTTTTTTTCGGGCGAAT‑ATTCCTTAAAAAAAAACCT‑C‑TGAGAAAGTAAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑a                                                         <  1:711157/150‑1 (MQ=255)                                                                                                   ttttttttttttCGGGGGGATAA‑TCCTTAAAAAAAAA‑CC‑CTTCGGAGAGGGAAATAAG‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑aaaaaa                                                    <  1:893693/144‑1 (MQ=255)                                                                                                        cTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTg                                               >  2:848451/1‑150 (MQ=255)                                                                                                         ttttttCTGGCGAATATTTCTTAAAAAAAAAA‑CT‑C‑TGTGACAGTGAAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGa                                              <  2:735094/150‑1 (MQ=255)                                                                                                         ttttttCGGGCGAAAAATTCTTAAAAAAAAAC‑CT‑C‑GGGAACAGTGAAATAAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTg                                               <  2:1014390/149‑1 (MQ=255)                                                                                                                 ggCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGt                                      >  2:565051/1‑150 (MQ=255)                                                                                                                 ggCGAATATTTCTTAAAAAAAAAC‑CT‑C‑TGTGACAGTGAAATAAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTa                                     <  1:360191/151‑1 (MQ=255)                                                                                                                 ggCGAATAATT‑CTTAAAAAAAAA‑CC‑C‑GGGGAAAAGTGAATTAACCGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAATAGTGACTGaaa                                        <  1:966946/148‑1 (MQ=255)                                                                                                                    ggATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaaa                                   >  2:1224181/1‑150 (MQ=255)                                                                                                                     gaaaaTTCCTAAATAATAAA‑CC‑C‑TGTGAGAGTGAAAATAG‑CGCCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaaaa                                  <  1:462481/147‑1 (MQ=255)                                                                                                                      aTAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAaat                                >  2:643630/1‑151 (MQ=255)                                                                                                                      aTAATTCCTAAATAAAAAC‑CC‑C‑TGTAACAGGGAAATTAC‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTaaa                                   <  2:915367/148‑1 (MQ=255)                                                                                                                            tcTTAAAAAAAAA‑CTAC‑TGTAAAAGTGAAATTAG‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATc                            <  2:900602/149‑1 (MQ=255)                                                                                                                             cTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCtt                          >  1:616360/1‑150 (MQ=255)                                                                                                                              ttAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                         >  1:667152/1‑150 (MQ=255)                                                                                                                              ttAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTa                         >  2:1020890/1‑150 (MQ=255)                                                                                                                                 ataataaA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                     >  1:571205/1‑151 (MQ=255)                                                                                                                                 ataataaA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTc                     >  2:1222916/1‑151 (MQ=255)                                                                                                                                 aaaaaaaa‑CT‑C‑TGGGAAAGGGAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTaa                        <  2:87351/143‑1 (MQ=255)                                                                                                                                   aataaA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                   >  2:394988/1‑151 (MQ=255)                                                                                                                                   aaaaaacct‑c‑tGGGACAGTAAAATAAC‑CCGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                    <  1:647787/144‑1 (MQ=255)                                                                                                                                   aaaaaaaaacgttatgACAGGGAAATAAA‑CGCCACCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAAtt                      <  1:1221637/137‑1 (MQ=255)                                                                                                                                   aaaaaa‑CT‑C‑TGTGACAGTGAAATAAC‑GGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCGAAAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCg                    <  2:336766/148‑1 (MQ=255)                                                                                                                                    ataaC‑CT‑C‑TGTAACGGTGGAATCAC‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGc                   <  1:1227709/150‑1 (MQ=255)                                                                                                                                      aaa‑CT‑C‑TGTAACAGGAAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTa                <  1:487354/151‑1 (MQ=255)                                                                                                                                               tgtgACAGTGAAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTc             <  1:64554/148‑1 (MQ=255)                                                                                                                                                gtgACAGTGAAATAAA‑CGCCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCt            <  2:1259433/148‑1 (MQ=255)                                                                                                                                                      gTGAAATCAA‑CGGCGCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCa    <  2:433861/150‑1 (MQ=255)                                                                                                                                                       tGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAc  >  1:1176264/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       TATCCAACCCCAGCTAAAGTTATCACAAATATAACCGCCCAATATCGGCCCACATATCGGAGCGATAATCACGGTCATTGACCATAATGCCAGAGCAAATGTTCTTTTTTCTGGCGGATAATTCCTTAATAATAAA‑CT‑C‑TGTGACAGTGGAATTAA‑CGGCCCCGCCATTAACCCCTGAACGACTCTAAAAAATATCAGCACATCAAGATTGGTCGATAGGCTACACATTAATGAAGACAGCG‑AAAAAAAAGTGACTGAAAGTAAAAATAATCTTAATTCGCCTATTCTTTGTGCCAAC  >  NC_000913/2481655‑2481952

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