Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,038	0	A	G	43.1%	79.5 / 13.6	58	intergenic (+209/+527)	ycdU/serX	putative inner membrane protein/tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (6/27);  new base G (25/0);  total (31/27)
Fisher's exact test for biased strand distribution p-value = 4.74e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGT  >  NC_000913/1096889‑1097183                                                                                                                                                      |                                                                                                                                                           gATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATCCCGGTGAGGGGTGGa                                                                                                                                                            <  1:91171/150‑1 (MQ=255)  aTTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAtt                                                                                                                                                          <  1:40064/151‑1 (MQ=255)   ttgttgACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGtt                                                                                                                                                          >  2:894362/1‑150 (MQ=255)    tgttgACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTg                                                                                                                                                         <  1:498850/150‑1 (MQ=255)         aCTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGGGGGGGGGAGTGAc                                                                                                                                                       >  1:1169630/1‑147 (MQ=255)          cTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcg                                                                                                                                                    <  2:73358/149‑1 (MQ=255)             aGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcg                                                                                                                                               <  1:1025921/151‑1 (MQ=255)               tAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACtcgcttcgc                                                                                                                                              <  1:806994/150‑1 (MQ=255)                aaCCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACtcgctttgc                                                                                                                                              >  2:924268/1‑149 (MQ=255)                aaCCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTtcgctc                                                                                                                                            <  2:193809/151‑1 (MQ=255)                 aCCAGACAGTTGATGTGGACTATCTCCCCTCGTCCTAACATACGTGGGCGGGGTCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTGTAAATAGGGCGGTGGGGGGGGGATTGACtcgcttcgct                                                                                                                                             <  1:356966/149‑1 (MQ=18)                   cAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTGACTCCCCTCCCCCGcc                                                                                                                                         >  1:1247738/1‑151 (MQ=21)                     gACAGTTGATGTGCACGATTTCCCCTCGCCATAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCt                                                                                                                                       <  2:723720/151‑1 (MQ=255)                       cAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTACTCGCTTTGGCCGccc                                                                                                                                        >  2:786661/1‑148 (MQ=255)                            gATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCGGCTCGCTCGCCCCCCggg                                                                                                                                  >  2:997689/1‑149 (MQ=18)                               gtgCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGc                                                                                                                              <  1:471444/150‑1 (MQ=255)                                   aCGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCGCTTCGCCCCCCCCGCCGGCAGc‑cccc                                                                                                                         >  2:971506/1‑148 (MQ=18)                                       tttCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTAATCCCTTCGCTCCCCCCGCGGGCCGc‑cccc                                                                                                                         >  2:640752/1‑144 (MQ=18)                                        ttCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACt                                                                                                                    <  2:1139969/150‑1 (MQ=255)                                         tCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTTGAAGCGCCTCGCCCCCCCCGCGGGCGGC‑CCCCTCAc                                                                                                                     >  1:171561/1‑148 (MQ=17)                                                      aCAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGGGGGGGGGGGTGACTCGCCTCCCCCGCCCTGCGGGCAGC‑CCCCGCCCTGCg                                                                                                                 >  2:1061076/1‑139 (MQ=14)                                                              gggCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑Tg                                                                                            <  1:1092731/150‑1 (MQ=38)                                                                gCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCt                                                                                         <  2:1111025/151‑1 (MQ=35)                                                                                  tgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑cc                                                                  <  1:704091/150‑1 (MQ=31)                                                                                    cAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTTAATCGGCTCGCTCTCCCTTATGGGAGG‑CCGGTGACTGGGGTT‑GGGGTCCGTGCCAC‑TTGGTGACGCCA‑G‑GTG‑T‑G‑GA‑ccccc                                                               >  1:307258/1‑151 (MQ=2)                                                                                     aaTATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑ccccc                                                               <  1:1101773/150‑1 (MQ=31)                                                                                      atacACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCCCGCCCTGCGGGGCAC‑CCCCTCCCCGCGGTT‑CCCGCCTGTCCCACTTGCCTGT‑GCCA‑GATTG‑T‑C‑AA‑ccccc                                                               >  2:697878/1‑150 (MQ=2)                                                                                       tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCTCCCCCCGCCGGCAGC‑CCGCCCACTGCGTTCTACGG‑CTGTCCCAC‑TCGCTGTCGCCA‑G‑GTT‑T‑C‑GG‑ccccccg                                                             >  2:1100597/1‑148 (MQ=11)                                                                                        ataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTc                                                           <  1:71029/151‑1 (MQ=31)                                                                                           cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCggg                                                        <  1:982654/151‑1 (MQ=31)                                                                                           cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCgg                                                         <  1:895347/150‑1 (MQ=31)                                                                                             aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCGCTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGc                                                      <  1:421945/151‑1 (MQ=14)                                                                                               aaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCCCTTCCCCCGCCCCCCGGGCAGG‑CCCCCCACTGCGTTC‑ACCGTCTGGCCAAC‑TGGGTGTCGCCA‑G‑GTG‑T‑C‑GA‑CCCCCCGTCGGGGCtt                                                    >  2:294512/1‑151 (MQ=11)                                                                                                  aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCCGCGGGCAGC‑CCCCCCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGGTGTCGCCC‑G‑TTT‑T‑C‑CA‑CCCCCCGTCGGGGGtt                                                    >  2:426395/1‑148 (MQ=11)                                                                                                  aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGAGGGGGGGGTTGTATCGCCTCGCTCGCCCTTCGGGCGAC‑CCCCCCACCGCGGTT‑ACCGTCTCTTCACC‑CGGCTCGTGCCA‑C‑ATG‑T‑C‑Gc‑ccccccgggggggggtttt                                                 >  1:330670/1‑137 (MQ=2)                                                                                                     gCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGAGGGGGGGGTTGACTCTCTTCCCTCGCCCCGGGGGCGAC‑CCGCCCCCCCGGGTC‑ACCGTCTTTTCAAC‑CGGCTGTTGCCACG‑TTT‑T‑C‑AA‑CCCCCgccgggggcgtccctc                                               >  2:1073084/1‑136 (MQ=2)                                                                                                       ccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACCGCGTTC‑ACCGGCTGTCCAAC‑CTGCTGTCGCCC‑G‑GTG‑T‑T‑GA‑CCCCCCGGCGGGGGTTCTCATccc                                            >  2:1141411/1‑151 (MQ=11)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCGCTTCGCCCGCCCTGCGGGCAGG‑CCCCCCACCGCGTTT‑CCCGTCTTTCCCAC‑TGGCTGTCCCCC‑G‑GTG‑T‑T‑GACCCCCCGGCCGGGGC‑TCTCCTtccccc                                         >  2:1278149/1‑150 (MQ=2)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATcccccc                                         <  1:140167/150‑1 (MQ=31)                                                                                                            cTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGGGGGGGGGGTTGACGCCCCTCCCTCGCCCTCCCGGGAGC‑CCCCCCCCCGCGGTT‑CCGGCCTTTTCCCC‑CGGGGGTCGCCC‑G‑CGGTT‑G‑GA‑CCCCCGG‑CGGGGCGTCTCCTccccccc                                        >  1:140252/1‑149 (MQ=0)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCCCCTCGCTCGCCCTCCGGGGAGC‑CCGCCCACAGCGTTC‑ACCGTCCGTCCACC‑TGGCCGGCGCCC‑G‑GTT‑T‑C‑CT‑CCCCCCGTCGGGGGCTTTCCTcccccccctg                                     >  2:388104/1‑147 (MQ=2)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGgt                                      <  2:464879/150‑1 (MQ=31)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGAATTAATCCCTTCGCTCGCCCCGCGGGGAGC‑CCCCTCACTTCGTTC‑ACCGTCTGTCCCAC‑TGGCTGTCGCCC‑G‑TTG‑T‑C‑GA‑CCCCCCGGCGGGGGTTCTCCTCCCCCCCgggtg                                   >  2:1182749/1‑151 (MQ=11)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGAATCGCTTCCCTCGCCCTTCTGGCGGG‑CCCCCCCCTCCGTTT‑AACGTCTTTCCAAC‑ATGCTGGTGCCC‑G‑GTG‑T‑CGGA‑CCCCCGG‑CGGGGGGTCTCATcccccccg                                       >  2:30812/1‑144 (MQ=0)                                                                                                                    cGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCa                               <  2:1019261/151‑1 (MQ=31)                                                                                                                    cGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCa                               <  1:976293/151‑1 (MQ=31)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTGACTCCCTTCCCTCCCCCTTCGGGGGGG‑CCCCTCCACGGGGTC‑TCGGTCTGTCCACC‑CGGGCGTCGCCC‑GATGG‑T‑C‑GA‑CCCCC‑GTCGGGGCGTTTCTTCCCCCCCgggtgt                                  >  1:95624/1‑147 (MQ=2)                                                                                                                     gAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCaa                              <  2:557119/151‑1 (MQ=31)                                                                                                                       gCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATa                            <  1:603628/151‑1 (MQ=31)                                                                                                                            tcttTAAATATGGCGGTGAGGGGGGGGGTGAATCGCTTCCCTCCCCCCCCGGGCAGC‑CCCCCCCCCGCGTTT‑ACCGTCTGTCCCAC‑CGGCTGTCGCCC‑G‑GTG‑T‑T‑GA‑CCCCCCGCCGGGGGTTCCCCTTCCCCCCgcgtggggcaaaaa                          >  1:909338/1‑139 (MQ=2)                                                                                                                             cttTAAATATGGCGGTGAGGGGGGGGTTGACTCGCCTCGCTCGCCCTGCGGGGCGC‑CCGCCCCCCCGGTTC‑ACGGGCTTTTCAAC‑CGGCCGGCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGCTCGGGGGGTCTCCTCCCCCCGgggtgggggaaata                          >  2:1243363/1‑138 (MQ=2)                                                                                                                              tttAAATATGGCGGTGAGGGGGGGGGTGAATCGCTTCGCTCGCCCTGCGGGGAGCGCCGC‑CACCGCGGTT‑ACCGGCTGTCCCAA‑TGGCGGGCGCCA‑G‑TTT‑TCC‑GC‑CCCCCGGCCGGGGCTT‑TCCTCCCCCCGGTGGGTGCAATACACGaa                      >  2:404180/1‑150 (MQ=2)                                                                                                                               ttAAATATGGCGGTGAGGGGGGGGGTGACTCGCTTCCCTCGCCCCGCGGGGAGC‑CCCCCCCCTGCGTTT‑ACCGTCTGTCCCAC‑TTGCTGTCGCCC‑G‑TTG‑T‑C‑GA‑CCCCCCGGCGGGGGGTCTCCTCCCCCCCGTGTGGGCCATACACGaaa                     >  2:690163/1‑150 (MQ=2)                                                                                                                                tAAATATGGCGGTGAGGGGTGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGTCTTCTCATCCCCCCGGTGTGTTCAATACACGaaaa                    <  2:1305473/150‑1 (MQ=16)                                                                                                                                tAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGaaaa                    <  2:371200/150‑1 (MQ=33)                                                                                                                                   atatGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAGTACACGgaaaaga                 >  1:1145107/1‑148 (MQ=16)                                                                                                                                     atGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGGCGCCA‑G‑TTG‑T‑C‑GA‑ACCCCGGTCGGGGCATCGTATCTTCCCGGCGTGGGGAGGACATGcaaaaaaac               >  1:516492/1‑149 (MQ=1)                                                                                                                                        gCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCg           <  1:584014/151‑1 (MQ=33)                                                                                                                                                    ggATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGt  >  1:1196750/1‑148 (MQ=33)                                                                                                                                                      |                                                                                                                                                           GATTGTTGACTGAGTAACCAGACAGTTGATGTGCACGATTTCCCCTCGCCCTAACAGACGTGGGCGGGGGCACCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGATTGACTCGCTTCGCTCGCCCTGCGGGCAGC‑CCGCTCACTGCGTTC‑ACGGTCTGTCCAAC‑TGGCTGTCGCCA‑G‑TTG‑T‑C‑GA‑CCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATACACGAAAAAAAAGCCCGTACTTTCGT  >  NC_000913/1096889‑1097183

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