Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,086,469	C→T	6.0%	A198A (GCG→GCA)	plaP ←	putrescine importer, low affinity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,086,469	0	C	T	6.0%	220.5 / 2.8	67	A198A (GCG→GCA)	plaP	putrescine importer, low affinity
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base T (2/2);  total (30/37)
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.11e-01

GCCGCCGATCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATT  >  NC_000913/2086325‑2086614                                                                                                                                                 |                                                                                                                                                  gccgccGATCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCa                                                                                                                                                >  3:132905/1‑148 (MQ=255)      cgATCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGcc                                                                                                                                           <  1:153132/148‑1 (MQ=255)         tCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGTGCATCGccaga                                                                                                                                        >  4:183902/1‑148 (MQ=255)         tCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccaga                                                                                                                                        >  2:51753/1‑148 (MQ=255)              gcgGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGTGCATCGccagaccaga                                                                                                                                   <  3:183902/148‑1 (MQ=255)              gcgGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccagaccaga                                                                                                                                   <  4:132905/148‑1 (MQ=255)                  tCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATgg                                                                                                                               <  2:4913/148‑1 (MQ=255)                  tCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATgg                                                                                                                               >  4:155036/1‑148 (MQ=255)                  tCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATgg                                                                                                                               <  3:312130/148‑1 (MQ=255)                    aGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGa                                                                                                                              <  2:360778/147‑1 (MQ=255)                    aGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGAc                                                                                                                             <  1:40536/148‑1 (MQ=255)                     ggAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGAc                                                                                                                             >  3:265873/1‑147 (MQ=255)                     ggAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACg                                                                                                                            <  8:235415/148‑1 (MQ=255)                          aTCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGc                                                                                                                       <  5:367199/148‑1 (MQ=255)                            cGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGcc                                                                                                                                           >  6:275769/1‑126 (MQ=255)                            cGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGcc                                                                                                                                           <  5:275769/126‑1 (MQ=255)                              cACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGc                                                                                                                   >  3:188570/1‑148 (MQ=255)                                       tCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGcc                                                                                                                                           >  2:225177/1‑115 (MQ=255)                                       tCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGcc                                                                                                                                           <  1:225177/115‑1 (MQ=255)                                               tctGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGcctt                                                                                                  >  8:160757/1‑148 (MQ=255)                                                ctGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGcctt                                                                                                  >  3:124112/1‑147 (MQ=255)                                                   cATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGccttcac                                                                                               <  4:265873/147‑1 (MQ=255)                                                   cATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGccttca                                                                                                <  1:51753/146‑1 (MQ=255)                                                   cATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGccttca                                                                                                <  4:234195/146‑1 (MQ=255)                                                      cTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGccttcacctt                                                                                            >  7:241043/1‑147 (MQ=255)                                                         tGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACATGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGccttcacctt                                                                                            <  2:140133/144‑1 (MQ=255)                                                         tGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCaa                                                                                         <  3:155036/147‑1 (MQ=255)                                                         tGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCaa                                                                                         <  8:241043/147‑1 (MQ=255)                                                             ttcttcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccagaccaga                                                                                                                                   <  2:54374/101‑1 (MQ=255)                                                             ttcttcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccagaccaga                                                                                                                                   >  1:54374/1‑101 (MQ=255)                                                             ttcttcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccagaccaga                                                                                                                                   >  5:205148/1‑101 (MQ=255)                                                             ttcttcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGccagaccaga                                                                                                                                   <  6:205148/101‑1 (MQ=255)                                                             ttcttcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACg                                                                                    >  2:83887/1‑148 (MQ=255)                                                                 tcCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCa                                                                                 >  8:290374/1‑147 (MQ=255)                                                                         gTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAAcc                                                                        <  4:188570/148‑1 (MQ=255)                                                                          ttGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCa                                                                       >  6:206278/1‑148 (MQ=255)                                                                               gATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGcc                                                                   >  6:326507/1‑147 (MQ=255)                                                                                  gCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACg                                                                                                                            <  1:286447/87‑1 (MQ=255)                                                                                  gCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACg                                                                                                                            >  2:286447/1‑87 (MQ=255)                                                                                       cAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       >  1:36333/1‑119 (MQ=255)                                                                                       cAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       <  2:36333/119‑1 (MQ=255)                                                                                       cAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCAc                                                          >  4:15280/1‑148 (MQ=255)                                                                                         aaGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCCGCGTACCAGCGCCTTCACCTTCAACTACGCCATAACCCACCCTGCCCAGAATCACcg                                                        >  4:246109/1‑148 (MQ=255)                                                                                                tAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAg                                                 <  5:206278/148‑1 (MQ=255)                                                                                                  aaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGTGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGc                                                <  2:297288/147‑1 (MQ=255)                                                                                                  aaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       <  1:287936/108‑1 (MQ=255)                                                                                                  aaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       >  2:287936/1‑108 (MQ=255)                                                                                                  aaGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGc                                                <  7:160757/147‑1 (MQ=255)                                                                                                     gagaAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCaaa                                                                                        <  6:225401/104‑1 (MQ=255)                                                                                                     gagaAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCaaa                                                                                        >  5:225401/1‑104 (MQ=255)                                                                                                       gaAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATc                                                   >  2:295134/1‑139 (MQ=255)                                                                                                         aGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGTGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACg                                                                                    >  7:154687/1‑104 (MQ=255)                                                                                                         aGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGTGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACg                                                                                    <  8:154687/104‑1 (MQ=255)                                                                                                           cacaGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGc                                                <  5:48584/138‑1 (MQ=255)                                                                                                           cacaGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGc                                                >  6:48584/1‑138 (MQ=255)                                                                                                             caGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAac                                    <  3:15280/148‑1 (MQ=255)                                                                                                                         cccGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACa                        <  1:83887/148‑1 (MQ=255)                                                                                                                               aaTCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       <  5:42745/79‑1 (MQ=255)                                                                                                                               aaTCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       >  6:42745/1‑79 (MQ=255)                                                                                                                                atcatcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGa                 <  3:246109/148‑1 (MQ=255)                                                                                                                                  catcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       >  3:53778/1‑76 (MQ=255)                                                                                                                                  catcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAAt                                                                                       <  4:53778/76‑1 (MQ=255)                                                                                                                                    tcGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCACGATCAGc                                                <  7:204221/113‑1 (MQ=255)                                                                                                                                        gATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCt         <  7:290374/148‑1 (MQ=255)                                                                                                                                             cATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCt         <  3:339903/143‑1 (MQ=255)                                                                                                                                             cATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCt         >  4:339903/1‑143 (MQ=255)                                                                                                                                               tGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGAtt  >  7:321147/1‑148 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GCCGCCGATCAGCGCGGTCAGGAAAATCGCACGCGGGATCACGCGCTCTGCATCTTTGGTTTCTTCCGACAGGTTGCTGATGCCGTCAAAGCCGGTAAAGGAGAAGCACAGGATTGTCGCCCCGGTAATCATCGGGATAACATGCGCATCGCCAGACCAGAATGGACGAGTGCTCGCCAGCGTACCAGCGCCTTCACCTTCAAATACGCCATAAACAACCATGCCCAGAATCACCGCGATCAGCACTACCTGCAACACGACGATTACAGTGTTGAAGTTCGCTACGGATT  >  NC_000913/2086325‑2086614

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