Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,305,427	A→G	Y391H (TAC→CAC)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,305,427	0	A	G	96.9%	213.5 / NA	65	Y391H (TAC→CAC)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base G (23/39);  minor base C (2/0);  total (26/39)
Fisher's exact test for biased strand distribution p-value = 1.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGT  >  NC_000913/4305292‑4305561                                                                                                                                        |                                                                                                                                       cgtcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGccc                                                                                                                                    >  1:671638/1‑140 (MQ=255)  gtcgtAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCt                                                                                                                                   >  2:154578/1‑140 (MQ=255)     gtAATAGCCAAGATAGATGTTATACACCGCGTGGGGGTTCTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCATAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGt                                                                                                                                <  2:277840/140‑1 (MQ=255)         tAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGc                                                                                                                            <  2:462658/140‑1 (MQ=255)         tAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGc                                                                                                                            <  1:852770/140‑1 (MQ=255)           gCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCa                                                                                                                          <  1:154578/140‑1 (MQ=255)               aGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTg                                                                                                                      <  1:656682/140‑1 (MQ=255)                   aGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAgg                                                                                                                  >  2:444986/1‑140 (MQ=255)                          ataCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGt                                                                                                           >  2:8862/1‑140 (MQ=255)                          ataCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGACGATTTCATTCAGGGTGTGGCCCTGGTTCGACAGGTGCAGGGTCTGGt                                                                                                           >  2:281077/1‑140 (MQ=255)                           taCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCTGATTTCATGCATAGTGCGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTc                                                                                                          >  2:321046/1‑140 (MQ=255)                           taCACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTc                                                                                                          <  1:67572/140‑1 (MQ=255)                               ccGCGCGGGCGTTGTGGGTGACCGAACCGTACTAGACGCCGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGa                                                                                                      <  2:312374/140‑1 (MQ=255)                                cgcgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGaa                                                                                                     >  2:648070/1‑140 (MQ=255)                                 gcgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAAt                                                                                                    >  1:579553/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGGTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATg                                                                                                   >  2:673154/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGTAGGGTCTGGTCGTGAATg                                                                                                   >  2:125701/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATg                                                                                                   >  1:812415/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGGTCGACACGTGCAGCGTCTGGTCGTGAATg                                                                                                   >  1:321103/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTCGCCCTGGTTCGACAGGTGCAGGGTCTGGTCGTGAATg                                                                                                   >  1:346411/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTCGCCCTGGATCGGCAAGTGCAGCGTCTGGTCGTGAATg                                                                                                   >  1:312374/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGATCGTCATGTGCAGGGTCTGGTCGTGAATg                                                                                                   >  1:303464/1‑140 (MQ=255)                                  cgcgGGCGTTGTGGCGGACAGAACCCTAAGAGCCGCGTCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATg                                                                                                   <  1:374737/140‑1 (MQ=255)                                    cgGGCGTTGTGGCTGACTGAACCGTAATAGCTGCGCCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTa                                                                                                 <  2:289335/140‑1 (MQ=255)                                         gTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGa                                                                                            <  1:648070/140‑1 (MQ=255)                                         gTTGTGGCTGACAGAACCGTAAAAGCCGCGGATGGTCCATTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGa                                                                                            <  2:346411/140‑1 (MQ=255)                                         gTTGTGGCTGACAGAACCGTAAAAGCCGCAGCTGGCCTAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGa                                                                                            <  1:274765/140‑1 (MQ=255)                                          ttGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCGGGATCTGGTAGTGAATGTACTTGAt                                                                                           >  1:234359/1‑140 (MQ=255)                                                 tGACAGAACCGTAATAGCCGCGACTGGCCCAGTTAATGGAAAGTGCAGCCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcg                                                                                    <  2:312919/140‑1 (MQ=255)                                                   aCAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcgcg                                                                                  <  2:636200/140‑1 (MQ=255)                                                   aCAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATcgcg                                                                                  <  1:805026/140‑1 (MQ=255)                                                      gAACCGTTGTGGCCGGAGCTGGGCCAGTTATGGTGCAGTGCAGTCGGCGGCTTAATCATGTCGCCGATTTCATTCATGGTGTTGCCCTGGTTCGCCAGGTGGAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTa                                                                               <  2:389738/140‑1 (MQ=255)                                                        aCCGTAATAGCCGCGGCTGGCCCAGCTATAGGCAAGTCCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTAtt                                                                             <  1:379835/140‑1 (MQ=255)                                                          cGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTa                                                                           <  2:73385/140‑1 (MQ=255)                                                           gTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTAc                                                                          <  1:66789/140‑1 (MQ=255)                                                             aaTAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCTGGGGGTGGACGTGACTGTCCTAGATGGTTTCGCGGTATTTACCa                                                                        >  1:419581/1‑140 (MQ=255)                                                                aGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAata                                                                     <  1:483261/140‑1 (MQ=255)                                                                    gcgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAAt                                                                 >  2:526330/1‑140 (MQ=255)                                                                    gcgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAAt                                                                 <  1:424391/140‑1 (MQ=255)                                                                    gcgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAAt                                                                 >  2:424391/1‑140 (MQ=255)                                                                    gcgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAAt                                                                 <  2:812415/140‑1 (MQ=255)                                                                     cgGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTTATGGTAGCGCGGTATTTACCAATATAATc                                                                >  1:299478/1‑140 (MQ=255)                                                                      ggCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCa                                                               <  2:671638/140‑1 (MQ=255)                                                                       gCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                    <  1:790711/134‑1 (MQ=255)                                                                       gCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAatat                                                                    >  2:790711/1‑134 (MQ=255)                                                                              cAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATg                                                       <  2:579553/140‑1 (MQ=255)                                                                                  tATTGGCAGGTGCAGGGTGCAGCTTTATCACGTCGACGAGTTCATTCAGGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                   <  2:419581/140‑1 (MQ=255)                                                                                  tATTGGCAAGTGCAGACGGCAGCTTAATCTTGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTa                                                   <  1:344134/140‑1 (MQ=255)                                                                                   aTTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTAt                                                  <  1:552434/140‑1 (MQ=255)                                                                                   aTTGGCAAGTGCAGGCGACAGCGTAATCATGGCGCCGATCTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTAt                                                  <  2:299478/140‑1 (MQ=255)                                                                                           gTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAg                                          >  2:132830/1‑140 (MQ=255)                                                                                           gTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAg                                          >  2:691756/1‑140 (MQ=255)                                                                                            tGCAGGCGGCAGCTTTATCATGTCGTCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                         <  1:349668/140‑1 (MQ=255)                                                                                            tGCAGGCGGCAGCTTTATCATGTCGTCGACTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                         <  1:321046/140‑1 (MQ=255)                                                                                            tGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGTAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                         <  1:125701/140‑1 (MQ=255)                                                                                            tGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                         <  1:673154/140‑1 (MQ=255)                                                                                            tGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGa                                         <  1:526330/140‑1 (MQ=255)                                                                                              cAGGCGGCAGCTTAATCATGTCGCCAATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGAcc                                       >  1:838002/1‑140 (MQ=255)                                                                                                ggcggcAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCgg                                     <  1:444986/140‑1 (MQ=255)                                                                                                 gcggcAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGc                                    >  1:129726/1‑140 (MQ=255)                                                                                                  cggcAGCTTAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGcc                                   <  1:8862/140‑1 (MQ=255)                                                                                                          tAATCATGTCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGc                           <  2:129726/140‑1 (MQ=255)                                                                                                                  tCGCCGATTTCATTCATGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATaaa                   >  1:213503/1‑140 (MQ=255)                                                                                                                           tcattcatGGTGTGGCCCTGGTACGCCAGGTCCAGTGTCCGGTCGTGACTGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTc          <  2:234359/140‑1 (MQ=255)                                                                                                                            cattcatGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTcc         <  1:691756/140‑1 (MQ=255)                                                                                                                                catGGTGTGGCCCTGGTACGCCGGGTGCAGGATCTAGTCCTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgt     <  2:303464/140‑1 (MQ=255)                                                                                                                                  tGGTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGcgtcg   <  1:132830/140‑1 (MQ=255)                                                                                                                                   ggTGTGGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCATGTGTGCGGGATAAGACGCACTTCCGcgtcgt  >  2:383608/1‑140 (MQ=255)                                                                                                                                        |                                                                                                                                       CGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCGCCGATTTCATTCATGGTGTAGCCCTGGTTCGCCAGGTGCAGGGTCTGGTCGTGAATGTACTTGATGGTATCGCGGTATTTACCAATATAATCATTGATATGCTTATTGCCCCAGACCGGCCAGGTGTGCGGCATAAACAGCACTTCCGCGTCGT  >  NC_000913/4305292‑4305561

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