Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,409,928	A→C	7.6%	intergenic (‑494/‑42)	PP_3022 ← / → PP_3024	AraC family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,409,928	0	A	C	7.6%	186.7 / 5.1	63	intergenic (‑494/‑42)	PP_3022/PP_3024	AraC family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (33/25);  new base C (3/2);  total (36/27)
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

GGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATGTGAA  >  NC_002947/3409783‑3410073                                                                                                                                                  |                                                                                                                                                   ggCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAAcc                                                                                                                                                 >  3:239325/1‑149 (MQ=255)   ccATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAAccc                                                                                                                                                >  3:315850/1‑148 (MQ=255)      tCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGc                                                                                                                                            <  1:263311/149‑1 (MQ=255)        aaGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGa                                                                                                                                          >  4:338285/1‑149 (MQ=255)                    gcgGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAcccag                                                                                                                                              >  5:98537/1‑131 (MQ=255)                    gcgGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAcccag                                                                                                                                              <  6:98537/133‑3 (MQ=255)                      ggCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAAc                                                                                                                            <  4:239325/149‑1 (MQ=255)                          tACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTa                                                                                                                        >  3:230928/1‑149 (MQ=255)                            ccTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAg                                                                                                                       <  7:114372/148‑1 (MQ=255)                                  cTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAg                                                                                                                <  3:338285/148‑1 (MQ=255)                                   tACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAggg                                                                                                              >  3:204486/1‑149 (MQ=255)                                    aCATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAggg                                                                                                              >  4:181512/1‑148 (MQ=255)                                     cATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCACCCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGt                                                                                                             >  7:217750/1‑148 (MQ=255)                                      aTCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTa                                                                                                            >  4:114417/1‑148 (MQ=255)                                       tCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAccca                                                                                                                                               >  1:308011/1‑112 (MQ=255)                                       tCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTTGCCGGCCGCCAccca                                                                                                                                               <  2:308011/113‑2 (MQ=255)                                          ccTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTg                                                                                                       <  2:225991/149‑1 (MQ=255)                                          ccTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTg                                                                                                       >  1:249728/1‑149 (MQ=255)                                                       gcgcgcACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTc                                                                                           <  5:135664/148‑1 (MQ=255)                                                               ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGc                                                                                      >  8:314388/1‑145 (MQ=255)                                                               ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGc                                                                                      <  7:314388/145‑1 (MQ=255)                                                               ccGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCAtt                                                                                   >  1:155054/1‑148 (MQ=255)                                                                 gTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGAc                                                                                <  2:249728/149‑1 (MQ=255)                                                                 gTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGAc                                                                                >  4:175796/1‑149 (MQ=255)                                                                  tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                               <  8:147669/149‑1 (MQ=255)                                                                  tGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACt                                                                               >  2:130749/1‑149 (MQ=255)                                                                    aGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTc                                                                              <  2:346915/148‑1 (MQ=255)                                                                      cGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTcagc                                                                           <  4:204486/149‑1 (MQ=255)                                                                        agaCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTcagcag                                                                         <  3:175796/149‑1 (MQ=255)                                                                          aCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAgg                                                                        <  4:230928/148‑1 (MQ=255)                                                                          aCCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGc                                                                       >  8:268422/1‑149 (MQ=255)                                                                             gtggtgGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                    >  6:322749/1‑149 (MQ=255)                                                                                 tggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                    <  6:7852/145‑1 (MQ=255)                                                                                 tggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGaa                                                                    >  5:7852/1‑145 (MQ=255)                                                                                  ggCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACagag                                                               >  2:22937/1‑149 (MQ=255)                                                                                   gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCCACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                              >  5:393064/1‑149 (MQ=255)                                                                                   gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                              >  1:208116/1‑149 (MQ=255)                                                                                   gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                              >  5:64767/1‑149 (MQ=255)                                                                                   gCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAgg                                                              >  2:275088/1‑149 (MQ=255)                                                                                     aaGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCCTCCATTGACTCAGCATGCGAACAGAGGc                                                             >  8:384243/1‑148 (MQ=255)                                                                                        cGATCACCCAGGCGTCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAAt                                                         <  6:64767/149‑1 (MQ=255)                                                                                              gccAGGCGCCAAGTGCCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCAc                                                   >  7:76303/2‑149 (MQ=255)                                                                                               ccAGGCGGCACTTGGCGCCTGGCGCGACAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACg                                                  >  2:118693/1‑149 (MQ=255)                                                                                                ctggcggcCCTTTGCGCCTGGCGCGCAAGGAGTTCTGGGTTGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGc                                                 <  4:315850/147‑1 (MQ=255)                                                                                                     ggcACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTc                                            >  7:104516/1‑149 (MQ=255)                                                                                                            ggCGCCTGGCTCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTcgtacgta                                     >  4:21279/1‑149 (MQ=255)                                                                                                            ggCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACt                                       <  3:114417/147‑1 (MQ=255)                                                                                                             gcgcCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTa                                     >  1:194223/1‑148 (MQ=255)                                                                                                                    gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCt                                                                                            <  8:6337/86‑1 (MQ=255)                                                                                                                    gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCt                                                                                            >  7:6337/1‑86 (MQ=255)                                                                                                                    gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATc                               <  2:155054/147‑1 (MQ=255)                                                                                                                    gcgcgcCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCa                              <  6:393064/148‑1 (MQ=255)                                                                                                                             gAGTTCTGGGTGGCGGCCGGCACCCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCCCGCTGGTCGTACTTTCAAATCACCTGCGCagt                    >  4:283578/1‑147 (MQ=255)                                                                                                                             gAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCAtt                    <  1:22937/149‑1 (MQ=255)                                                                                                                                   tGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCTTGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATgg              <  7:384243/149‑1 (MQ=255)                                                                                                                                   tGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATg               >  4:350340/1‑148 (MQ=255)                                                                                                                                     ggTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGCATTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCACATTGAAt                >  1:82905/1‑146 (MQ=255)                                                                                                                                     ggTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCACATTGAAt                <  2:82905/145‑1 (MQ=255)                                                                                                                                     ggTGGCGGCAGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGat            <  8:104516/149‑1 (MQ=255)                                                                                                                                       tGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGata           >  2:133317/1‑148 (MQ=255)                                                                                                                                           ggccggcAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAAAAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCAtg      >  6:123467/1‑149 (MQ=255)                                                                                                                                             ccggcAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCAtgtg    <  1:133317/149‑1 (MQ=255)                                                                                                                                               ggcAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGCCTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATGTGaa  >  5:290894/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                   GGCCATCAAGGTGGTGGGCGCGGCCTACCTGTTCTACATCGCCTGGCAACTGTGGCGCGCACCCGTGAGCGAGACCGTGGTGGCAAGCGATCACCCAGGCGGCACTTGGCGCCTGGCGCGCCAGGAGTTCTGGGTGGCGGCCGGCAACCCTAGCGACTGACAATATAAACTCTAGC‑TTCAGGGTAAATTGGTATCAACCCTCCATGCATTGACTCAGCAGGCGAACAGAGGCAAATTGCCACGCTTGTCGTACTTACAAATCACCTGCGCATTGAATGGATATCATGTGAA  >  NC_002947/3409783‑3410073

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