Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,622,250	0	A	G	44.4%	87.8 / 25.9	72	E161G (GAG→GGG)	yjjJ	putative protein kinase
Reads supporting (aligned to +/- strand):  ref base A (14/26);  new base G (32/0);  total (46/26)
Fisher's exact test for biased strand distribution p-value = 1.14e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCT  >  NC_000913/4622102‑4622399                                                                                                                                                      |                                                                                                                                                       ccGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCggggagg                                                                                                                                                      <  2:974683/151‑1 (MQ=255) ccGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCggggagg                                                                                                                                                      <  1:904481/151‑1 (MQ=255) ccGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCggggag                                                                                                                                                       <  1:35426/150‑1 (MQ=255)  cGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCggggaggg                                                                                                                                                     <  2:266147/151‑1 (MQ=255)  cGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGGCgggggggg                                                                                                                                                     >  1:1071315/1‑151 (MQ=255)   gCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCggggagggg                                                                                                                                                    <  2:1060823/151‑1 (MQ=255)       ggCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGG‑TTGgggggggggggg                                                                                                                                                    >  1:1127348/1‑147 (MQ=255)        gCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTG‑TTGGTCGGGGGGGGGaaat                                                                                                                                                >  1:590930/1‑148 (MQ=255)         cTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTCGGGGGGGGGAATTAt                                                                                                                                              >  2:452998/1‑151 (MQ=255)               aGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGCGGGGGGGGGGATTATCAGC‑g                                                                                                                                        >  1:486553/1‑150 (MQ=255)                gggCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGCGGGGGGGGGGATTATCAGC‑g                                                                                                                                        >  1:1153740/1‑149 (MQ=255)                    gcgcgTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑CTGGTCGGGGGGGGGGATTATCAGc                                                                                                                                          >  1:1303028/1‑144 (MQ=255)                        gTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTG‑GTTGTCGGGGGGGGGGATTATccacggagg                                                                                                                                     >  1:57100/1‑137 (MQ=255)                        gTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACt                                                                                                                              <  2:504095/151‑1 (MQ=255)                         tGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGCTTG‑GTGGGCGGGGGGGGGAATTATGAGG‑Gaccggtttct                                                                                                                              >  2:994898/1‑141 (MQ=255)                               aGGAAGTTAGCCGCGCAACTGAGTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGG‑GTGGTTGGGGGGGGGAATTATCCGC‑GATGGGTTACTggccaa                                                                                                                        >  2:1021039/1‑145 (MQ=255)                                    gTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCt                                                                                                                  <  1:719943/151‑1 (MQ=255)                                       aGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCg                                                                                                               <  1:801899/151‑1 (MQ=255)                                          cgcgCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGaa                                                                                                            <  1:1212287/151‑1 (MQ=255)                                           gcgcAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGaaa                                                                                                           <  2:851999/151‑1 (MQ=255)                                              cAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTc                                                                                                        <  2:620682/151‑1 (MQ=255)                                                  tGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGCGGGGGGGGGAATTATCAGC‑GATGGATTACTGCCCAACACCCTGCGGAAATTCctct                                                                                                    >  2:191875/1‑151 (MQ=255)                                                    aaTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTgg                                                                                                  <  1:874390/151‑1 (MQ=255)                                                    aaTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTCGGGGGGGGGGATTATCAGC‑GATGGATTACTGCACAACAACCCGCGGGAAATCCTCTgg                                                                                                  >  1:287162/1‑151 (MQ=255)                                                    aaTCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGGGGGGGGGGGAATTATCAGC‑GATGGGTTAATGCACAACACCCTGCGGAAATTCCTCTgg                                                                                                  >  2:334717/1‑151 (MQ=255)                                                      tCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGGTGG‑GTGGTGGGGGGGGGGGATTATCACC‑GGTGGGTTATTACCCCACCCCCCGCGGGAAtt                                                                                                         >  1:983639/1‑142 (MQ=255)                                                       ctgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTGGGGGGGGGGAATTATCAGC‑GAAGGAGTACTGCACAACACCCTGCGGGAAATCCTCTGGATc                                                                                               >  1:1248943/1‑151 (MQ=255)                                                        tgactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTCGGGGGGGGGGATTATTAGC‑GATGGATTACTGGACAACACCCCGCGGGAAATCCTCTgg                                                                                                  >  1:650360/1‑147 (MQ=255)                                                         gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGa                                                                                                 <  1:257104/147‑1 (MQ=255)                                                         gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATc                                                                                               <  2:685346/149‑1 (MQ=255)                                                         gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCaa                                                                                             <  1:351113/151‑1 (MQ=255)                                                         gactgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGGCGGGGAGGGGAATTATTCGC‑GATGGATTTCTTCACAACACCCTGCGGAAATTCCTCTGGATCa                                                                                              >  2:911486/1‑150 (MQ=255)                                                          actgaTGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGGCGGGGAGGGGGATTATACTC‑GATGGATTACTGCGCAACACACTCCCGGAAATCCTCTTGATCaa                                                                                             >  1:601829/1‑150 (MQ=255)                                                            tgatgaTATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTCGGGGGGGGGGATTATCAGC‑GATGGGTTACTGCACAACACCCTGCGGAAAATCCTCTGGGTCaaaa                                                                                           >  2:1069476/1‑150 (MQ=255)                                                             gatgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGAATTATATGCGGAGGGATT‑CTCCCCAAAACCCCGCGGCAATTCCTCTGGATTaaaaa                                                                                          >  2:772732/1‑150 (MQ=255)                                                              atgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTCGGGGGGGGGGATTATCAGC‑GATGGGTTACTGCACAACACCCCGCGGGAATTCCTCTGGGTCaaaaa                                                                                          >  2:429476/1‑149 (MQ=255)                                                               tgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTGGGGGGGGGGGATTATCAGC‑GATGGATTTATGCACAACACCCCGCGGGAATTCCTCTTGATCaaaaa                                                                                          >  2:839212/1‑148 (MQ=255)                                                               tgatATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGCGGGGAGGGGGATTTTCAGC‑GATGGATTACTTCACCAAAACCCGCGGGAAATTCTCTTGGTCaaaaaat                                                                                        >  1:343548/1‑148 (MQ=255)                                                                  tataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGGATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCAcc                                                                                    >  2:1150027/1‑151 (MQ=255)                                                                   ataCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGCTTG‑TTGGTTGGGGAGGGGGAATATCAGC‑GAAGGGTTATTTCACAACACCCAGCGGAAATTACTCTGGACAAAAAAATCAacc                                                                                   >  1:1189002/1‑151 (MQ=255)                                                                     aCGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGAATTATCAGC‑GATGGATTAATGGACAACACCCTCCGGGAATTCCTCTGGGTTCAAAAATCAccc                                                                                   >  1:463393/1‑149 (MQ=255)                                                                      cGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGTGGGGGGGGGGAAATTTTAGG‑GAAGGGATATTGGCCAAAACCCCGCGGGAAATAATCTGGGTGAAAAAATCAccc                                                                                   >  1:1342187/1‑148 (MQ=255)                                                                      cGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑GTGGGTGGGGGGGGGGATTATCAGC‑GATGGTATAATGCACAAAACCCTGCGGGAATTCCTCTTGGTCAAAAAATCACCCAtt                                                                                >  2:565154/1‑151 (MQ=255)                                                                       gTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGATGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGGATTATTAGC‑GATGGATTACTGCACAACACCCCGCGGAAAATCCTCTGGATCAAAAAATCACCCAtt                                                                                >  1:1212859/1‑150 (MQ=255)                                                                           ctGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTAc                                                                              <  2:108960/148‑1 (MQ=255)                                                                             ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACa                                                                         <  1:1324891/151‑1 (MQ=255)                                                                             ggCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACa                                                                         <  1:543735/151‑1 (MQ=255)                                                                                 ggAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTg                                                                     <  1:1352851/151‑1 (MQ=255)                                                                                   aagaagATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGGCGGGGGGGGGAATTATCAGC‑GATGGATTACTGCACAACAACCTGCGGGAATTCCTCCTGATCCAAAACTCACCCATTACGAACAGCTgg                                                                    >  2:1267453/1‑150 (MQ=255)                                                                                    agaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑ATGGGCGGGGAGGGGGATTATTAGC‑GATGGGTTTCTGCACAAAACCCTGCGGGAATTTCTCTGGATCAAAAAATCACCCATTACGAACAGCCgg                                                                    >  1:1351420/1‑149 (MQ=255)                                                                                     gaagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTTG‑TTTGTCCGGGGGGGGGATTATCAAC‑GATGGATTACTGCACCACAACCCGCGGGAATTTCTCTTGATCAAAAAATCACCCCTTACGaa                                                                           >  1:814252/1‑141 (MQ=255)                                                                                      aagaTGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGGTTACTGCACAACACCCTGCGGGAATTCCTCTGGATCAAAAAATCACCCATTACGAACAGCTgg                                                                    >  2:139387/1‑147 (MQ=255)                                                                                           gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGGTTTGGTGGGGGGGGGGAA‑TATTAGC‑GGTGGATTACTTCACCACACCCTGCGGGAATTTCTTTGGGTCAAAAAATCACCCATTACGAACAGCTGGCaa                                                                 >  2:964313/1‑145 (MQ=255)                                                                                           gtgCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGAt                                                            <  2:1244384/150‑1 (MQ=255)                                                                                            tgctaaaagCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGGATTATCAGC‑GATGGGTTACTGCCCAACACCCTGCGGGAATTTCTCCGGGTCAAAAAATCACCCCTTACGAACAGCCGGCaa                                                                 >  2:539592/9‑144 (MQ=255)                                                                                              ctctACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATTAGC‑GATGGATTACTGCACAAAACCCCGCGGAAATTCCTCTGGATCAAAAAATCACCCATTACGAAAAGCTGGCAAAGGATGCAc                                                        >  2:697338/1‑151 (MQ=255)                                                                                                   cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGAATTATTAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTTCTCTGGGTCAAAAAATCACCCATTACGAACAGCTGGCAAGTGATGCACTgg                                                     >  2:585176/1‑149 (MQ=255)                                                                                                   cGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCa                                                   <  1:671766/151‑1 (MQ=255)                                                                                                    gCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAg                                                  <  2:178078/151‑1 (MQ=255)                                                                                                         gACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGgag                                               <  2:929901/149‑1 (MQ=255)                                                                                                         gACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGaa                                             <  1:628429/151‑1 (MQ=255)                                                                                                            cGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTGGGGGGGGGGAATTATCAGC‑GGTGGATTACTGCACAAAACCCCGCGGGAATTTCTCTGGATCAAAAACTCACCCATTACCAACAGCCGGCAACAGATGCACAGGCCGGAGTAAtt                                          >  1:1168272/1‑151 (MQ=255)                                                                                                                tttAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATtgtg                                       <  2:1342199/150‑1 (MQ=255)                                                                                                                      ggTGAATATACTGGCGGTTGG‑TTTGTCGGGGGGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGGAATTCCTCTTGATCAAAAACTCACCCATTACGAACAGCTGGCCAGTGGTGCACTGGCAGGAGAAAATGTggggtttt                                >  1:1140105/1‑146 (MQ=255)                                                                                                                         gAATATACTGGCGGTTGG‑TTGGTCGGGGGGGGGAATTTTCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAAAAGCTGGCAAGTGATGCACTGGCAGGGGGAATTGTGGGttctt                                >  1:818753/1‑148 (MQ=255)                                                                                                                         gAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTtcttct                              <  2:486719/150‑1 (MQ=255)                                                                                                                          aaTATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgc                            <  2:1063045/151‑1 (MQ=255)                                                                                                                           atatACTGGCGGTTGG‑TTTGTCGGGGGGGGGGATTATCAGC‑GATGGATTACTGGACAACACCCCGCGGGAATTTCTCTGGATCAAAAACTCACCCCTTACAAACAACTGGCAAGTGATGCACTGGCCGGAGGAAATGTggggtttt                                >  1:229181/1‑141 (MQ=255)                                                                                                                                cTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggcgg                      >  1:881948/1‑151 (MQ=255)                                                                                                                                 tGGCGGTTGG‑TTGGACGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTgcgggcgg                      >  2:240272/1‑150 (MQ=255)                                                                                                                                  ggCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg                    <  2:1126139/151‑1 (MQ=255)                                                                                                                                  ggCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg                    <  2:304387/151‑1 (MQ=255)                                                                                                                                    cGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGgcggcg                    >  2:581932/1‑149 (MQ=255)                                                                                                                                     ggttgg‑ttGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGagc                 >  2:113541/1‑151 (MQ=255)                                                                                                                                         gg‑ttggtCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGcc             >  2:487312/1‑151 (MQ=255)                                                                                                                                              ggtCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAGTTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCaaa          >  2:666764/1‑150 (MQ=255)                                                                                                                                               gtCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAt        <  2:49259/151‑1 (MQ=255)                                                                                                                                                 cGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt      >  2:1038999/1‑151 (MQ=255)                                                                                                                                                 cGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt      <  2:1170323/151‑1 (MQ=255)                                                                                                                                                  ggggaggggaATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAAttt      >  2:359938/1‑150 (MQ=255)                                                                                                                                                   gggaggggaATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTAc    <  2:1127357/151‑1 (MQ=255)                                                                                                                                                   gggaggggaATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTAc    <  1:430928/151‑1 (MQ=255)                                                                                                                                                      aggggaATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACct  >  2:559003/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                       CCGCAGGGCTTTTTAGGGCGCGCGTGGGGCAGGAAGTTAGCCGCGCAACTGAATCTGACTGATGATATACGTCTCTGGCAGGAAGAAGATGTGCTCTACGCCCTGACCGTATTTAACGGTGAATATACTGGCGGTTGG‑TTGGTCGGGGAGGGGAATTATCAGC‑GATGGATTACTGCACAACACCCTGCGGAAATTCCTCTGGATCAAAAACTCACCCATTACGAACAGCTGGCAAGTGATGCACTGGCAGGAGAAATTGTGGGTTCTTCTGCGGGCGGCGAGCAGCCAAAATTTACCT  >  NC_000913/4622102‑4622399

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