Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,559,306	A→C	I202S (ATC→AGC)	ddpC ←	D,D‑dipeptide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,559,306	0	A	C	100.0%	231.8 / NA	66	I202S (ATC→AGC)	ddpC	D,D‑dipeptide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (38/28);  total (38/28)

ACCATCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAA  >  NC_000913/1559162‑1559446                                                                                                                                                 |                                                                                                                                             aCCATCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAg                                                                                                                                          <  1:547248/149‑1 (MQ=255)  ccATCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATg                                                                                                                                       <  2:730040/151‑1 (MQ=255)   cATCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGGAGGGAGTTATGTTAAAAAAGCCAGGTGCTCAGAGGc                                                                                                                                      >  1:589865/1‑151 (MQ=255)     tCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGcc                                                                                                                                     >  2:751237/1‑150 (MQ=255)     tCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTACGATGcc                                                                                                                                     >  1:986293/1‑150 (MQ=255)      cgcgCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGGGTTACGTAAAATATGACAGTTGCTCAAATGCCa                                                                                                                                    >  1:442642/1‑150 (MQ=255)      cgcgCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCa                                                                                                                                    >  2:474932/1‑150 (MQ=255)       gcgcCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCCGATGCCaa                                                                                                                                   >  1:367808/1‑150 (MQ=255)       gcgcCCCATTCAGCACGCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCa                                                                                                                                    <  1:322863/149‑1 (MQ=255)        cgcCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAAc                                                                                                                                  >  1:771009/1‑150 (MQ=255)        cgcCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAAc                                                                                                                                  <  1:307654/150‑1 (MQ=255)         gcCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACg                                                                                                                                 >  1:567729/1‑150 (MQ=255)          ccccATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACgag                                                                                                                               <  1:127667/151‑1 (MQ=255)           cccATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACg                                                                                                                                 >  1:529609/1‑148 (MQ=255)                 cAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACgcgc                                                                                                                         <  2:162362/150‑1 (MQ=255)                  aGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCa                                                                                                                       <  1:520624/151‑1 (MQ=255)                    cACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCaaa                                                                                                                     <  1:393850/151‑1 (MQ=255)                      cTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAgg                                                                                                                   >  2:623372/1‑151 (MQ=255)                        cGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGtt                                                                                                                 >  1:651716/1‑151 (MQ=255)                          gttgttGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGtttt                                                                                                               >  2:719843/1‑151 (MQ=255)                                    aCCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGACACGCGCCAAAGGGTTTCGCCGCCTg                                                                                                      >  1:480337/1‑150 (MQ=255)                                      ccAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTCGGGAGGTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGCTTTCGCCGCCt                                                                                                       >  1:238347/1‑147 (MQ=255)                                        aGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACAt                                                                                                 <  2:842613/151‑1 (MQ=255)                                        aGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACAt                                                                                                 >  2:630571/1‑151 (MQ=255)                                          gCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAg                                                                                               <  2:368695/151‑1 (MQ=255)                                            cAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGATTTCGCCGCCTGAACATAgg                                                                                              >  1:829904/1‑150 (MQ=255)                                                  aacccgaagTGGCGGCCATTAAAAACGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGGTACGTAAAATATGCCCGTTGCTCAGATTCCAACGAGACCCGCCAAAGGGTTTCGCCCCCTGAACATAggggtac                                                                                         >  1:636197/9‑144 (MQ=255)                                                  aTCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTg                                                                                       <  1:495266/151‑1 (MQ=255)                                                         cGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTAc                                                                                 >  2:672651/1‑150 (MQ=255)                                                           tGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAAc                                                                              >  2:1263648/1‑151 (MQ=255)                                                               ggcCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTaa                                                                           >  1:410560/1‑150 (MQ=255)                                                                 ccATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGc                                                                        <  1:909085/151‑1 (MQ=255)                                                                      aaaaTCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTAAGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGcc                                                                   >  2:49185/1‑151 (MQ=255)                                                                             cgcTACCGATATCCAGCGATGCCTGCCCGATCAGGGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAAAATAGGTATACTGGCGTGCAACTAATGCTTGCCCCcgcgt                                                            <  1:619155/151‑2 (MQ=255)                                                                             cgcTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCgcgc                                                            <  2:595188/151‑1 (MQ=255)                                                                             cgcTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCgcgc                                                            >  2:945007/1‑151 (MQ=255)                                                                              gcTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGcc                                                           >  1:1037496/1‑151 (MQ=255)                                                                                tACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGcc                                                           >  2:733634/1‑149 (MQ=255)                                                                                       atCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGcac                                                   <  2:604151/150‑1 (MQ=255)                                                                                       atCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGcac                                                   >  2:438875/1‑150 (MQ=255)                                                                                            gCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATaaaa                                             <  2:309493/151‑1 (MQ=255)                                                                                            gCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGAAGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATaaa                                              >  2:1264637/1‑150 (MQ=255)                                                                                             cGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATaaa                                              >  2:1099556/1‑149 (MQ=255)                                                                                                tGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAggg                                          >  2:563700/1‑150 (MQ=255)                                                                                                tGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAggg                                          >  1:1187523/1‑150 (MQ=255)                                                                                                  ccTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAAt                                       <  2:892354/151‑1 (MQ=255)                                                                                                    tGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTc                                     >  2:1081130/1‑151 (MQ=255)                                                                                                       aCGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGca                                  <  1:926954/151‑1 (MQ=255)                                                                                                         gATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGcaca                                <  1:865073/151‑1 (MQ=255)                                                                                                            cAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaata                             <  1:1265990/151‑1 (MQ=255)                                                                                                             aGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaatag                            <  1:671712/151‑1 (MQ=255)                                                                                                               gggcggAATGGAGTTACGTAGAATATTCCAGTTGCGCAGACGCCAACGAGACGCCACAAAGGTTTTCGCCGCCTGAACGTTGGGATACTGGCGTACAACTAATGCTTGGCACCCCGCCAGGCGCACATAAAAGGGAATTCGCCcaatag                            >  1:1245837/2‑149 (MQ=255)                                                                                                                  cggTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaatagcaat                        >  1:431767/1‑150 (MQ=255)                                                                                                                   ggTAGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCAcaatagcaatag                      <  1:751203/151‑1 (MQ=255)                                                                                                                      aGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGcc                    >  2:94164/1‑150 (MQ=255)                                                                                                                      aGGGAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGcc                    >  2:774007/1‑150 (MQ=255)                                                                                                                         gAGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGc                 >  1:396081/1‑150 (MQ=255)                                                                                                                          aGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATATGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCAt               <  1:832859/151‑1 (MQ=255)                                                                                                                          aGTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCAt               <  2:994647/151‑1 (MQ=255)                                                                                                                           gTTACGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATg              <  1:78128/151‑1 (MQ=255)                                                                                                                              aCGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCg           >  1:1229340/1‑151 (MQ=255)                                                                                                                              aCGTAAAATATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTGCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGc            <  1:519307/150‑1 (MQ=255)                                                                                                                                  aaaaTATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTaa       <  2:651752/151‑1 (MQ=255)                                                                                                                                   aaaTATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTaa       >  2:363937/1‑150 (MQ=255)                                                                                                                                    aaTATGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTaaac     >  2:1066722/1‑151 (MQ=255)                                                                                                                                       atGCCAGTTGCTCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTaaacaaa  <  2:897042/151‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                             ACCATCGCGCCCCATTCAGCACTCGGTTGTTGAGCACCCAGGCCAATAAATCCCAACGTGGCGGCCATTAAAATCGCGCTACCGATATCCAGCGATGCCTGCACGATCAGCGGCGGTAGGGAGTTACGTAAAATATGCCAGTTGATCAGATGCCAACGAGACGCGCCAAAGGTTTTCGCCGCCTGAACATAGGTATACTGGCGTACAACTAATGCTTGCCCCCGCGCCAGGCGCACATAAAAGGGAATTCGCACAATAGCAATAGCCAGCATGGCGTTAAACAAA  >  NC_000913/1559162‑1559446

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