Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,004,063	A→C	L39R (CTG→CGG)	rarD ←	putative chloramphenical resistance permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,004,063	0	A	C	100.0%	232.6 / NA	63	L39R (CTG→CGG)	rarD	putative chloramphenical resistance permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (27/36);  total (27/36)

TAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAA  >  NC_000913/4003915‑4004204                                                                                                                                                     |                                                                                                                                               tAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCgg                                                                                                                                              >  2:225376/1‑151 (MQ=255)  aGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCgg                                                                                                                                              <  2:377565/150‑1 (MQ=255)  aGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCgg                                                                                                                                              >  2:476375/1‑150 (MQ=255)     cAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGAttt                                                                                                                                          >  1:231500/1‑151 (MQ=255)       gTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCa                                                                                                                                        <  2:582256/151‑1 (MQ=255)          gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATcg                                                                                                                                     >  2:595331/1‑151 (MQ=255)          gCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATc                                                                                                                                      <  2:710641/150‑1 (MQ=255)                aaTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATcggcgg                                                                                                                                 >  1:387144/1‑149 (MQ=255)                aaTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGc                                                                                                                               <  1:323612/151‑1 (MQ=255)                 aTCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCgggg                                                                                                                               >  2:136561/1‑149 (MQ=255)                   cAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCAc                                                                                                                             <  1:428766/150‑1 (MQ=255)                    aGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACgt                                                                                                                           <  1:217762/151‑1 (MQ=255)                       aCGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACgta                                                                                                                          >  1:728749/1‑149 (MQ=255)                            agagACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTaaa                                                                                                                     >  2:920859/1‑149 (MQ=255)                             gagaCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCa                                                                                                                  <  2:413874/151‑1 (MQ=255)                             gagaCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCa                                                                                                                  <  2:405389/151‑1 (MQ=255)                                aCTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACt                                                                                                               <  2:741077/151‑1 (MQ=255)                                     cAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAg                                                                                                          <  2:491883/151‑1 (MQ=255)                                      aGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGGAGTAAAACAAATTTAAg                                                                                                          >  1:761183/1‑150 (MQ=255)                                        cATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTa                                                                                                        >  2:616316/1‑150 (MQ=255)                                          tAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAAACAACTTTAAGTa                                                                                                        >  1:738061/1‑148 (MQ=255)                                           aaaaaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCt                                                                                                    <  1:341078/151‑1 (MQ=255)                                              aaTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgc                                                                                                 <  1:570778/151‑1 (MQ=255)                                               aTTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGgcg                                                                                                <  1:51046/151‑1 (MQ=255)                                                       tgGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGGGCGTGGTAAATATACTTGTAGCACgacggccccggaacc                                                                                        >  2:256325/1‑137 (MQ=255)                                                                 aTCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGGACGTAGGAAATCAAATTGAAGTACGCTGGCGCTAT‑AACCCaaataaaa                                                                              >  1:764637/1‑150 (MQ=255)                                                                   cAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAaataaaataa                                                                           <  1:406000/151‑1 (MQ=255)                                                                       gTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCgg                                                                       <  2:358607/151‑1 (MQ=255)                                                                           ttAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCaa                                                                    <  1:675884/150‑1 (MQ=255)                                                                            tAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCCCGTAGGAAATCCACATGAAGTACGCCTGCGCCAT‑ACCCCAAATAAAATAAACGGCaa                                                                    >  1:770317/1‑149 (MQ=255)                                                                             aaaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCATGTAGTAAATAAACTTGAAGGAAGCCTGCGCCCTAACCCCAAAAAAAAAAACCGCCa                                                                     >  1:906706/1‑148 (MQ=255)                                                                             aaaTAGGACCACTGGCGGCAAATGCTCAGCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAaga                                                                  <  1:836951/150‑1 (MQ=255)                                                                              aaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGc                                                                <  1:126210/151‑1 (MQ=255)                                                                              aaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTAAGCTGGCGCTAT‑ACCCCCAATAAAATAAACGGCaa                                                                    >  2:857927/1‑147 (MQ=255)                                                                              aaTAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGAAAGAAATCAAACTGAAATAAGCTGGCGCCAT‑ACCCCAAATAAAAAAAAAGGCaa                                                                    >  2:776483/1‑147 (MQ=255)                                                                                 aGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGGACCTAGTAAATCAACTTGAAAGAAGCCGGCGGTAT‑ACCCCCAATAAAAAAAGACGCaaaa                                                                  >  2:662962/1‑144 (MQ=255)                                                                                      cACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGGAGTAAAACAACATGAAGTACGCCGGGGCTAT‑ACCCCAAAAAAAAAAAACGGCAAGAGCGAgga                                                           >  1:107692/1‑146 (MQ=255)                                                                                        cTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCCAATAAAAAAAAAGGCAAAAGCGAGGAATa                                                        >  2:48022/1‑149 (MQ=255)                                                                                        cTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGGAGTAAATCAAATTGAAATAAGCCGGCGCTAT‑ACCCCAAATAAAAAAAAAGGCAAGAGCGAGTaaaa                                                        >  1:362525/1‑147 (MQ=255)                                                                                          ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                                     <  1:774401/150‑1 (MQ=255)                                                                                          ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGc                                                     <  1:622345/150‑1 (MQ=255)                                                                                          ggcggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAACGGCAAGAGCGAGTAATACGcc                                                    >  2:265832/1‑151 (MQ=255)                                                                                            cggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                                  <  1:232954/151‑1 (MQ=255)                                                                                            cggcAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCt                                                  <  2:582092/151‑1 (MQ=255)                                                                                                     cTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtt                                          <  2:715381/150‑1 (MQ=255)                                                                                                       catcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgttt                                         <  1:199075/149‑1 (MQ=255)                                                                                                        atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCgtttg                                        <  1:142639/149‑1 (MQ=255)                                                                                                        atcaGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGGAGTAAATCAACTTGAAGTAAGCTGGCGCTAT‑ACCCCCAATAAAATAAACGGCAAGAGCGAGTAATACGCCCTGCCGCgtttgtt                                      >  2:326655/1‑151 (MQ=255)                                                                                                            gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGc                                  <  1:828505/151‑1 (MQ=255)                                                                                                            gCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATAAACTTTAAGTACGCTGGCGCTAT‑ACCCCAAATAAAAAAAACGGCAAGAGCGAGTAATACGCCCTGCCGCGGTTGtttt                                    >  2:239805/1‑149 (MQ=255)                                                                                                               ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTTAAGGAAGCCGGGGCTAT‑ACCCCCAATAAAAAAAACGGCAAGAGCGAGGAATACGCCCCGCCCCGTTTGTTTTTca                                 >  2:770019/1‑149 (MQ=255)                                                                                                               ccaccaTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatc                               <  2:381975/151‑1 (MQ=255)                                                                                                                  ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                             <  2:144172/150‑1 (MQ=255)                                                                                                                  ccaTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGcatcca                             <  1:505189/150‑1 (MQ=255)                                                                                                                    aTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg                          <  2:813861/151‑1 (MQ=255)                                                                                                                    aTAAAGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCg                          <  2:595123/151‑1 (MQ=255)                                                                                                                       aaGAAAAACGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGa                       <  1:379866/151‑1 (MQ=255)                                                                                                                          aaaaaCGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCAAGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCGCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATa                    <  1:417322/151‑1 (MQ=255)                                                                                                                            aaaCGACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACt                  <  2:799010/151‑1 (MQ=255)                                                                                                                                gACCAGATCACGCGATGCGTCCGGATTTCATCGGCGGGCCCGTAGTAAAACAACTTGAAGTAAGCTGGCGCGAT‑ACCCCAAATAAAATAAACGGCAAGAGCGAGTAATAAGCCCCGCCTCGGTTGTTTTTCATACCTCCGGAAAACTCAtt              >  2:926800/1‑151 (MQ=255)                                                                                                                                       tCACGCGATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAAAGCGAGTAATACGCCCTGCCGCGTTTGTTTTTCATCCCTCGGGGATACTCATTTTTTAtt       >  1:463438/1‑151 (MQ=255)                                                                                                                                           gcgATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTa   <  1:575838/151‑1 (MQ=255)                                                                                                                                            cgATGCGTCCGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTaa  <  2:738079/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               TAGCCAGTTGCCACCAATCAGCACGGCAGAGACTGCCAGCATAAAAATTTTCTGTGGCGTCTGAATCAGCGTTTTTAAATAGGACCACTGGCGGCAAATGCTCATCAGCACCACCATAAAGAAAAACGACCAGATCACGCGATGCGTCAGGATTTCATCGGCGGGCACGTAGTAAATCAACTTGAAGTACGCTGGCGCTAT‑ACCCCAAATAAAATAAGCGGCAAGAGCGAGTAATACGCCCTGCCGCGTTTGTTTTGCATCCATCGGGAATACTCATTTTTAATTGGTAA  >  NC_000913/4003915‑4004204

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