Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,643,960	T→C	Y408C (TAC→TGC)	prlC ←	oligopeptidase A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,643,960	0	T	C	98.6%	248.8 / ‑4.7	69	Y408C (TAC→TGC)	prlC	oligopeptidase A
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (32/36);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 4.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTAA  >  NC_000913/3643824‑3644095                                                                                                                                         |                                                                                                                                        aTTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGa                                                                                                                                      <  1:1143848/140‑1 (MQ=255)           gggTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACcgcg                                                                                                                            <  1:242925/139‑1 (MQ=255)           cGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACcgcg                                                                                                                            <  1:1128691/140‑1 (MQ=255)            ggTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCgcgc                                                                                                                           <  2:1028847/140‑1 (MQ=255)                  aGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCg                                                                                                                     <  2:83772/140‑1 (MQ=255)                  aGTTAAAAGTCAAATGCGCGACCGGTTTTTGCAAAGACCCATCAGCTTTACTCGTCTGGCCTACGCATTCATCCATCCCCGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCg                                                                                                                     <  2:1004174/140‑1 (MQ=255)                      aCAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTAt                                                                                                                 <  1:476053/140‑1 (MQ=255)                        aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTAttt                                                                                                               >  1:142724/1‑140 (MQ=255)                        aaGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTAttt                                                                                                               <  2:286303/140‑1 (MQ=255)                              aaTACGCGACCGGTTTTTGTAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCa                                                                                                         >  2:643476/1‑140 (MQ=255)                              aaTACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCa                                                                                                         >  1:964960/1‑140 (MQ=255)                                tACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATa                                                                                                       <  2:1071823/140‑1 (MQ=255)                                  cgcgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACa                                                                                                     >  2:167253/1‑140 (MQ=255)                                    cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGt                                                                                                   >  1:75957/1‑140 (MQ=255)                                    cgACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGt                                                                                                   >  2:688887/1‑140 (MQ=255)                                     gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCTTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGtt                                                                                                  <  2:765664/140‑1 (MQ=255)                                     gACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCTAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGtt                                                                                                  >  1:321193/1‑140 (MQ=255)                                                 cAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTa                                                                                      <  1:98163/140‑1 (MQ=255)                                                     gAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATc                                                                                  <  2:428376/140‑1 (MQ=255)                                                       aCCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCg                                                                                <  1:688887/140‑1 (MQ=255)                                                            cAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGc                                                                           <  2:783268/140‑1 (MQ=255)                                                            cAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGc                                                                           <  2:142724/140‑1 (MQ=255)                                                             aGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGcc                                                                          <  2:428583/140‑1 (MQ=255)                                                                tttACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATATAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGa                                                                       <  2:54513/140‑1 (MQ=255)                                                                tttACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGa                                                                       >  2:746337/1‑140 (MQ=255)                                                                    cGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATc                                                                   >  2:722608/1‑140 (MQ=255)                                                                         cTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACAt                                                              >  1:78800/1‑140 (MQ=255)                                                                         cTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACAt                                                              >  1:549859/1‑140 (MQ=255)                                                                           ggCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCt                                                            <  2:444043/140‑1 (MQ=255)                                                                                aCGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACg                                                       <  2:772673/140‑1 (MQ=255)                                                                                aCGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACg                                                       <  1:582686/140‑1 (MQ=255)                                                                                 cGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGc                                                      >  1:456057/1‑140 (MQ=255)                                                                                 cGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGc                                                      >  2:564050/1‑140 (MQ=255)                                                                                 cGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGc                                                      >  1:648678/1‑140 (MQ=255)                                                                                 cGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGc                                                      >  1:53406/1‑140 (MQ=255)                                                                                  gCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGct                                                     <  2:964960/140‑1 (MQ=255)                                                                                     gTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCtt                                                  <  1:167253/140‑1 (MQ=255)                                                                                      tCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACg                                                       <  2:862284/134‑1 (MQ=255)                                                                                      tCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACg                                                       >  1:862284/1‑134 (MQ=255)                                                                                      tCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCttt                                                 >  1:804710/1‑140 (MQ=255)                                                                                       catccatcCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTa                                                >  2:787697/1‑140 (MQ=255)                                                                                          ccatccaCGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCg                                             <  2:804710/140‑1 (MQ=255)                                                                                              ccaCGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGa                                         <  1:499882/140‑1 (MQ=255)                                                                                                aCGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATg                                       >  2:155234/1‑140 (MQ=255)                                                                                                    cccGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGt                                   >  1:463289/1‑140 (MQ=255)                                                                                                      cgccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaa                                 <  2:549859/140‑1 (MQ=255)                                                                                                      cgccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaa                                 <  2:648678/140‑1 (MQ=255)                                                                                                       gccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaaa                                >  1:1038747/1‑140 (MQ=255)                                                                                                       gccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaaa                                >  2:1070788/1‑140 (MQ=255)                                                                                                       gccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaaa                                >  2:82396/1‑140 (MQ=255)                                                                                                       gccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaaa                                <  1:418473/140‑1 (MQ=255)                                                                                                       gccgcGCTTGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTaaa                                <  2:1147795/140‑1 (MQ=255)                                                                                                         cgcgCTTGTTTTCACGGACATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATa                              <  2:1026279/140‑1 (MQ=255)                                                                                                              ttGTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCtt                         <  1:186530/140‑1 (MQ=255)                                                                                                                gTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTaa                       <  1:787697/140‑1 (MQ=255)                                                                                                                gTTTTCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTaa                       <  2:75957/140‑1 (MQ=255)                                                                                                                 ttttCACGGGCATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAAc                      <  1:643476/140‑1 (MQ=255)                                                                                                                           cATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            >  2:771447/1‑140 (MQ=255)                                                                                                                           cATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            >  1:153408/1‑140 (MQ=255)                                                                                                                           cATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            >  1:870354/1‑140 (MQ=255)                                                                                                                           cATACAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            >  2:66235/1‑140 (MQ=255)                                                                                                                              aCAGATCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAgg         <  2:78800/140‑1 (MQ=255)                                                                                                                                  aTCGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGt     >  1:885692/1‑140 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            <  2:640366/131‑1 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAAc            >  1:640366/1‑131 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTa   >  2:368152/1‑140 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTa   >  1:65265/1‑140 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTa   >  2:1127409/1‑140 (MQ=255)                                                                                                                                    cGAGGCAGAAGCTACCGCGCAGTTAGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTa   >  2:360018/1‑140 (MQ=255)                                                                                                                                     gAGGCAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTaa  <  1:746337/140‑1 (MQ=255)                                                                                                                                         |                                                                                                                                        ATTTACCGGGCGGTTGAAGTTACAAGTCAAATACGCGACCGGTTTTTGCAAAGAACCATCAGCTTTACGCATCTGGCCTACGCAGTCATCCATCCACGCCCCGCCGCGCTTGTTTTCACGGGCATACAGATCGAGGTAGAAGCTACCGCGCAGTTCGTTATTTTCGTCATACAGTTCGAAGAAACGTACATCCGGATGCCAGACATCAACATCTTTACGCTCTTTAGCGGTGATGCCGTAAATACGCTTAACCACTTCAAACAGGCCGTTAA  >  NC_000913/3643824‑3644095

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