Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,599,290	G→C	11.4%	L173L (CTC→CTG)	PP_4071 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,599,290	0	G	C	11.4%	204.6 / 14.2	70	L173L (CTC→CTG)	PP_4071	membrane protein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base C (2/6);  total (34/36)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCG  >  NC_002947/4599147‑4599429                                                                                                                                                |                                                                                                                                            gCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAgcgc                                                                                                                                        <  2:353416/149‑1 (MQ=255)    tGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCCCCACCCCACGGGGTAGGCCAgcgcagg                                                                                                                                     >  8:218879/1‑146 (MQ=255)        ttCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCt                                                                                                                                 <  7:20943/149‑1 (MQ=255)          ccAGCCGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATAACCGATTAATTGCTTATCCACCACCCCAAGGGGTAGGCGCTCGCCTGCGCTgg                                                                                                                               >  8:155008/1‑149 (MQ=255)            aGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTggcc                                                                                                                             >  4:131362/1‑149 (MQ=255)            aGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCCAGCGCAGGCGCTcgcc                                                                                                                             <  2:354150/149‑1 (MQ=255)                ggTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCCAGCGCAGGCGCTcgcctacc                                                                                                                         <  7:151462/149‑5 (MQ=39)                 gTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGccc                                                                                                                        <  8:349527/149‑1 (MQ=255)                 gTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGccc                                                                                                                        <  3:181563/149‑1 (MQ=255)                      gAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCCAGCGCAGGCGCTCGCCTACCCCGTgg                                                                                                                   <  7:218879/149‑2 (MQ=38)                       aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGTGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGcc                                                                                                                  >  4:133088/1‑149 (MQ=255)                       aaGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTggg                                                                                                                  <  5:352507/149‑3 (MQ=38)                        aGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGCCGGCCAGCGCAGGCGCTCGCCTACCCCGTgggg                                                                                                                 <  4:77359/149‑4 (MQ=38)                            tCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             >  4:209013/1‑149 (MQ=255)                                 gTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGt                                                                                                        <  6:362429/149‑1 (MQ=255)                                 gTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGt                                                                                                        >  5:303184/1‑149 (MQ=255)                                  tCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTc                                                                                                       <  7:291561/149‑1 (MQ=255)                                        gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGc                                                                                                  <  8:275710/148‑1 (MQ=255)                                        gAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCg                                                                                                 >  4:141432/1‑149 (MQ=255)                                         aGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATGATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCgg                                                                                                <  2:174596/149‑1 (MQ=255)                                            cACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCCAgcgcag                                                                                                                                      <  3:185080/108‑3 (MQ=255)                                            cACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCCAgcgcag                                                                                                                                      >  4:185080/1‑106 (MQ=255)                                             aCTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                            <  5:228173/133‑1 (MQ=255)                                             aCTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGcc                                                                                                            >  6:228173/1‑133 (MQ=255)                                                 cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                          <  3:141432/147‑1 (MQ=255)                                                 cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAAc                                                                                          <  4:171953/147‑1 (MQ=255)                                                 cccAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCATTGTCGATGCGGGTCAAc                                                                                          <  4:34489/147‑1 (MQ=255)                                                   cAGCGAATCTGCGGGGCCGTTCTGGCTGATAAACGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCCCCTGCGCTGCCCGGCCCCGTGCCTGCGCCTTGTCGATGGGGGTCAACgag                                                                                       >  8:117106/1‑148 (MQ=255)                                                    aGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACgaga                                                                                      >  1:163712/1‑148 (MQ=255)                                                     gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             >  4:271227/1‑124 (MQ=255)                                                     gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTgcgc                                                                                                             <  3:271227/124‑1 (MQ=255)                                                     gCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGaaa                                                                                    <  7:328141/149‑1 (MQ=255)                                                        aaTCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACtt                                                                                 >  7:11059/1‑149 (MQ=255)                                                              cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATc                                                                            <  3:209013/148‑1 (MQ=255)                                                              cGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCg                                                                           <  5:300512/149‑1 (MQ=255)                                                                      ttCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGt                                                                    >  3:358683/1‑148 (MQ=255)                                                                           gCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATa                                                               >  3:316276/1‑148 (MQ=255)                                                                           gCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATa                                                               <  2:163712/148‑1 (MQ=255)                                                                              gATACATGCACCGGGCTGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGTCTAGCGCCTGCGCTGGCCGGCCCAGTTCCTGCACCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAg                                                            <  7:155008/148‑1 (MQ=255)                                                                                  aaGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAgctgc                                                        >  5:59704/1‑148 (MQ=255)                                                                                   aGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGctc                                                      >  4:82521/1‑149 (MQ=255)                                                                                     gCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCtt                                                    >  5:251780/1‑149 (MQ=255)                                                                                       aCCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGc                                                  >  6:118638/1‑149 (MQ=255)                                                                                         cgggcgggTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCg                                                >  3:334971/1‑149 (MQ=255)                                                                                            gcgggTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcc                                             >  6:73029/1‑149 (MQ=255)                                                                                                gTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgcc                                         >  1:293457/1‑149 (MQ=255)                                                                                                 tCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgcct                                        <  4:316276/149‑1 (MQ=255)                                                                                                   gatgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                       <  5:73029/148‑1 (MQ=255)                                                                                                   gatgatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGcctgcctg                                       <  3:82521/148‑1 (MQ=255)                                                                                                      gatCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                   >  4:129984/1‑149 (MQ=255)                                                                                                       atCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGcc                                  >  1:270838/1‑149 (MQ=255)                                                                                                        tCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                   >  6:167884/1‑147 (MQ=255)                                                                                                        tCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGgcgc                                   <  5:167884/147‑1 (MQ=255)                                                                                                         cACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgcct                                        >  1:5832/1‑141 (MQ=255)                                                                                                         cACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGgcctgcct                                        <  2:5832/141‑1 (MQ=255)                                                                                                            cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGt                                                                                                        <  6:94318/74‑1 (MQ=255)                                                                                                            cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGt                                                                                                        >  5:94318/1‑74 (MQ=255)                                                                                                            cGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCa                             <  4:358683/149‑1 (MQ=255)                                                                                                            cGCATAATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGCCCGGCCCCGTGCCTGCGCCTTGTCCATGCGGGTCAACGAGAAACTTTGATCGGTGCTGt                                                                    >  8:4334/1‑110 (MQ=255)                                                                                                               atcatTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGa                          <  8:11059/149‑1 (MQ=255)                                                                                                                  atTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGc                                                  >  8:127827/1‑122 (MQ=255)                                                                                                                  atTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGc                                                  <  7:127827/122‑1 (MQ=255)                                                                                                                   ttGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTa                      <  5:118638/149‑1 (MQ=255)                                                                                                                       tGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCAc                  >  3:193833/1‑149 (MQ=255)                                                                                                                       tGATCCACCAACCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCAc                  >  1:298267/1‑149 (MQ=255)                                                                                                                         aTCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACgg                <  4:334971/149‑1 (MQ=255)                                                                                                                             accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGCAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAAt            >  8:292366/1‑149 (MQ=255)                                                                                                                             accaccCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAAt            <  3:131362/149‑1 (MQ=255)                                                                                                                                    cACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTa       >  3:214398/1‑147 (MQ=255)                                                                                                                                       ggggTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGACCAGCACGAGGTACCACGGCAATGGGTACAgag  >  8:154618/1‑147 (MQ=255)                                                                                                                                                |                                                                                                                                            GCTTGGGTTCCAGACGGTACCGAAGACTCAACGTCAGTGGAGTCACTGCCCAGCGAATCTGCGGGGCCGTTCTGGCTGATAAAGGCACCGGGCGGGTCGATGATCACCGCATCATTGCTGATCCACCACCCCACGGGGTAGGCGAGCGCCTGCGCTGGCCGGCCCCGTGCCTGCGCCTTGTCGATGCGGGTCAACGAGAAACTTTGATCGGTGCAGTCAATAAAGCTGCTCTTGCCGGCCTGCCTGGCGCCCAGCACGAGGTACCACGGCAATCGGTACAGCG  >  NC_002947/4599147‑4599429

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