Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,403,484	G→A	6.0%	R33C (CGC→TGC)	PP_3882 ←	terminase small subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,403,484	0	G	A	6.0%	193.5 / 2.0	66	R33C (CGC→TGC)	PP_3882	terminase small subunit
Reads supporting (aligned to +/- strand):  ref base G (31/31);  new base A (2/2);  total (33/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.06e-01

GCAAGGATTTCGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGTTTTGCGCACGGAGGG  >  NC_002947/4403336‑4403631                                                                                                                                                     |                                                                                                                                                     gCAAGGATTTCGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATgcg                                                                                                                                                      >  7:181961/1‑149 (MQ=255)   aaGGATTTCGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCt                                                                                                                                                    >  8:291398/1‑149 (MQ=255)       aTTTCGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGGCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgc                                                                                                                                                <  7:135568/149‑1 (MQ=255)           cGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgctggc                                                                                                                                            >  6:47012/1‑149 (MQ=255)            ggAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgctggct                                                                                                                                           >  4:324579/1‑149 (MQ=255)             gAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgctggct                                                                                                                                           >  3:195367/1‑148 (MQ=255)              aaTAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGTGGCGCTGGGAATGCGCTTAGCTGGCTGa                                                                                                                                         <  1:20326/149‑1 (MQ=255)              aaTAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGa                                                                                                                                         <  4:40542/149‑1 (MQ=255)                  gCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGACTGAGGCt                                                                                                                                     >  5:140318/1‑149 (MQ=255)                   cAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgct                                                                                                                                               >  6:87968/1‑138 (MQ=255)                   cAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAgct                                                                                                                                               <  5:87968/138‑1 (MQ=255)                      tCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGtt                                                                                                                                 >  1:47492/1‑149 (MQ=255)                      tCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGtt                                                                                                                                 >  8:29498/1‑149 (MQ=255)                          acaGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑ttg                                                                                                                            <  8:181961/149‑1 (MQ=255)                           caGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑ttgt                                                                                                                           <  8:230491/149‑1 (MQ=255)                                 cTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCg                                                                                                                     <  6:56530/149‑1 (MQ=255)                                    aaGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCCGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCgg                                                                                                                    <  2:108773/147‑1 (MQ=255)                                      gcgcGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGgcgc                                                                                                                >  3:160484/1‑149 (MQ=255)                                        gcgAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCg                                                                                                                     >  6:275932/1‑142 (MQ=255)                                        gcgAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCg                                                                                                                     <  5:275932/142‑1 (MQ=255)                                              ccATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATc                                                                                                                      <  7:50898/135‑1 (MQ=255)                                              ccATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATc                                                                                                                      >  8:50898/1‑135 (MQ=255)                                              ccATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCgg                                                                                                         >  6:56065/1‑148 (MQ=255)                                               cATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGa                                                                                                       <  3:324579/149‑1 (MQ=255)                                                      cTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCctc                                                                                                <  4:87298/149‑1 (MQ=255)                                                         ggTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCctct                                                                                               <  7:291398/147‑1 (MQ=255)                                                         ggTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGa                                                                                             <  7:243047/149‑1 (MQ=255)                                                            aaGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTAACCGGGATTCCCTCTGAc                                                                                            >  1:343752/1‑147 (MQ=255)                                                              gCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGc                                                                                        <  4:195367/149‑1 (MQ=255)                                                              gCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGATGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGc                                                                                        >  4:183481/1‑149 (MQ=255)                                                               cacaCCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCAGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGc                                                                                        >  5:130821/1‑148 (MQ=255)                                                                acacCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCt                                                                                                                                                    >  7:79566/1‑88 (MQ=255)                                                                acacCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCt                                                                                                                                                    <  8:79566/88‑1 (MQ=255)                                                                  agccATGCTGTCAATCAGCACCGTGAGCCGCCCCCACGCCACATGGCCATCTTAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTTGCTGAGTCTCGTTC‑TTGTTGATCGGTCGCTTACCg                                                                                                          <  6:165688/125‑1 (MQ=255)                                                                  acCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCa                                                                                    <  4:160484/149‑1 (MQ=255)                                                                  acCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCCCTTACCg                                                                                                          >  5:165688/1‑127 (MQ=255)                                                                   cccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAg                                                                                   >  1:38621/1‑149 (MQ=255)                                                                   cccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCACTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCctct                                                                                               >  1:331738/1‑137 (MQ=255)                                                                   cccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCACTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCctct                                                                                               <  2:331738/137‑1 (MQ=255)                                                                    ccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCCCGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAAGGGGGGGAGCTGGCTGAGGCGCGGGC‑TGGTGGAGGGGGGGGTTACCGGGTTTCCGTGTGACAGGGTTCCGc                                                                                  >  8:95159/1‑149 (MQ=255)                                                                    ccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAg                                                                                   >  1:294948/1‑148 (MQ=255)                                                                    ccATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAg                                                                                   <  2:294948/148‑1 (MQ=255)                                                                          ggTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAg                                                                                   <  2:293042/142‑1 (MQ=255)                                                                          ggTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGATC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAg                                                                                   >  1:293042/1‑142 (MQ=255)                                                                           gTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCg                                                                           <  5:56065/149‑1 (MQ=255)                                                                                 cagcaCCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTgg                                                                     <  5:47012/149‑1 (MQ=255)                                                                                    caCCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCa                                                                                    >  1:161696/1‑131 (MQ=255)                                                                                    caCCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCa                                                                                    <  2:161696/131‑1 (MQ=255)                                                                                    caCCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGACTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGgtc                                                                  <  6:140318/149‑1 (MQ=255)                                                                                      ccGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTc                                                                                                                                    <  2:353624/82‑1 (MQ=255)                                                                                      ccGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTc                                                                                                                                    >  1:353624/1‑82 (MQ=255)                                                                                             ccgccCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGCTCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCgg                                                         >  7:95941/1‑149 (MQ=255)                                                                                                cccccACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCACTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCa                                                      >  8:60228/1‑149 (MQ=255)                                                                                                  cccACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCAtt                                                    >  3:137081/1‑149 (MQ=255)                                                                                                   ccACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCAttt                                                   >  2:200985/1‑149 (MQ=255)                                                                                                   ccACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGAGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGTGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCAttt                                                   >  2:195395/1‑149 (MQ=255)                                                                                                    cACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCAtttt                                                  <  6:162961/149‑1 (MQ=255)                                                                                                               gCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCaaa                                       <  7:298348/149‑1 (MQ=255)                                                                                                                  aTCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCaaaaa                                     >  6:71529/1‑148 (MQ=255)                                                                                                                      gAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTCTTTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCaaaaaaaa                                  >  1:279112/1‑148 (MQ=255)                                                                                                                            aaGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAAtcatttc                           <  2:47492/148‑1 (MQ=255)                                                                                                                            aaGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAAtcatttc                           <  7:29498/148‑1 (MQ=255)                                                                                                                             aGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAAtcatttcat                         >  6:356474/1‑149 (MQ=255)                                                                                                                                   cGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTcgcg                   >  1:81363/1‑149 (MQ=255)                                                                                                                                    gggggCGCAGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGTGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTcgcg                   <  6:130821/148‑1 (MQ=255)                                                                                                                                        gcgcTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGtttt              <  2:38621/149‑1 (MQ=255)                                                                                                                                           cTGGGAATGCACTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGTTTTgcg           <  7:60228/149‑1 (MQ=255)                                                                                                                                                    cgcTTAGCTGGCTGAGGCGCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGTTTTGCGCACGGAggg  <  1:200985/149‑1 (MQ=255)                                                                                                                                                     gcTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGTTTTGCGCACGGAggg  <  1:195395/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCAAGGATTTCGGAATAGCAGTCGCACAGGCGCTCAAGCGCGAAGCCATCTGCCTCGGTAAGCACACCCATGCGGTCAAGCAGCACCGTGAGCCGCCCCCACGCCACCTGGCCATCTGAGCTCAAGTGGTCGGGGGCGCTGGGAATGCGCTTAGCTGGCTGAGGCTCGTTC‑TTGTTGATCGGGCGCTTACCGGGATTCCCTCTGACAAGCTTCAGCTCCGTCGGTGTGGGTCGTCTCCCGGCCATTTTTTGCTCTCAAAAAAAATCATTTCATTTCGCGGTTTTGCGCACGGAGGG  >  NC_002947/4403336‑4403631

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