Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,068,440	G→A	G279G (GGC→GGT)	yihQ ←	sulfoquinovosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,068,440	0	G	A	98.5%	209.3 / ‑4.2	68	G279G (GGC→GGT)	yihQ	sulfoquinovosidase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (37/30);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

ATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCCAGCAGGG  >  NC_000913/4068305‑4068576                                                                                                                                        |                                                                                                                                         aTTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcg                                                                                                                                      >  2:524184/1‑140 (MQ=255) aTTCCACTGC‑TAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcgc                                                                                                                                     <  1:461671/140‑1 (MQ=255)  ttCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcgc                                                                                                                                     <  1:624858/140‑1 (MQ=255)  ttCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcgc                                                                                                                                     <  2:407944/140‑1 (MQ=255)  ttCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTCTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcgc                                                                                                                                     <  1:216759/140‑1 (MQ=255)   tCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCcgcg                                                                                                                                    <  2:557726/140‑1 (MQ=255)      aCTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTa                                                                                                                                 >  2:696713/1‑140 (MQ=255)      aCTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTa                                                                                                                                 >  1:607429/1‑140 (MQ=255)      aCTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTa                                                                                                                                 >  1:598700/1‑140 (MQ=255)        tGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACTGGACCAGTGCTGGGCCTAGATGCCGATGACTTTCACACCCGCGTTACg                                                                                                                               >  2:284039/1‑140 (MQ=255)        tGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACg                                                                                                                               >  2:655476/1‑140 (MQ=255)          cTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTACAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCa                                                                                                                             <  2:615941/140‑1 (MQ=255)           ttAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCAt                                                                                                                            >  1:407458/1‑140 (MQ=255)             aaTGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACGAGTCCTGCGCCCAGATGCCGTTGACGTTCAAACCCGCGTTACGCATgg                                                                                                                          >  2:193225/1‑140 (MQ=255)                gCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGgtgt                                                                                                                       <  1:696713/140‑1 (MQ=255)                 cGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTc                                                                                                                      <  2:545518/140‑1 (MQ=255)                  gTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTcc                                                                                                                     >  1:55611/1‑140 (MQ=255)                    gAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACATTCGCACCCGCGTTACGCATGGTGTCCAg                                                                                                                   >  1:247788/1‑140 (MQ=255)                    gAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGAAGCCGTTGACCTTCACATCCGCTATACGCATGGTGTCCAg                                                                                                                   >  2:319816/1‑140 (MQ=255)                    gAATCCAGTTGCGGGGAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAg                                                                                                                   >  2:103196/1‑140 (MQ=255)                     aaTCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGt                                                                                                                  <  1:60011/140‑1 (MQ=255)                         cAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTAACCTTCACACCAGCATTACGCATGGTGTCCAGTttct                                                                                                              >  1:348119/1‑140 (MQ=255)                               cGGGTAGTTTTCGCTTTTCCACTTCCAGTTCCACATCTCGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGc                                                                                                        <  1:278472/140‑1 (MQ=255)                                gggTAGTTTTCGCTGTTCCACTTCCTGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGca                                                                                                       <  1:284039/140‑1 (MQ=255)                                      ttttCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTc                                                                                                 >  2:459925/1‑140 (MQ=255)                                          cGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTc                                                                                             <  1:655476/140‑1 (MQ=255)                                          cGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTc                                                                                             <  2:55611/140‑1 (MQ=255)                                          cGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTc                                                                                             <  1:123181/140‑1 (MQ=255)                                                tCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGccc                                                                                       >  1:605988/1‑140 (MQ=255)                                                        cAGTTCCACATCACGTGTTTGCAAAAAGAGGTCATACCAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGcc                                                                               <  1:258814/140‑1 (MQ=255)                                                           ttCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAg                                                                            <  2:446406/140‑1 (MQ=255)                                                               acaTCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGtt                                                                        >  2:64976/1‑140 (MQ=255)                                                               acaTCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGtt                                                                        >  2:440272/1‑140 (MQ=255)                                                                 aTCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTAc                                                                      >  2:88958/1‑140 (MQ=255)                                                                    aCGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTcc                                                                   <  2:550485/140‑1 (MQ=255)                                                                        gTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCa                                                               >  2:596794/1‑140 (MQ=255)                                                                           tGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGTCCTGAAGGCCGAGGGTTACTCCGTCATaa                                                            >  1:226319/1‑140 (MQ=255)                                                                           tGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTAAAAGCCGAGGGTTACTCCGTCATaa                                                            >  1:314782/1‑140 (MQ=255)                                                                                      gTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCggg                                                 <  1:440272/140‑1 (MQ=255)                                                                                      gTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCAGAATGCCAAGCTTTGCTCCGTAATAAATCCAGTCggg                                                 >  1:353683/1‑140 (MQ=255)                                                                                               aTACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTctgg                                        >  1:549229/1‑140 (MQ=255)                                                                                                tACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAg                                                      >  1:10164/1‑125 (MQ=255)                                                                                                tACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAg                                                      <  2:10164/125‑1 (MQ=255)                                                                                                 aCCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTctggct                                      <  2:607429/140‑1 (MQ=255)                                                                                                 aCCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTctggct                                      >  1:638807/1‑140 (MQ=255)                                                                                                    ggACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCtggctggc                                   >  1:6332/1‑140 (MQ=255)                                                                                                     gACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCggg                                                 >  1:619230/1‑125 (MQ=255)                                                                                                     gACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCggg                                                 <  2:619230/125‑1 (MQ=255)                                                                                                        cAGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTcc                               >  1:599559/1‑140 (MQ=255)                                                                                                         aGTCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAGTGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTccc                              >  1:560554/1‑140 (MQ=255)                                                                                                           tCCTGCGCCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCcag                            >  1:94897/1‑140 (MQ=255)                                                                                                                cgcCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTGCTCCGTCACAAAACCAAGCGGGCAGATCTGACTGGCGTCTCAGCAggg                       >  2:358906/1‑140 (MQ=255)                                                                                                                cgcCCAGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAggg                       >  2:705010/1‑140 (MQ=255)                                                                                                                  cccAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCg                     >  1:478603/1‑140 (MQ=255)                                                                                                                     aGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGTCAGTTCTGGCTGGCGTCCCAGCAGGGCGGtt                  >  2:279656/1‑140 (MQ=255)                                                                                                                     aGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGtt                  >  1:604876/1‑140 (MQ=255)                                                                                                                     aGATGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGtt                  >  1:424646/1‑140 (MQ=255)                                                                                                                       aTGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTaa                >  1:723367/1‑140 (MQ=255)                                                                                                                       aTGCCGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGCTCTAGCAGGCGTCTCAGCAGGGCGGTTaa                >  1:343765/1‑140 (MQ=255)                                                                                                                           cGTTGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAAtttt            <  1:88958/140‑1 (MQ=255)                                                                                                                             ttGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTc          <  2:407458/140‑1 (MQ=255)                                                                                                                             ttGACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTc          <  1:459925/140‑1 (MQ=255)                                                                                                                               gACCTTCACACCCGCTTTACGCATGTTGTACAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCca        <  1:358906/140‑1 (MQ=255)                                                                                                                               gACCTTCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCca        <  2:604876/140‑1 (MQ=255)                                                                                                                                 ccTTCACATCCGCGTTGCGCATTGTGTCCATTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCcagc      <  2:343765/140‑1 (MQ=255)                                                                                                                                  cTTCACACCCGCGTTACACATGGTGTCCAGTGTGTCCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCcagca     <  2:348119/140‑1 (MQ=255)                                                                                                                                   ttCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCcagcag    <  2:599559/140‑1 (MQ=255)                                                                                                                                   ttCACACCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCcagcag    <  2:605988/140‑1 (MQ=255)                                                                                                                                     cacaCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCCAGCAggg  <  1:524184/140‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         ATTCCACTGCTTAATGCGTGAATCCAGTTGCGGGTAGTTTTCGCTGTTCCACTTCCAGTTCCACATCACGCGTTTGCCAAAAGAGGTCATACGAATACCGGACCAGTCCTGCGCCCAGATGCCGTTGACCTTCACGCCCGCGTTACGCATGGTGTCCAGTTTCTTCTGGCACACTTCCGTCCCGCCCTGAATGCCGAGCGTTACTCCGTCATAAATCCAGTCGGGCAGTTCTGGCTGGCGTCCCAGCAGGGCGGTTAATTTTTCCAGCAGGG  >  NC_000913/4068305‑4068576

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