Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,507,501	A→C	5.9%	I241L (ATT→CTT)	ydcP →	putative peptidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,507,501	0	A	C	5.9%	281.1 / 6.9	82	I241L (ATT→CTT)	ydcP	putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (39/37);  new base C (2/3);  total (42/40)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAA  >  NC_000913/1507353‑1507649                                                                                                                                                     |                                                                                                                                                      ttGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGa                                                                                                                                                       <  2:698534/149‑1 (MQ=255)      cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                                  <  1:235408/149‑1 (MQ=255)      cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                                  <  1:812987/149‑1 (MQ=255)      cGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGaa                                                                                                                                                  <  3:547312/149‑1 (MQ=255)       gCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAg                                                                                                                                                 <  1:363594/149‑1 (MQ=255)           gCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCg                                                                                                                                             <  2:473392/149‑1 (MQ=255)            cGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGt                                                                                                                                            <  1:897257/149‑1 (MQ=255)              tGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTa                                                                                                                                          >  2:688383/1‑149 (MQ=255)                  gTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaaa                                                                                                                                      <  4:696206/149‑1 (MQ=255)                  gTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaaa                                                                                                                                      <  3:300579/149‑1 (MQ=255)                    ttGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGa                                                                                                                                    <  2:24150/149‑1 (MQ=255)                     tGCCATACCCAGTGAAAGACGGTCAGGGGCGGGTGGTTTCATATGACAAACATCTGCTGTCGGTGAAACATAACGATAAGACTACCAACCTCGGCGCGCTATATGATGCTGGTGTACGCTACTTTAAGATTGAAGGGCGTTACAAAGat                                                                                                                                   <  4:312756/149‑1 (MQ=255)                       ccATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGatat                                                                                                                                 <  4:376384/149‑1 (MQ=255)                         aTACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                               >  2:878877/1‑149 (MQ=255)                          tACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAg                                                                                                                              >  3:160773/1‑149 (MQ=255)                             acaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaaa                                                                                                                                      <  3:153485/138‑1 (MQ=255)                             acaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACaaa                                                                                                                                      >  4:153485/1‑138 (MQ=255)                             acaTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCAAACCTCGGCGAGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTa                                                                                                                           >  4:240397/1‑149 (MQ=255)                              caTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTAc                                                                                                                          <  3:484230/149‑1 (MQ=255)                              caTTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTAc                                                                                                                          <  3:779983/149‑1 (MQ=255)                                ttGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGt                                                                                                                        <  4:862540/149‑1 (MQ=255)                                  gAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTATCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGATACAAAGATATGAGCTACGTg‑a                                                                                                                     >  1:283762/1‑149 (MQ=255)                                   aaaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTg‑aa                                                                                                                    <  2:851537/149‑1 (MQ=255)                                    aaGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTg‑aag                                                                                                                   >  3:222165/1‑149 (MQ=255)                                     aGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTg‑aaga                                                                                                                  >  3:327555/1‑149 (MQ=255)                                         gATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAatat                                                                                                              >  3:241378/1‑149 (MQ=255)                                         gATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAatat                                                                                                              <  4:243269/149‑1 (MQ=255)                                              ggggCGGGGGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCg                                                                                                         >  4:76891/1‑149 (MQ=255)                                               gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGc                                                                                                        >  4:105685/1‑149 (MQ=255)                                               gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGc                                                                                                        >  1:884222/1‑149 (MQ=255)                                               gggCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGc                                                                                                        >  1:272612/1‑149 (MQ=255)                                                    ggtggtTTCTTATGAGAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGCTTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCAtt                                                                                                   <  2:157721/149‑1 (MQ=255)                                                     gtggtTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGCCTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTa                                                                                                  <  2:388815/149‑1 (MQ=255)                                                          ttCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGcc                                                                                             <  2:272612/149‑1 (MQ=255)                                                              tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                               >  4:654262/1‑112 (MQ=255)                                                              tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGa                                                                                                                               <  3:654262/112‑1 (MQ=255)                                                              tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGAt                                                                                         >  2:117366/1‑149 (MQ=255)                                                              tATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGAt                                                                                         >  1:35145/1‑149 (MQ=255)                                                                  aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCtt                                                                                     <  4:25987/149‑1 (MQ=255)                                                                  aaaaaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCtt                                                                                     >  1:635541/1‑149 (MQ=255)                                                                     aaCATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGAt                                                                                  <  1:265010/149‑1 (MQ=255)                                                                        aTCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGcc                                                                               <  4:129426/149‑1 (MQ=255)                                                                              tgTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCattatt                                                                         <  3:277575/149‑1 (MQ=255)                                                                               gTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTg                                                                        >  2:643280/1‑149 (MQ=255)                                                                                  gATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaag                                                                     >  2:421307/1‑149 (MQ=255)                                                                                   aTGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTgaaga                                                                    >  2:54996/1‑149 (MQ=255)                                                                                    tGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGTTGGTTCCGGGTTTCCGCTCTTTAAAGTTTGAAGGCGGTTCCAAAG‑TATGAACTTAGTGAAAGAATATAACCGCCCTTTTTCCCCAATTCCTTGTTGCCATTTTTgaaaaa                                                                   >  3:538498/1‑149 (MQ=255)                                                                                    tGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAg                                                                                                                                                 <  1:705991/72‑1 (MQ=255)                                                                                    tGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAATCTTGAAg                                                                                                                                                 >  2:705991/1‑72 (MQ=255)                                                                                     gAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAAc                                                                  >  4:260569/1‑149 (MQ=255)                                                                                         gATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                              >  2:762870/1‑149 (MQ=255)                                                                                         gATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTgg                                                              >  1:83367/1‑149 (MQ=255)                                                                                                aTCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTg                                                       <  1:117366/149‑1 (MQ=255)                                                                                                    gACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgc                                                   <  2:79863/149‑1 (MQ=255)                                                                                                          cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCat                                                                             >  1:688769/1‑117 (MQ=255)                                                                                                          cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCat                                                                             <  2:688769/117‑1 (MQ=255)                                                                                                          cAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCtt                                             >  2:37638/1‑149 (MQ=255)                                                                                                              cTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatc                                         >  4:599325/1‑149 (MQ=255)                                                                                                               tCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatca                                        <  3:691839/149‑1 (MQ=255)                                                                                                                cGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAg                                       <  1:643280/149‑1 (MQ=255)                                                                                                                   cgcgcTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTc                                    <  4:272219/149‑1 (MQ=255)                                                                                                                      gcTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTa                                 <  3:535896/149‑1 (MQ=255)                                                                                                                      gcTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTa                                 >  2:587961/1‑149 (MQ=255)                                                                                                                        tgattgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACt                               <  4:222165/149‑1 (MQ=255)                                                                                                                          attgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatca                                        <  2:7563/138‑1 (MQ=255)                                                                                                                          attgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTtcatca                                        >  1:7563/1‑138 (MQ=255)                                                                                                                          attgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGa                             >  1:50782/1‑149 (MQ=255)                                                                                                                            tgatGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAAc                           >  3:844620/1‑149 (MQ=255)                                                                                                                               tGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACAtt                        >  4:686153/1‑149 (MQ=255)                                                                                                                                 cTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTc                      >  3:210151/1‑149 (MQ=255)                                                                                                                                  tGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                     <  4:210151/149‑1 (MQ=255)                                                                                                                                  tGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCt                     <  2:635541/149‑1 (MQ=255)                                                                                                                                   ggTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCtt                    >  2:224740/1‑149 (MQ=255)                                                                                                                                         acgcccTTCAATCTTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgcgcg                                                <  2:54956/111‑1 (MQ=255)                                                                                                                                         acgcccTTCAATCTTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGgcgcgcg                                                >  1:54956/5‑115 (MQ=255)                                                                                                                                          gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCa             >  4:130156/1‑149 (MQ=255)                                                                                                                                          gCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCa             >  3:107034/1‑149 (MQ=255)                                                                                                                                           cTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCAt            <  2:35145/149‑1 (MQ=255)                                                                                                                                               ttCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGAc        >  1:339913/1‑149 (MQ=255)                                                                                                                                                  aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGaa                            >  4:11279/1‑126 (MQ=255)                                                                                                                                                  aaGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGaa                            <  3:11279/126‑1 (MQ=255)                                                                                                                                                   aGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGaa    >  4:522818/1‑149 (MQ=255)                                                                                                                                                    gATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGaaa   >  4:836032/1‑149 (MQ=255)                                                                                                                                                     aTTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTATTGAACATTTCTTTGTTCCATCTACGGaaaa  >  2:730485/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TTGCTCGCAGGCGTGCCGTTTGCCATACACATTGAAAGACGATCAGGGGCGGGTGGTTTCCTATGAAAAACATCTGCTGTCGATGAAAGATAACGATCAGACTGCCAACCTCGGCGCGCTGATTGATGCTGGTGTACGCTCCTTCAAGATTGAAGGGCGTTACAAAGATATGAGCTACGTG‑AAGAATATCACCGCCCATTATCGCCAGATGCTTGATGCCATTATTGAAGAACGTGGCGATCTGGCGCGCGCTTCATCAGGTCGTACTGAACATTTCTTTGTTCCATCGACGGAAAA  >  NC_000913/1507353‑1507649

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑