Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,019,927	C→T	6.3%	G43S (GGT→AGT)	ompA ←	outer membrane protein A (3a;II*;G;d)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,019,927	0	C	T	6.3%	203.3 / 3.4	63	G43S (GGT→AGT)	ompA	outer membrane protein A (3a;II*;G;d)
Reads supporting (aligned to +/- strand):  ref base C (28/31);  new base T (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

GTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTAT  >  NC_000913/1019784‑1020070                                                                                                                                                |                                                                                                                                                gTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAgtg                                                                                                                                             >  5:258399/1‑148 (MQ=255)        aCGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                                                                                                                                               <  2:217394/139‑1 (MQ=255)        aCGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAg                                                                                                                                               >  1:217394/1‑139 (MQ=255)              ccTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACTAGTGTCATGGTACTggg                                                                                                                                >  8:79239/1‑148 (MQ=255)              ccTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTggg                                                                                                                                <  6:255097/148‑1 (MQ=255)              ccTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTggg                                                                                                                                <  7:260502/148‑1 (MQ=255)                 ttGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAAACACCAAAAGCCCCAGCGCCCAGGTGGTTATCATGGGTCGGGCCAGTGTTGGTGATGAAACCAGTGTCATGGTACTGGGAcc                                                                                                                             >  1:193219/1‑148 (MQ=255)                    tACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATTCGGGTTAACCTGGTCACCACCAAAAGCACCCGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGCTGAAACCAGTGGCATGGTACTGGGACCAg                                                                                                                           <  5:47006/147‑1 (MQ=255)                     aCGGCGTACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGcc                                                                                                                         >  8:301750/1‑148 (MQ=255)                      cGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGcc                                                                                                                         <  4:30731/147‑1 (MQ=255)                            aCGACCTAACCAGTCGTAACCCATTTCAAAGCCATCATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGTGCCATTGTTGTTGATGAAAcc                                                                                                                                                 <  2:72275/117‑1 (MQ=255)                            aCGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAAcc                                                                                                                                                 >  1:72275/1‑117 (MQ=255)                               aCCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCCCCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCa                                                                                                               >  6:287040/1‑148 (MQ=255)                                 cTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCAcc                                                                                                             <  2:203718/148‑1 (MQ=255)                                    ccgaGTGGTAACCCATTTCAATGCAACCTTACGGGTTAACGTGTTAACCACAAAAAGCACCAGGGCCCTGTTGGTTTTCATGGGT‑GGGCCATTGTTGTTGATGATACCTGTGTCTTGGTACTGGGACCAGCGCTGTTTAGCACCAgtg                                                                                                         <  8:181660/145‑1 (MQ=255)                                     ccAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACTAgtgt                                                                                                                                            >  2:156192/1‑113 (MQ=255)                                     ccAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACTAgtgt                                                                                                                                            <  1:156192/113‑1 (MQ=255)                                       aGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAg                                                                                                           >  2:305866/1‑144 (MQ=255)                                       aGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAg                                                                                                           <  1:305866/144‑1 (MQ=255)                                          cGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTAc                                                                                                      >  5:184779/1‑146 (MQ=255)                                            tAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTGGGGTTTCATGGGGCGGGCCATTGTTGTTGATGACACCAGTGTCATGGTCCTGGGACCCGCCCAGTTTAGCACCAGTGTACAAgg                                                                                                  >  4:156927/1‑148 (MQ=255)                                                    ttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCGCCAGTGTACCAGGTGTTATCt                                                                                          >  6:136882/1‑148 (MQ=255)                                                    ttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCt                                                                                          >  2:84557/1‑148 (MQ=255)                                                    ttCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCt                                                                                          >  4:138078/1‑148 (MQ=255)                                                      cAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCttt                                                                                        >  5:312079/1‑148 (MQ=255)                                                       aaaGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCttt                                                                                        <  6:312079/147‑1 (MQ=255)                                                       aaaGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTc                                                                                       <  3:138078/148‑1 (MQ=255)                                                       aaaGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTc                                                                                       <  7:301750/148‑1 (MQ=255)                                                          gCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCgg                                                                                     <  7:45040/147‑1 (MQ=255)                                                           ccAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATAGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTcc                                                                                      >  8:45040/1‑144 (MQ=255)                                                                 tACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCt                                                                             >  1:188292/1‑148 (MQ=255)                                                                   cGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACTAGTGTCATGGTACTGGGGCCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTg                                                                            <  7:79239/147‑1 (MQ=255)                                                                             tGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCt                                                                        <  1:208076/141‑1 (MQ=255)                                                                             tGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCt                                                                        >  2:208076/1‑141 (MQ=255)                                                                                  accaccAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGaaa                                                            <  6:184779/148‑1 (MQ=255)                                                                                  accaccAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGaa                                                             <  5:65028/147‑1 (MQ=255)                                                                                      ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAATCCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccag                                                        <  8:299726/148‑1 (MQ=255)                                                                                      ccaAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccag                                                        <  5:287040/148‑1 (MQ=255)                                                                                        aaaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagcc                                                      <  1:84557/148‑1 (MQ=255)                                                                                        aaaaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagcc                                                      <  2:193219/148‑1 (MQ=255)                                                                                          aaGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagccag                                                    <  3:251423/148‑1 (MQ=255)                                                                                            gCACCAGCGCCCAGGGGGTTTTCATGGGTCGGGCCCTTGTTGGTGCTGAACCCAGTGTCATGGGACTGGGACCAGCCCAGTTTAGCACCAGTGTACCCGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTg                                                  <  3:156927/148‑1 (MQ=255)                                                                                              aCCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGcc                                                <  2:55823/148‑1 (MQ=255)                                                                                                        aGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGcc                                                >  7:267522/1‑138 (MQ=255)                                                                                                        aGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGcc                                                <  8:267522/138‑1 (MQ=255)                                                                                                         gTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagcc                                                      <  2:274018/131‑1 (MQ=255)                                                                                                         gTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAAccagcc                                                      >  1:274018/1‑131 (MQ=255)                                                                                                          ttGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATcgc                                     >  1:143260/1‑147 (MQ=255)                                                                                                           tGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGa                                   >  5:42259/1‑148 (MQ=255)                                                                                                                   aTGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGtt                                                                                              >  5:211893/1‑81 (MQ=255)                                                                                                                   aTGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGtt                                                                                              <  6:211893/81‑1 (MQ=255)                                                                                                                      ggTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCttt                        >  8:226311/1‑148 (MQ=255)                                                                                                                           ggCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCAt                   <  7:226311/148‑1 (MQ=255)                                                                                                                             ccATTGTTGTTGATGAAACCAGAGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTAACTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCCCTGCAAACGCGAAAGCAGTCTTTTTCAttt                 >  2:267012/1‑148 (MQ=255)                                                                                                                                   ttgttgATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTgcg           <  5:136882/148‑1 (MQ=255)                                                                                                                                      ttgATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCt        >  7:244401/1‑148 (MQ=255)                                                                                                                                      ttgATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCt        <  6:42259/148‑1 (MQ=255)                                                                                                                                      ttgATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCt        >  7:80006/1‑148 (MQ=255)                                                                                                                                      ttgATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCt        <  8:244401/148‑1 (MQ=255)                                                                                                                                         atgaAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCATCGt     >  8:203389/1‑148 (MQ=255)                                                                                                                                          tgaAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGtt    <  2:143260/148‑1 (MQ=255)                                                                                                                                          tgaAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGt     >  5:105543/1‑147 (MQ=255)                                                                                                                                            aaaCCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCg                                                               >  8:213048/1‑87 (MQ=255)                                                                                                                                            aaaCCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCg                                                               <  7:213048/87‑1 (MQ=255)                                                                                                                                            aaaCCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTat  >  7:117791/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                                GTTTTCAACGCTGCCTTTGTACGGCATACGACCTAACCAGTCGTAACCCATTTCAAAGCCAACATACGGGTTAACCTGGTAACCACCAAAAGCACCAGCGCCCAGTTGGTTTTCATGGGTCGGGCCATTGTTGTTGATGAAACCAGTGTCATGGTACTGGGACCAGCCCAGTTTAGCACCAGTGTACCAGGTGTTATCTTTCGGAGCGGCCTGCGCTACGGTAGCGAAACCAGCCAGTGCCACTGCAATCGCGATAGCTGTCTTTTTCATTTTTTGCGCCTCGTTAT  >  NC_000913/1019784‑1020070

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