Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,961,726	C→T	R350R (CGC→CGT)	rep →	DNA helicase and single‑stranded DNA‑dependent ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,961,726	0	C	T	100.0%	251.0 / NA	67	R350R (CGC→CGT)	rep	DNA helicase and single‑stranded DNA‑dependent ATPase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (32/35);  total (32/35)

GTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGA  >  NC_000913/3961582‑3961872                                                                                                                                                 |                                                                                                                                                   gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGt                                                                                                                                                 <  2:603533/148‑1 (MQ=255) gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTa                                                                                                                                                <  2:487711/149‑1 (MQ=255) gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAAc                                                                                                                                              <  2:664402/151‑1 (MQ=255) gTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAAc                                                                                                                                              >  1:372189/1‑151 (MQ=255)      ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCa                                                                                                                                            <  2:54152/148‑1 (MQ=255)      ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATc                                                                                                                                          <  2:445523/150‑1 (MQ=255)      ccGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCa                                                                                                                                         <  2:603672/151‑1 (MQ=255)        gAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGt                                                                                                                                       <  1:448514/151‑1 (MQ=255)         aaCTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTc                                                                                                                                      >  1:640385/1‑151 (MQ=255)           cTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTcgc                                                                                                                                    >  1:422957/1‑151 (MQ=255)                ttATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGgtg                                                                                                                               >  2:349872/1‑151 (MQ=255)                 tATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGgtgt                                                                                                                              >  2:483990/1‑151 (MQ=255)                  aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGgtgt                                                                                                                              <  2:783565/150‑1 (MQ=255)                  aTGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGtt                                                                                                                             >  1:50268/1‑151 (MQ=255)                    ggCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGttt                                                                                                                            <  2:167436/150‑1 (MQ=255)                       gcgGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGaa                                                                                                                         >  1:468532/1‑150 (MQ=255)                          gAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAAt                                                                                                                     <  2:78860/151‑1 (MQ=255)                            gCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTc                                                                                                                   >  2:316340/1‑151 (MQ=255)                            gCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTc                                                                                                                   >  2:283116/1‑151 (MQ=255)                                aaaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTg                                                                                                                >  2:141282/1‑150 (MQ=255)                                 aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGAt                                                                                                              <  2:640427/151‑1 (MQ=255)                                 aaGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGAt                                                                                                              <  1:373511/151‑1 (MQ=255)                                  aGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATg                                                                                                             >  1:524124/1‑151 (MQ=255)                                    tATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGc                                                                                                            >  1:162344/1‑150 (MQ=255)                                       tAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCaaa                                                                                                         >  1:88691/1‑150 (MQ=255)                                         aGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAAcc                                                                                                      >  1:384553/1‑151 (MQ=255)                                          gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCTATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCg                                                                                                     <  2:83400/151‑1 (MQ=255)                                          gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAAcc                                                                                                      >  1:617761/1‑150 (MQ=255)                                          gcgcGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCg                                                                                                     <  2:61467/151‑1 (MQ=255)                                            gcgAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGc                                                                                                    >  2:302190/1‑150 (MQ=255)                                                    cGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTAc                                                                                           <  2:481827/151‑1 (MQ=255)                                                             tGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCATGATGCAAAACCGCATCCCGTACAAAATATCt                                                                                  >  2:213258/1‑151 (MQ=255)                                                                ggCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggt                                                                               <  2:221864/151‑1 (MQ=255)                                                                 gCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggtg                                                                              <  2:168552/151‑1 (MQ=255)                                                                 gCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCtggtg                                                                              <  1:828495/151‑1 (MQ=255)                                                                         cgTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcgtc                                                                       <  2:577783/150‑1 (MQ=255)                                                                          gTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTAcgtcgt                                                                     <  2:775570/151‑1 (MQ=255)                                                                            tACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGtt                                                                    <  1:698293/150‑1 (MQ=255)                                                                             aCTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGttt                                                                   <  2:328253/150‑1 (MQ=255)                                                                              cTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGttttt                                                                 >  2:851940/1‑151 (MQ=255)                                                                                ggCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGttttt                                                                 <  2:598214/149‑1 (MQ=255)                                                                                ggCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTtct                                                               >  1:267269/1‑151 (MQ=255)                                                                                ggCGAGCTGAGCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGttttt                                                                 <  2:496953/149‑1 (MQ=255)                                                                                   gAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTtctctc                                                            <  2:468568/151‑1 (MQ=255)                                                                                   gAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGACGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTtctctc                                                            >  1:188840/1‑151 (MQ=255)                                                                                           cGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGa                                                     <  1:101121/150‑1 (MQ=255)                                                                                             cccATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAAt                                                  <  2:642733/151‑1 (MQ=255)                                                                                              ccATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATc                                                 >  1:633759/1‑151 (MQ=255)                                                                                                       cGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACtt                                         <  1:831946/150‑1 (MQ=255)                                                                                                       cGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACt                                          >  1:357028/1‑149 (MQ=255)                                                                                                       cGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTg                                        <  2:204294/151‑1 (MQ=255)                                                                                                         tCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCt                                      >  1:16889/1‑151 (MQ=255)                                                                                                         tCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCt                                      >  2:653203/1‑151 (MQ=255)                                                                                                          cAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTg                                     >  2:780286/1‑151 (MQ=255)                                                                                                            aTAAAACGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGc                                   >  1:737934/1‑151 (MQ=255)                                                                                                                aaCGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTAt                               <  1:660997/151‑1 (MQ=255)                                                                                                                aaCGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTAt                               <  2:342179/151‑1 (MQ=255)                                                                                                                aaCGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTAt                               >  2:689234/1‑151 (MQ=255)                                                                                                                 aCGCAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATc                              <  1:727694/151‑1 (MQ=255)                                                                                                                    cAGTACAAAGATTACGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTgc                           >  1:421289/1‑151 (MQ=255)                                                                                                                                 aCGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTa               >  1:51855/1‑150 (MQ=255)                                                                                                                                  cGCCATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAAc             <  2:465813/151‑1 (MQ=255)                                                                                                                                    ccATTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAAccc           >  1:565964/1‑151 (MQ=255)                                                                                                                                      aTTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCgg         <  2:617801/151‑1 (MQ=255)                                                                                                                                      aTTCTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCgg         >  1:503365/1‑151 (MQ=255)                                                                                                                                         cTTTATCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACg      >  2:687951/1‑151 (MQ=255)                                                                                                                                             aTCGTGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGa  <  1:617223/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                   GTTCTCCGAACTGGGTTATGGCGCGGAGCTAAAAGTATTAAGCGCGAATAACGAAGAACATGAGGCTGAGCGCGTTACTGGCGAGCTGATCGCCCATCACTTCGTCAATAAAACGCAGTACAAAGATTACGCCATTCTTTATCGCGGTAACCATCAGTCGCGGGTGTTTGAAAAATTCCTGATGCAAAACCGCATCCCGTACAAAATATCTGGTGGTACGTCGTTTTTCTCTCGTCCTGAAATCAAGGACTTGCTGGCTTATCTGCGCGTGCTGACTAACCCGGACGATGA  >  NC_000913/3961582‑3961872

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