Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,369,961	A→C	I151S (ATT→AGT)	yjeH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,369,961	0	A	C	95.5%	223.1 / 1.6	66	I151S (ATT→AGT)	yjeH	putative transporter
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base C (34/29);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 5.97e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-02

CAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGG  >  NC_000913/4369813‑4370108                                                                                                                                                     |                                                                                                                                                    cAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTa                                                                                                                                                   <  2:1026001/151‑1 (MQ=255) cAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTa                                                                                                                                                   <  2:372522/151‑1 (MQ=255)  aGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTaa                                                                                                                                                  <  2:824521/151‑1 (MQ=255)    aCCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACg                                                                                                                                                <  1:363338/151‑1 (MQ=255)     ccGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACgg                                                                                                                                               >  1:321422/1‑151 (MQ=255)       gACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGtt                                                                                                                                             <  1:331719/151‑1 (MQ=255)         cAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTg                                                                                                                                           >  1:47237/1‑151 (MQ=255)          aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGt                                                                                                                                          <  1:638348/151‑1 (MQ=255)          aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGt                                                                                                                                          >  2:207509/1‑151 (MQ=255)          aaaaCACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGt                                                                                                                                          <  1:1049786/151‑1 (MQ=255)             acacCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGa                                                                                                                                       <  1:1126504/151‑1 (MQ=255)                ccAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTa                                                                                                                                    <  2:105012/151‑1 (MQ=255)                  aGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGc                                                                                                                                  >  2:1025166/1‑151 (MQ=255)                  aGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGc                                                                                                                                  <  1:848281/151‑1 (MQ=255)                         cACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGa                                                                                                                            >  1:417182/1‑150 (MQ=255)                              aTAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGc                                                                                                                      <  2:1093933/151‑1 (MQ=255)                               tAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGc                                                                                                                      >  2:776353/1‑150 (MQ=255)                                aaCGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCAc                                                                                                                    >  1:1064338/1‑151 (MQ=255)                                    cagcaAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACcgcg                                                                                                                <  1:32599/151‑1 (MQ=255)                                        aaaTAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGt                                                                                                             >  2:330702/1‑150 (MQ=255)                                        aaaTAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATGCGCAGGTTTGATATCGCCCGCCCACCAGATGGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTa                                                                                                            <  2:193591/151‑1 (MQ=255)                                         aaTAACCCGGTAAGTTCGATATTACTAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTAc                                                                                                           >  2:373535/1‑151 (MQ=255)                                               ccGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGatat                                                                                                     <  2:625561/151‑1 (MQ=255)                                                    aaGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACcac                                                                                                <  2:1123680/151‑1 (MQ=255)                                                    aaGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACcac                                                                                                <  1:477423/151‑1 (MQ=255)                                                    aaGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACca                                                                                                 <  1:762311/150‑1 (MQ=255)                                                        tCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAcc                                                                                             >  2:1252494/1‑150 (MQ=255)                                                        tCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAcc                                                                                             >  2:1148361/1‑150 (MQ=255)                                                          gATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAccagc                                                                                          >  1:1244509/1‑151 (MQ=255)                                                             atTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAccagcgcc                                                                                       >  2:654867/1‑151 (MQ=255)                                                             atTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCCCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAccagcgcc                                                                                       >  2:1171632/1‑151 (MQ=255)                                                              ttACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACAccagcgcc                                                                                       >  1:721362/1‑150 (MQ=255)                                                                   aGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTa                                                                                 <  2:1059266/151‑1 (MQ=255)                                                                     gTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTAcc                                                                               >  1:971008/1‑151 (MQ=255)                                                                     gTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTAcc                                                                               >  2:1036327/1‑151 (MQ=255)                                                                         cGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGt                                                                           >  2:1001000/1‑151 (MQ=255)                                                                         cGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGt                                                                           >  2:746227/1‑151 (MQ=255)                                                                            aaaGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCt                                                                        <  1:1017965/151‑1 (MQ=255)                                                                            aaaGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCt                                                                        <  1:610578/151‑1 (MQ=255)                                                                                gggATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGcc                                                                     >  1:202223/1‑150 (MQ=255)                                                                                   atatTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCcaaca                                                                 <  2:821149/151‑1 (MQ=255)                                                                                   atatTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCcaac                                                                  >  1:233083/1‑150 (MQ=255)                                                                                   atatTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGCTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGc                                                                      >  1:11160/1‑146 (MQ=255)                                                                                         gCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGtt                                                           >  1:422376/1‑151 (MQ=255)                                                                                          cAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTg                                                          >  2:442842/1‑151 (MQ=255)                                                                                             gTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCa                                                       <  2:1107674/151‑1 (MQ=255)                                                                                               ttGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAg                                                      <  2:519700/150‑1 (MQ=255)                                                                                                tGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCt                                                    >  2:907036/1‑151 (MQ=255)                                                                                                  atatCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTAt                                                  >  2:1096166/1‑151 (MQ=255)                                                                                                    atCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATgc                                                >  2:646638/1‑151 (MQ=255)                                                                                                                ccaGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACAt                                     <  2:1077106/150‑1 (MQ=255)                                                                                                                ccaGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATc                                    >  2:988962/1‑151 (MQ=255)                                                                                                                 caGATAGCGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCg                                   <  2:321428/151‑1 (MQ=255)                                                                                                                        cGACAATCAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGcc                            >  2:586331/1‑151 (MQ=255)                                                                                                                               cAGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcg                     <  1:189008/151‑1 (MQ=255)                                                                                                                               cAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcg                     <  2:136409/151‑1 (MQ=255)                                                                                                                                aGCGCGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcgg                    >  2:847584/1‑151 (MQ=255)                                                                                                                                 gcgcGACGATAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcggc                   <  1:1023924/151‑1 (MQ=255)                                                                                                                                  cgcgACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCcggcg                  <  2:612147/151‑1 (MQ=255)                                                                                                                                       cgaTAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAAt             >  2:621010/1‑151 (MQ=255)                                                                                                                                        gaTAAGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAAtt            >  1:294297/1‑151 (MQ=255)                                                                                                                                           aaGTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCGATTTg          <  1:1035659/150‑1 (MQ=255)                                                                                                                                             gTCCGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTg          >  1:1006949/1‑148 (MQ=255)                                                                                                                                               ccGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTg     >  2:943744/1‑151 (MQ=255)                                                                                                                                                cGGCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGc    >  2:668367/1‑151 (MQ=255)                                                                                                                                                  gCACTAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCgg  <  1:998406/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAGACCGACAAAACACCAGAACATCACTGATAACGCAGCAAATAACCCGGTAAGTTCGATATTACCAGGTGCCGGAAAGGGGATATTCGCAGGTTTGATATCGCCCGCCCACCAGATAGCGACAATCAGCGCGACGATAAGTCCGGCAATAACGGTTTGTAGATTAGCACTGGAACTGGCACCGCGAGTACCGATATACCACACCAGCGCCAGCGTACCGAGTTCTGCCAACAACAGTTGCCAGCTATGCCAGCCAAACATCGCCTGGCCGAACCCGGCGGCAATTTGTAGTGCGG  >  NC_000913/4369813‑4370108

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