Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,500	0	A	G	40.3%	102.3 / 23.0	68	pseudogene (1087/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (6/34);  new base G (27/0);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTC  >  NC_000913/4573354‑4573647                                                                                                                                                     |                                                                                                                                                         aaGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGgatgat                                                                                                                                                      <  2:776825/151‑1 (MQ=255)  aGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGAtt                                                                                                                                                     <  1:654623/151‑1 (MQ=255)    aTGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGTGAttg‑t                                                                                                                                                  >  2:1233336/1‑151 (MQ=255)        tGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTTGCGTTAGTGGCTCTGGTCT‑GGAGGGTGGGTTG‑TTgg                                                                                                                                               >  1:1243729/1‑151 (MQ=255)        tGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGGTTGTCGTTAGTGGCTC‑GGGGT‑GGGGGGGTGGTTG‑TTggg                                                                                                                                              >  1:1182285/1‑151 (MQ=255)          aGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGGGGGGTGGTTG‑TTggg                                                                                                                                              >  2:32710/1‑149 (MQ=255)          aGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAg                                                                                                                                            <  2:1295376/151‑1 (MQ=255)            gAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGTGT‑GGGGGGGGGGTTG‑TTggg                                                                                                                                              >  2:499663/1‑147 (MQ=255)             aaGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑At                                                                                                                                        <  2:83777/151‑1 (MQ=255)               gTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTa                                                                                                                                      <  1:172113/151‑1 (MQ=255)                   ctGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTaa                                                                                                                                     <  2:729847/148‑1 (MQ=255)                     ggCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGGT‑GGGGGGGTGGTTG‑TTTGGGGC‑ATTtaa                                                                                                                                    >  2:390435/1‑147 (MQ=255)                     ggCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GCGGT‑GGAGGGGTGGTTG‑GTGGGGGG‑ATTTAAAtt                                                                                                                                 >  2:1171722/1‑150 (MQ=255)                      gCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTtga                                                                                                                               <  2:388682/151‑1 (MQ=255)                       ccTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTtgat                                                                                                                              <  2:600388/151‑1 (MQ=255)                           gggCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACa                                                                                                                          <  2:1256778/151‑1 (MQ=255)                             gCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGAGGTTG‑TTGGGGGG‑ATTaaaaattgtgacag                                                                                                                         >  2:1206573/1‑139 (MQ=255)                                 gCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGaaa                                                                                                                   <  1:106335/151‑1 (MQ=255)                                     aTGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑g                                                                                                               <  2:899565/150‑1 (MQ=255)                                            tctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGGT‑GGGGGGGGGGATG‑TTGGGGGG‑CTTTAACTTGATGGCAGGTGAAA‑T‑GCCAtttg                                                                                                        >  1:274564/1‑148 (MQ=255)                                               ggAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTAAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑Ga                                                                                                    <  1:569522/150‑1 (MQ=255)                                                  aaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTg                                                                                                <  2:509632/151‑1 (MQ=255)                                                    acagacagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGcc                                                                                              <  1:1250549/151‑1 (MQ=255)                                                        acagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTa                                                                                          <  2:716652/151‑1 (MQ=255)                                                          agGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGGT‑GGGGGGGGGGTTT‑TTTGGGGG‑ATTTAAATTGATGAAAGGGGAAATT‑CCCTTTgggggggggc                                                                                                >  2:273965/1‑137 (MQ=255)                                                             cAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCC‑GGGGGGGGGGTTG‑TTGGGGGG‑GTTTAAATTTATGAAA‑GTGGAAAT‑GCCCTTGGG‑GAGGTGCCGTTAACCtt                                                                                     >  2:317150/1‑151 (MQ=255)                                                              aaCAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGCTAGTGGCCC‑GTGCG‑GGAGGGGGGGGTT‑CTGGGAGG‑ATTTCACTTGTTGACA‑GCGTAAAT‑GTCCTTTGT‑GGGGTGGCgcttaaac                                                                                      >  1:125225/1‑142 (MQ=255)                                                               aCAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGGGGATT‑TTGGGGGG‑CTTTAACTTTTTGACA‑GTGGAAAT‑GCCATTTGTGGGGGGGCCGTTAActt                                                                                      >  1:906957/1‑147 (MQ=255)                                                                         tACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGGGGTTG‑TTTGGGGC‑ATTAAAATTGATGAAA‑GGGAAAAT‑GCCATTTGT‑GAGGGGCCGGTAAACTTTATCGtt                                                                              >  2:987797/1‑146 (MQ=255)                                                                          aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGc                                                                         <  2:973907/150‑1 (MQ=255)                                                                             tAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GCGGT‑GGGGGGGGGGTTG‑TTGGGGGG‑CTTAAACTTGGTGACA‑GGGAAATTGGCCTTTGGG‑GGGGGCCCGTTAACCTCCAtg                                                                                 >  2:732691/1‑139 (MQ=255)                                                                                 aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                                                 <  1:824065/151‑1 (MQ=255)                                                                                 aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAAt                                                                 <  1:308347/151‑1 (MQ=255)                                                                                 aGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGGTGGGGGGGGGGTTTT‑TTGGGGGG‑ATTAAA‑TTTATTGCA‑GGGGAAAT‑TTCATTTGT‑GGGGTGCCGCTAACCCTTCTTGGTGATGCAGAACaa                                                                  >  1:769878/1‑150 (MQ=255)                                                                                  gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                                                <  1:1235602/151‑1 (MQ=255)                                                                                    gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATgg                                                               <  2:677175/150‑1 (MQ=255)                                                                                    gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGc                                                              <  2:182530/151‑1 (MQ=255)                                                                                    gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGGGGGTG‑TTGGGGGG‑ATTTAACTTGATGACA‑GGGGAATT‑TCCCTTTGT‑GAGGGGCCGTTAAACCTTATCGTTGATTCAGATCAATgg                                                               >  2:917316/1‑150 (MQ=255)                                                                                      aGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTg                                                            <  2:543440/151‑1 (MQ=255)                                                                                         gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                                          <  2:1111435/150‑1 (MQ=255)                                                                                               tGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGGGGGGTGGTTT‑TTGGGGGCAATTAAAC‑TGATGACAGGGGAAATT‑G‑CATTGGG‑GGGGGGGCGGTAAACTTCATTGTTGTTGCAGATCAATTGCTggttaaggg                                                    >  2:490316/1‑142 (MQ=255)                                                                                                ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGGT‑GGGGGGGTGGTTT‑TTGGGGGC‑GTTTAACTTGATGAGA‑GAGGAAAT‑GCCCCTTGG‑GAGGTGGCGGTAAACTTCCTCGTTGGTGCCAATCAATGGCTGttttagg                                                     >  2:380227/1‑143 (MQ=255)                                                                                                   aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGGT‑GGGGGGGGGATTG‑TTGGGGGG‑ATTTAAATTTGTGACA‑GGGGAAAT‑GTCATTGGT‑GGGGGGGCGTTTACCTTTATCGTTGGTGGAGATCAATGGGTggttcaag                                                     >  1:7144/1‑138 (MQ=255)                                                                                                   aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGTGGGTG‑TTTGGGGC‑CATAAACTTGATGGCA‑GTGGAAAT‑GGCAATGGG‑GGGGGGCCGTTAACCCTCATCGTTTGTGCAGATTAATGGGTGTTTCAGGGGCAgg                                               >  2:525826/1‑151 (MQ=255)                                                                                                       aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGCT‑GGGGGGGTGGTTG‑TTGGGGGG‑AATTAAATTGATGACA‑GTGGAAAT‑GCCCTTTGG‑GGGGGGGCGGTAAACTTCATCGTTGGTGCAGATCAAAGGGTGTTTCAGGGGCAGGctt                                            >  2:763625/1‑150 (MQ=255)                                                                                                            ccATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGc                                       <  1:756261/150‑1 (MQ=255)                                                                                                             cATCGGCAGGTTTGTCGTTAGTGGCTC‑GGGGT‑GGGGGGGTGATTG‑GTTGGGGG‑ATTAAAAATTATGACA‑GTGGAAAT‑GGCATTTGG‑GGGGGGGCGTTAAACTTTCTCGTTTATGCAGATCAAATGGTGTTTCAGGGAGAGGgttctt                                         >  1:1290092/1‑147 (MQ=255)                                                                                                                  gCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGGTGTTTGTTTGGGGGA‑ATTAAATTTG‑TGCCA‑GCGAAATT‑GCCATTTGG‑GGGGGGCCGTTAAACCTCATCTATGGTGCGGAACAATGGCTGTTGTAGGGACAAGCTTCTTGCGGagtgg                                >  1:490582/1‑150 (MQ=255)                                                                                                                  gCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGc                                <  1:683034/151‑1 (MQ=255)                                                                                                                   cAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCg                               <  1:1016990/151‑1 (MQ=255)                                                                                                                     ggTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCGTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa                             <  2:665930/151‑1 (MQ=255)                                                                                                                      gTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAag                            <  1:754840/151‑1 (MQ=255)                                                                                                                        ttGTCGTTAGTGGCTC‑GTGCT‑GGGGGGGTGGTTG‑TTGGGAGC‑ATTTAACTTGTTGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACACGCTTCTTGCCGggtggg                               >  1:924716/1‑144 (MQ=255)                                                                                                                        ttGTCGTTAGTGGCTC‑GTGCT‑GGGGGGATGATTG‑TTGGGGGC‑ATTTAACTTGATGACA‑GTGAAATT‑GCCATTGGG‑GAGGGGCCGTTAACCTTCATCGGTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCgga                             >  2:1118054/1‑146 (MQ=255)                                                                                                                         tGTCGTTAGTGGCTC‑GTGCT‑GGAGGGGTGATTG‑TTGGGGGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGGGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGaa                             >  2:904324/1‑147 (MQ=255)                                                                                                                          gTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagcagc                        <  2:1180087/151‑1 (MQ=255)                                                                                                                            cGTTAGTGGCTC‑GGGCT‑GGAGGGGTGTTTG‑TTGGGGGC‑ATTAAACTTGATGACA‑GTGGAATT‑GCCATTGGG‑GGGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGGaggcgcg                       >  1:13022/1‑145 (MQ=255)                                                                                                                                aGTGGCTC‑GTGCT‑GGGGGGGTGATTG‑TTGGGGGC‑CTTAAAATTGATGACA‑GTGAAATT‑GCCATTGGG‑GAGGTGGCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTTCAGGGACAGGCTTCTTGCAGGGTGGGGAGGAGCGATggt                  >  1:500743/1‑149 (MQ=255)                                                                                                                                aGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGtt                  <  2:956782/151‑1 (MQ=255)                                                                                                                                   ggCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTg               >  1:588601/1‑151 (MQ=255)                                                                                                                                   ggCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTg               >  2:722981/1‑151 (MQ=255)                                                                                                                                   ggCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTg               <  2:673324/151‑1 (MQ=255)                                                                                                                                          tGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgtg         >  2:165567/1‑151 (MQ=255)                                                                                                                                             cttgAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGCCA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGATAACCTTTATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCGTGCAGTGTGATAAGCGGTGATGTTATGATTGAGCt       >  1:661553/4‑151 (MQ=255)                                                                                                                                               ggAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTc      <  2:540563/150‑1 (MQ=255)                                                                                                                                                gAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTccc    <  2:1082066/151‑1 (MQ=255)                                                                                                                                                 aGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTccc    >  1:364935/1‑150 (MQ=255)                                                                                                                                                 aGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCt   <  1:757244/151‑1 (MQ=255)                                                                                                                                                  gggATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTc  <  1:474254/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                         AAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTC‑GTGCT‑GGAGGGATGATTG‑TTGGGAGC‑ATTAAACTTGATGACA‑GTGAAATT‑GCCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTC  >  NC_000913/4573354‑4573647

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