Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,447,423	T→C	5.9%	A1104A (GCT→GCC)	tamB →	translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,423	0	T	C	5.9%	208.0 / 2.5	68	A1104A (GCT→GCC)	tamB	translocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (36/28);  new base C (2/2);  total (38/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTCAGGACAAACAAGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCTACGT  >  NC_000913/4447277‑4447558                                                                                                                                                   |                                                                                                                                        gTTCAGTACAATCATTGTCTGGGTCTGATCGGGCAGATCAACATCCGTGAAGTGCTATTCAAAGCCTATGGTCAGGATCTGATTGTGCGTAATGGTCAGTTACTGATCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAATCt                                                                                                                                         <  4:5636/147‑1 (MQ=255)   tCAGGACAAACAAGGGCTGGGGCGGCACGGGCAGATCACCAGCCCTGAAGGGCGCTTCCATGCCTATCGTGAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTAtt                                                                                                                                      >  6:366715/1‑148 (MQ=255)    cAGGACAAACAAGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTAtt                                                                                                                                      <  4:375410/147‑1 (MQ=255)           aacaaGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAAccc                                                                                                                              >  5:150230/1‑148 (MQ=255)            acaaGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTAtt                                                                                                                                      <  1:261525/139‑1 (MQ=255)            acaaGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTAtt                                                                                                                                      >  2:261525/1‑139 (MQ=255)             caaGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCgg                                                                                                                            <  1:181134/148‑1 (MQ=255)                     gggTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGGCAGGATCTGATTGTGCCTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTa                                                                                                                      >  5:99054/1‑146 (MQ=255)                      ggTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACa                                                                                                                    <  5:271609/147‑1 (MQ=255)                                      tCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCCATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGc                                                                                                   >  5:142154/1‑148 (MQ=255)                                      tCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCAGCCAGATCAACCGTATCTTAATATTGAAGCCATTCGTAAACAGGATGCTACAGAAGACGACGTAATctc                                                                                                   >  5:136343/1‑146 (MQ=255)                                      acAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTTCGTAAAGGTGAGTTACTGTTCTATGGTCCGCCAGATCAACCGTATCTAAATATATAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGc                                                                                                   >  3:220261/2‑148 (MQ=255)                                               cTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTcg                                                                                          >  8:185016/1‑148 (MQ=255)                                                 gAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTcgcg                                                                                        >  8:312334/1‑148 (MQ=255)                                                      gcgcTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACt                                                                                   >  1:180047/1‑148 (MQ=255)                                                      gcgcTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACt                                                                                   <  6:150230/148‑1 (MQ=255)                                                      gcgcTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACAGTATATTAATATTGAAGCTATTCGTAACCCGGATGCTACaga                                                                                                                  <  5:130198/117‑1 (MQ=255)                                                       cgcTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGATACagaa                                                                                                                 >  6:130198/1‑117 (MQ=255)                                                         cTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGt                                                                                <  1:303829/148‑1 (MQ=255)                                                            ccATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTACTCGCCGGAGTTCGCGTCACTGGTCTg                                                                             >  6:265941/1‑148 (MQ=255)                                                             cATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCATCCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCTCGTCACTGGTCTgg                                                                            >  3:233609/1‑148 (MQ=255)                                                               tGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCg                                                                          >  7:330168/1‑148 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCgg                                                                                                                            >  2:184512/1‑97 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCgg                                                                                                                            <  1:184512/97‑1 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTgg                                                                                 >  1:286822/1‑140 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTgg                                                                                 <  2:286822/140‑1 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCg                                                                          >  7:142468/1‑147 (MQ=255)                                                                gCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCg                                                                          >  7:362052/1‑147 (MQ=255)                                                                    aTGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACg                                                                      >  3:9216/1‑147 (MQ=255)                                                                     tGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCCATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGaa                                                                    <  6:142154/148‑1 (MQ=255)                                                                     tGGTAATGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCCATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGaa                                                                    <  6:136343/148‑1 (MQ=255)                                                                                  aTTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGgaga                                                       <  2:180047/148‑1 (MQ=255)                                                                                  aTTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGgag                                                        <  5:265941/147‑1 (MQ=255)                                                                                   ttGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAt                                                      <  8:330168/148‑1 (MQ=255)                                                                                    tgtgCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAt                                                      <  4:54085/147‑1 (MQ=255)                                                                                    tgtgCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAt                                                      >  3:54085/1‑147 (MQ=255)                                                                                      tgCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtctt                                                   <  7:185016/148‑1 (MQ=255)                                                                                       gCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  >  6:127022/1‑148 (MQ=255)                                                                                       ccGTAAAGGGGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGACGCTATTCGTAACCCGGATGCTACAGAAGCCGACGTAATCGCCGGAGTACGCGTCACTGGGCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  <  8:413210/147‑1 (MQ=255)                                                                                           aaaGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTg                                              >  8:95032/1‑148 (MQ=255)                                                                                            aaGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGa                                             >  8:22745/1‑148 (MQ=255)                                                                                                tGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCg                                         >  4:29710/1‑148 (MQ=255)                                                                                                tGAGTTACTCTTCTCTGGTCCGACAGATCAACCGTATCTTAATATTGAATATATTCGTAACCAGTATGGTTCAGAAGACGACGTAATAGCCTGAGTTGGCGTCACTTGTCTGGCTGACGAAGCGAAATCGGATATCTTCTCTGAccc                                          <  4:220261/147‑1 (MQ=255)                                                                                                 gAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCgg                                        >  4:183861/1‑148 (MQ=255)                                                                                                    ttACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGa                                     <  8:362052/148‑1 (MQ=255)                                                                                                     tACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGCCCCGGCGAt                                    >  6:26064/1‑148 (MQ=255)                                                                                                       cTGGTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGAGGACGTATTCGCCTGAGTTCGCGTCAGTGGGCTGGCGGACGAACCGAAAGCGGTGATCTTCTCTGACCCGGCGATGt                                  <  6:99054/148‑1 (MQ=255)                                                                                                         gTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCg                                >  8:400482/1‑148 (MQ=255)                                                                                                            ctctGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGca                              <  5:127022/147‑1 (MQ=255)                                                                                                                 gTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  <  6:74218/122‑1 (MQ=255)                                                                                                                 gTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  >  5:74218/1‑122 (MQ=255)                                                                                                                    cGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTg                      <  7:22745/147‑1 (MQ=255)                                                                                                                     gCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGcc                    <  4:9216/148‑1 (MQ=255)                                                                                                                      ccAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCt                   >  6:45054/1‑148 (MQ=255)                                                                                                                          aTCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTc               >  7:274355/1‑148 (MQ=255)                                                                                                                           tCAACCGTATCTTAATATTGAAGCTATTTGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  >  1:170880/1‑112 (MQ=255)                                                                                                                           tCAACCGTATCTTAATATTGAAGCTATTTGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGAtcttc                                                  <  2:170880/112‑1 (MQ=255)                                                                                                                           tCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCt              >  2:201744/1‑148 (MQ=255)                                                                                                                              aCCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGc                     >  4:238451/1‑138 (MQ=255)                                                                                                                              aCCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGc                     <  3:238451/138‑1 (MQ=255)                                                                                                                                cGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTAttt         >  1:349242/1‑148 (MQ=255)                                                                                                                                cGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTAttt         >  3:165862/1‑148 (MQ=255)                                                                                                                                 gTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTg        >  4:135067/1‑148 (MQ=255)                                                                                                                                    tCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCt      <  5:45054/147‑1 (MQ=255)                                                                                                                                      ttAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATg                                   <  6:79577/116‑1 (MQ=255)                                                                                                                                      ttAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATg                                   >  5:79577/1‑116 (MQ=255)                                                                                                                                      ttAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCTACg   <  2:349242/148‑1 (MQ=255)                                                                                                                                       tAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCTACGt  >  7:258947/1‑148 (MQ=255)                                                                                                                                        aaTATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCTACGt  <  8:258947/147‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                        GTTCAGGACAAACAAGGGCTGGGTCTGAACGGGCAGATCAACATCCCTGAAGGGCGCTTCCATGCCTATGGTCAGGATCTGATTGTGCGTAAAGGTGAGTTACTGTTCTCTGGTCCGCCAGATCAACCGTATCTTAATATTGAAGCTATTCGTAACCCGGATGCTACAGAAGACGACGTAATCGCCGGAGTTCGCGTCACTGGTCTGGCGGACGAACCGAAAGCGGAGATCTTCTCTGACCCGGCGATGTCGCAACAAGCTGCCTTGTCTTATTTGCTACGT  >  NC_000913/4447277‑4447558

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: ‑