Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,987,443	G→A	7.7%	S233N (AGT→AAT)	PP_1778 →	lipopolysaccharide ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,987,443	0	G	A	7.7%	188.1 / 4.1	65	S233N (AGT→AAT)	PP_1778	lipopolysaccharide ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (31/29);  new base A (2/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTACCGT  >  NC_002947/1987300‑1987587                                                                                                                                                |                                                                                                                                                   ggTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAgg                                                                                                                                               <  3:183745/149‑1 (MQ=255) ggTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAgg                                                                                                                                               <  1:306812/149‑1 (MQ=255)     ttgttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCgg                                                                                                                                            <  7:45946/148‑1 (MQ=255)       gttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑tt                                                                                                                                         <  4:305788/148‑1 (MQ=255)        ttCTGGTTTACCCCTATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCGC‑ttta                                                                                                                                       <  6:41278/149‑3 (MQ=255)        ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑t                                                                                                                                      >  2:40850/1‑149 (MQ=255)        ttCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGACTCAAATCAGGCGC‑ttta                                                                                                                                       <  6:220072/149‑3 (MQ=255)             gTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCAt                                                                                                                                   <  1:207374/147‑1 (MQ=255)                 aCCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCg                                                                                                                              <  1:147946/148‑1 (MQ=255)                  ccccAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCgg                                                                                                                             >  6:363509/1‑148 (MQ=255)                      aaTCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGAttt                                                                                                                        <  5:363509/149‑1 (MQ=255)                           tttACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCttt                                                                                                                   <  2:304479/149‑1 (MQ=255)                                 tGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCtg                                                                                                             <  1:40850/149‑1 (MQ=255)                                           gATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTGGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTg                                                                                                   >  7:100114/1‑149 (MQ=255)                                              aCCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcc                                                                                                >  8:196117/1‑149 (MQ=255)                                                cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCtgt                                                                                                            <  6:170536/135‑1 (MQ=255)                                                cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCtgt                                                                                                            >  5:170536/1‑135 (MQ=255)                                                cTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCAAGACCGGGATTTCCTTTGGACTGTTGTTGCTTGGactg                                                                                              >  1:323823/1‑149 (MQ=255)                                                 tGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgt                                                                                             >  5:263474/1‑149 (MQ=255)                                                     cacCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCTTGATTTGAGTCAGGCGG‑TTG‑Tc                                                                                                                                     <  2:342351/105‑1 (MQ=255)                                                      acCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctg                                                                                        >  8:112259/1‑149 (MQ=255)                                                        cTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGAACGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgttcctgtt                                                                                      >  4:27464/1‑149 (MQ=255)                                                                cACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGt                                                                              <  7:321910/149‑1 (MQ=255)                                                                   ccTTAGTTTTAACCCAATCTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTC‑TCATGAACGTGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTCTTTCGTCGGTCtg                                                                           >  4:334603/1‑149 (MQ=255)                                                                             aaCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAt                                                                 >  1:343974/1‑149 (MQ=255)                                                                              aCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtg                                                                <  1:355189/149‑1 (MQ=255)                                                                              aCCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTGGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAt                                                                 <  8:100114/148‑1 (MQ=255)                                                                                   aTGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGa                                                           >  4:303449/1‑149 (MQ=255)                                                                                      tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgtt                                                                                            >  3:217444/1‑113 (MQ=255)                                                                                      tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGcctgtt                                                                                            <  4:217444/113‑1 (MQ=255)                                                                                      tATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGACTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATg                                                         <  6:145241/148‑1 (MQ=255)                                                                                           cATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGGTTTGATCATGAACGGGATAACATTTTTTATTTTTTTGATTGGCATGTTAATGTTTAGTCTTTCTGTGCCTTAAATGTTGTATGTACt                                                     >  1:298225/1‑149 (MQ=255)                                                                                           cATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATTACCGGGATTTGCTTTGGTCTCTTGTTGATTGGCCTGTTCCTGTTTCGTCGGTCAGTGCATGAAATGGTGGATGTACTAt                                                   >  7:128261/1‑149 (MQ=255)                                                                                            aTTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATg                                                  >  8:219295/1‑149 (MQ=255)                                                                                                  aCTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtt                                            >  7:89623/1‑149 (MQ=255)                                                                                                    tGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                          >  2:124882/1‑149 (MQ=255)                                                                                                      ccTTTCAGGACGTACTCCTGTACGCTTATGCTCCTGATTTTAGTCAGGCGG‑TTG‑TCTTGAGGGGGTTTTGCTTTGGTCTTTTTTTGCTTGGCCTGTTCCTTTTTCGTCGGTCTGTGCGTGTAATGGTGGATGTACTATGAGCCtttt                                          <  3:334603/147‑1 (MQ=255)                                                                                                      ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                              >  1:245738/1‑127 (MQ=255)                                                                                                      ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAAtggt                                                              <  2:245738/127‑1 (MQ=255)                                                                                                      ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                          <  7:196117/147‑1 (MQ=255)                                                                                                      ccTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCtttt                                          >  6:149408/1‑147 (MQ=255)                                                                                                          tCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAgt                                      <  3:93859/147‑1 (MQ=255)                                                                                                            acgacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCa                                  <  3:315117/149‑1 (MQ=255)                                                                                                            acgacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGAGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTAGGCCTGTTCATGTTTCGTCGGTCTGTGCGTGAAATGGTGGATGTACTATGAGACTTTTGAGTGTCa                                  >  4:107084/1‑149 (MQ=255)                                                                                                              gacgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGt                                >  8:165777/1‑149 (MQ=255)                                                                                                               acgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATTACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGc                                               <  4:309156/133‑1 (MQ=255)                                                                                                               acgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGc                                               >  3:309156/1‑133 (MQ=255)                                                                                                                cgTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGa                              <  6:263474/149‑1 (MQ=255)                                                                                                                 gTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGAt                             >  8:147895/1‑149 (MQ=255)                                                                                                                     tCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTgg                         <  4:251182/149‑1 (MQ=255)                                                                                                                       cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTggg                        >  7:37416/1‑148 (MQ=255)                                                                                                                       cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTggg                        >  8:81792/1‑148 (MQ=255)                                                                                                                       cTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGc                       <  5:149408/149‑1 (MQ=255)                                                                                                                            cGCTAAAGCGCCTGTTTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTTCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGc                  <  2:298225/149‑1 (MQ=255)                                                                                                                               tAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCt                >  5:2642/1‑148 (MQ=255)                                                                                                                                 aaccgcCTGACTCAAATCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                        >  4:267664/4‑122 (MQ=255)                                                                                                                                 aaccgcCTGACTCAAATCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGa                                        <  3:267664/119‑1 (MQ=255)                                                                                                                                 aaGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCt                <  6:2642/146‑1 (MQ=255)                                                                                                                                 aaGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTAc              <  3:303449/148‑1 (MQ=255)                                                                                                                                  aGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACa             >  7:81446/1‑148 (MQ=255)                                                                                                                                   ccgcCTGACTCAAATCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGACTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACaa            >  4:341863/2‑148 (MQ=255)                                                                                                                                      ccTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACg        >  6:342857/1‑149 (MQ=255)                                                                                                                                         gATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGCCAAAGCCTACAAAACGTAc     >  1:306098/1‑149 (MQ=255)                                                                                                                                          aTTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTAc     <  7:112259/148‑1 (MQ=255)                                                                                                                                            ttGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTACCGt  >  4:309434/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                   GGTGTTGTTCTGGTTTACCCCAATCGTTTACATGGTCAATGTGATACCTGAACACCTTCGGGCCACCCTTAGTTTTAACCCCATGTATCCCATTGTTACTGCCTATCACGACGTACTCCTGTACGCTAAAGCGCCTGATTTGAGTCAGGCGG‑TTG‑TCATGACCGGGATTTCCTTTGGTCTGTTGTTGCTTGGCCTGTTCCTGTTTCGTCGGTCTGTGCCTGAAATGGTGGATGTACTATGAGCCTTTTGAGTGTCAGTGATGTGGGCAAAGCCTACAAAACGTACCGT  >  NC_002947/1987300‑1987587

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