Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,377,792	A→T	6.8%	intergenic (‑70/+19)	blc ← / ← ampC	outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,377,792	0	A	T	6.8%	199.7 / 2.9	72	intergenic (‑70/+19)	blc/ampC	outer membrane lipoprotein cell division and growth lipocalin/penicillin‑binding protein; beta‑lactamase, intrinsically weak
Reads supporting (aligned to +/- strand):  ref base A (22/45);  new base T (2/3);  total (24/48)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

CCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCG  >  NC_000913/4377643‑4377937                                                                                                                                                      |                                                                                                                                                  ccacgcGGCGGCGTAGGAGCACTGCAGGCAACGACCAGAAATGCCGCTGTCGCTGCGGCACCGAGTGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGaaaa                                                                                                                                                   >  2:1266502/1‑150 (MQ=255)     gcggcggcgTAGGAGAACTGCAGGCAACGCCCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTAACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGATAATTCg                                                                                                                                               <  1:235668/150‑1 (MQ=255)     gcggcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCg                                                                                                                                               <  2:321899/150‑1 (MQ=255)        gcggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAc                                                                                                                                            <  2:1203724/150‑1 (MQ=255)         cggcgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAcc                                                                                                                                           <  1:462281/150‑1 (MQ=255)          gtcgTAGGAGAAGTGCAGGCAACGACCACAAATCCAGCTGTCGCTGCGGCAACGAGATGGAGCAGGCGCAAAAATGTTCCCTTACTGCTTATTTTCAAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                                                                          <  2:223124/148‑1 (MQ=255)          ggcgTAGGAGAACTGCACGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGCGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGGGTTTTTCTAAGTGTAGATGACAGCAAGCAAAAGCGGAGAAAAGGTCCGAAAATTCGGAccc                                                                                                                                          <  2:108585/150‑1 (MQ=255)            ggTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCGTCCTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttt                                                                                                                                        <  1:535347/149‑4 (MQ=255)            cgTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTATTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGa                                                                                                                                        >  2:1031565/1‑150 (MQ=255)             gTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAACGGTCCGAAAATTCGGACCCGAt                                                                                                                                       >  1:400023/1‑150 (MQ=255)              tAGGAGAACTTCAGGCAACGACCCTAAATGCAGCTGTCGCTGCGGCAACGAGACGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATg                                                                                                                                      <  2:352241/150‑1 (MQ=255)              tAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATg                                                                                                                                      <  1:1255225/150‑1 (MQ=255)              aaGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGGCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAATTTTCGGAccttttc                                                                                                                                      <  2:1021715/149‑6 (MQ=255)               aGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGCGAGGGAGCAGGCGCATAAATGGTTCCTTACTGGTTGTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGGCCGAAAATTCGGACCCGATgg                                                                                                                                     <  2:160999/150‑1 (MQ=255)               aGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATgg                                                                                                                                     <  1:229164/150‑1 (MQ=255)                 gagaACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGaa                                                                                                                                   >  2:556772/1‑150 (MQ=255)                  agaACTGCAGGCCACGACCATAAAAGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAACTGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAt                                                                                                                                  <  1:650249/150‑1 (MQ=255)                    aaCTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAAttt                                                                                                                                <  1:268003/150‑1 (MQ=255)                        gCAGTCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACt                                                                                                                            <  2:612432/150‑1 (MQ=255)                            gCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGAATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGTCAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTag                                                                                                                        <  1:237513/150‑1 (MQ=255)                             cAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTaga                                                                                                                       >  1:1170320/1‑150 (MQ=255)                              aaCGACCAGAAATGCAGCTGTCGCGGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTagag                                                                                                                      <  2:295501/150‑1 (MQ=255)                               aCGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGc                                                                                                                     >  2:1187899/1‑150 (MQ=255)                                 gACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGAATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGt                                                                                                                   <  1:1018031/150‑1 (MQ=255)                                    cAGAAATGCAGCTGTCGCTGCGGCAACGCGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTaa                                                                                                                <  2:1163530/150‑1 (MQ=255)                                         atggcgCTCTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCTAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCTATGGAATTTTACTGTAGAGCGTTAAGCATc                                                                                                           <  2:600194/145‑1 (MQ=255)                                           gCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTg                                                                                                         >  2:364810/1‑150 (MQ=255)                                               cTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCACGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTACAGCGTTAAGAATCTGCCAg                                                                                                     <  2:735634/150‑1 (MQ=255)                                                     cTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGATAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGgcggcg                                                                                               <  2:352955/150‑1 (MQ=255)                                                           cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                                                                         >  1:516657/1‑150 (MQ=255)                                                           cAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAct                                                                                         >  2:860532/1‑150 (MQ=255)                                                           cAACGAGAGGGAGCAGGCGCATAAATGTTACCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGA‑CCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGActc                                                                                        <  2:38127/150‑1 (MQ=255)                                                            aaCGTGAGGGAGCAGGCGCATAACTGTGTCCTTACTGGTTTTTTTCTAAGTGTAGATCACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGAgtc                                                                                        <  2:509125/150‑1 (MQ=255)                                                               gagaGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCg                                                                                     <  2:702598/150‑1 (MQ=255)                                                                 gagGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGGCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCt                                                                                   <  2:359251/150‑1 (MQ=255)                                                                 gagGGAGCAGGCGCATAAATGTTTCCTTAATGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCt                                                                                   <  1:368298/150‑1 (MQ=255)                                                                        cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTCGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                                                                            <  1:1231857/150‑1 (MQ=255)                                                                        cAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTgg                                                                            <  2:516657/150‑1 (MQ=255)                                                                           gcgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGAt                                                                         <  2:299010/150‑1 (MQ=255)                                                                            cgcATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATa                                                                        <  2:603457/150‑1 (MQ=255)                                                                              cATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGt                                                                      >  2:1041458/1‑150 (MQ=255)                                                                                       ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                                             >  2:1098553/1‑150 (MQ=255)                                                                                       ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                                             <  2:170571/150‑1 (MQ=255)                                                                                       ttCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                                             <  1:654173/150‑1 (MQ=255)                                                                                       taccTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTtttgtttg                                                             <  1:36897/148‑1 (MQ=255)                                                                                        tCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGc                                                            <  1:495284/150‑1 (MQ=255)                                                                                         ccTAACTGGTGTTTTTCTAAGTGGAGATGACAGCATGGAACAGCGGTGAAAAGGTCCGATAATTAGTACCCGATGGAATTTTACTGTATAGCGTTAATAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGcc                                                           <  2:63013/150‑1 (MQ=255)                                                                                            tACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGc                                                        <  1:136481/150‑1 (MQ=255)                                                                                                ggTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCa                                                    >  2:296987/1‑150 (MQ=255)                                                                                                 gtatttttCTAAGTGTAGATGACAGCAAGGTAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCAc                                                   >  1:779088/4‑150 (MQ=255)                                                                                                 gTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCAc                                                   <  1:1047889/150‑1 (MQ=255)                                                                                                      tttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAc                                              <  1:1193502/150‑1 (MQ=255)                                                                                                       ttCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCCGCATCACGATAcc                                             >  2:1275048/1‑150 (MQ=255)                                                                                                         cTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCa                                           >  1:453309/1‑150 (MQ=255)                                                                                                             gtgtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGctc                                       <  1:1029533/150‑1 (MQ=255)                                                                                                               gtAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCTTCGACTCTCGCTGGATTGGGAGAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCtt                                     >  1:925639/1‑150 (MQ=255)                                                                                                                tAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCttt                                    >  1:1255327/1‑150 (MQ=255)                                                                                                                   aTGACAGCAAGTAAAAGCTGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTc                                 <  1:720449/150‑1 (MQ=255)                                                                                                                    tGACAGCAAGGACAACCGCACACAACTGCCTAAAATTCGGACCCTAAGGAATTTTACTGTAGAGCTTTAAGAATCTGCCACGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCt                                <  2:992456/150‑1 (MQ=255)                                                                                                                    tGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCt                                <  1:1206376/150‑1 (MQ=255)                                                                                                                     gACAGCAAGGAAAAGCGGAGAAAAGGTCCTAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTg                               <  2:459199/150‑1 (MQ=255)                                                                                                                     gACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCAAGATACCCAGCTCTTTTTCtt                               >  2:240287/1‑149 (MQ=255)                                                                                                                       cAGCAAGGAAAAGCGGCGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGa                             <  1:375100/150‑1 (MQ=255)                                                                                                                               aaaaGCGGAGAAAAGGTCCGAAAATTCGGACCTGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATTGTTTTTGTTTGCCAGCATCACGCTACCCAGCTCTTTTTCGGGAATAAAcgc                     >  2:1167728/1‑150 (MQ=255)                                                                                                                                    cGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGAGTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACAt                >  1:1157661/1‑150 (MQ=255)                                                                                                                                        gaAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCGGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCt            >  2:777693/1‑150 (MQ=255)                                                                                                                                        catcgggTCCGAATTTTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCt            >  2:1214607/6‑150 (MQ=255)                                                                                                                                        caaatGGCCCCAAAATGCGCACCCGTTGTAATTTTACTGTCGAGCGTTAATTATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGATTGCCACCATCACGATACCCAGCTCTTTTTCTGGAAAAAACGCGACATAGCt            <  2:985128/149‑1 (MQ=255)                                                                                                                                             ggTCCGAATTTTCGGACCTGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCaa       >  2:393796/1‑150 (MQ=255)                                                                                                                                                 cGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATAACCAGCTCTTTTTCTGGAATAAACGCGACATAGCTAACAAAtac   >  1:344059/1‑148 (MQ=255)                                                                                                                                                  taaaaTCCGGACCCGATGGAATTTTACTGTAGCGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccg  <  1:1266502/149‑1 (MQ=255)                                                                                                                                                  gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGTGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGATTAAACGCGACATAGCTACCAAATccg  >  1:131684/1‑150 (MQ=255)                                                                                                                                                  gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccg  <  1:419010/150‑1 (MQ=255)                                                                                                                                                  gAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATccg  <  2:388857/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  CCACGCGGCGGCGTAGGAGAACTGCAGGCAACGACCAGAAATGCAGCTGTCGCTGCGGCAACGAGAGGGAGCAGGCGCATAAATGTTTCCTTACTGGTTTTTTTCTAAGTGTAGATGACAGCAAGGAAAAGCGGAGAAAAGGTCCGAAAATTCGGACCCGATGGAATTTTACTGTAGAGCGTTAAGAATCTGCCAGGCGGCGTCGACTCTCGCTGGATTGGGATAGTTTTTGTTTGCCAGCATCACGATACCCAGCTCTTTTTCTGGAATAAACGCGACATAGCTACCAAATCCG  >  NC_000913/4377643‑4377937

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