Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	639,838	T→A	6.4%	S303C (AGC→TGC)	PP_0551 ←	MORN domain‑containing protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	639,838	0	T	A	6.4%	185.6 / 2.1	62	S303C (AGC→TGC)	PP_0551	MORN domain‑containing protein
Reads supporting (aligned to +/- strand):  ref base T (28/30);  new base A (2/2);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

CTTTCGACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACG  >  NC_002947/639690‑639986                                                                                                                                                        |                                                                                                                                                     cTTTCGACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCt                                                                                                                                                      <  3:169396/149‑1 (MQ=255) cTTTCGACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCt                                                                                                                                                      <  3:43745/149‑1 (MQ=255)   ttCGACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGc                                                                                                                                                    >  8:90261/1‑149 (MQ=255)       aCCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATc                                                                                                                                                <  8:241592/149‑1 (MQ=255)       aCCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCAGCGAtg                                                                                                                                                >  1:259604/1‑148 (MQ=255)          ttGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGATGGGACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGc                                                                                                                                              >  1:111907/1‑149 (MQ=255)                 tGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAg                                                                                                                                       >  5:1984/1‑148 (MQ=255)                   ggCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCAGCGATGGCAGccactaca                                                                                                                                    <  2:259604/149‑7 (MQ=255)                    gCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAAc                                                                                                                                   >  3:133017/1‑149 (MQ=255)                      aGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTc                                                                                                                                 >  6:64018/1‑149 (MQ=255)                       gAATCAGCCGGCCATGGCCATGGTTGGCATCGTTGAGGAAACCGCCCTT‑GTACTTGGTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACATTCCGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTc                                                                                                                                 <  1:97565/149‑1 (MQ=255)                          tCAGCCGCCCATGGTCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCAGCGATGGCAGcca                                                                                                                                         >  4:329657/1‑136 (MQ=255)                          tCAGCCGCCCATGGTCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCAGCGATGGCAGcca                                                                                                                                         <  3:329657/137‑2 (MQ=255)                            aGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTc                                                                                                                           >  4:45511/1‑149 (MQ=255)                               cgccCATGGCCATTGTAGGCCTCTGTGAGGAAACAGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTTCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGcc                                                                                                                        <  7:131596/149‑1 (MQ=255)                                       tccatggTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCGCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                 <  3:309565/147‑1 (MQ=255)                                        ccatggTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGa                                                                                                               <  7:91795/149‑1 (MQ=255)                                           tggTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGcc                                                                                                            >  2:192247/1‑149 (MQ=255)                                                  cATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGAAGTATGTTACCTTGAAGTGGGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                      >  1:121645/1‑149 (MQ=255)                                                   aTCGTTGAGTAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGa                                                                                                    <  2:136908/149‑1 (MQ=255)                                                   aTCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTg                                                                                                     >  2:360726/1‑148 (MQ=255)                                                      gTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAAtt                                                                                                 >  4:246383/1‑149 (MQ=255)                                                          aGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCAc                                                                                              <  3:45511/148‑1 (MQ=255)                                                            gAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGc                                                                                                                         <  8:197028/119‑1 (MQ=255)                                                            gAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGc                                                                                                                         >  7:197028/1‑119 (MQ=255)                                                               aCCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                 <  2:86153/124‑1 (MQ=255)                                                               aCCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGc                                                                                                                 >  1:86153/1‑124 (MQ=255)                                                                 cgccgAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                      <  2:70303/149‑1 (MQ=255)                                                                 cgccgAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                      <  1:251674/149‑1 (MQ=255)                                                                 cgccgAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                      <  5:348950/149‑1 (MQ=255)                                                                 cgccgAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCcac                                                                                      <  2:111907/149‑1 (MQ=255)                                                                   ccgAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCacac                                                                                    <  2:116792/149‑1 (MQ=255)                                                                           aCTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCa                                                                             >  6:361347/1‑149 (MQ=255)                                                                             ttGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATc                                                                           >  4:46768/1‑149 (MQ=255)                                                                               gctgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCgg                                                                         >  1:81079/1‑149 (MQ=255)                                                                                ctgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑CACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCgg                                                                         >  5:306466/1‑148 (MQ=255)                                                                                 tgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                      >  8:192067/1‑118 (MQ=255)                                                                                 tgcCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCtt                                                                                                      <  7:192067/118‑1 (MQ=255)                                                                                        ggCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGt                                                                <  1:192247/149‑1 (MQ=255)                                                                                           cAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTc                                                                                 <  8:207871/129‑1 (MQ=255)                                                                                           cAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTc                                                                                 >  7:207871/1‑129 (MQ=255)                                                                                           cAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGc                                                             >  2:209431/1‑149 (MQ=255)                                                                                              ttGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggc                                                          <  3:188102/149‑1 (MQ=255)                                                                                              ttGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggc                                                          >  6:205648/1‑149 (MQ=255)                                                                                                        ccTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGcggcctcggcct                                                  <  3:46768/147‑1 (MQ=255)                                                                                                        ccTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                <  2:121645/149‑1 (MQ=255)                                                                                                         cTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                >  7:336574/1‑148 (MQ=255)                                                                                                          ttGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTc                                                >  4:262007/1‑147 (MQ=255)                                                                                                               ccGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                         <  3:246383/149‑1 (MQ=255)                                                                                                               ccGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                         <  1:269461/149‑1 (MQ=255)                                                                                                               ccGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCgg                                         <  6:1984/149‑1 (MQ=255)                                                                                                                    a‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTc                                    >  1:368538/1‑149 (MQ=255)                                                                                                                    a‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTc                                    <  5:361347/149‑1 (MQ=255)                                                                                                                        cGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCtt                                 <  4:224659/149‑1 (MQ=255)                                                                                                                        cGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCtt                                 <  7:90261/149‑1 (MQ=255)                                                                                                                                ggCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCCGCCTTCCAGGCGGTTGTCCTTGAATTCg                          <  1:209431/148‑1 (MQ=255)                                                                                                                                 gCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGcc                        >  6:354403/1‑149 (MQ=255)                                                                                                                                   gAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGa                      <  5:64018/149‑1 (MQ=255)                                                                                                                                           ccTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGc              >  3:95065/1‑149 (MQ=255)                                                                                                                                             ttGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCAGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCTTTTTCGTAGCGTCCTCGGCATTCCAGGTGGTTGTCCTTGAATTCGCCGATGTACTGCTc            >  5:272959/1‑149 (MQ=255)                                                                                                                                               gTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTc            <  2:185143/147‑1 (MQ=255)                                                                                                                                               gTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCg           <  3:262007/148‑1 (MQ=255)                                                                                                                                               gTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGc          >  4:259396/1‑149 (MQ=255)                                                                                                                                                  gT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCAt       >  4:5767/1‑149 (MQ=255)                                                                                                                                                        tGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACg  >  4:353873/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                     CTTTCGACCTTGCCACTGGGCAGAATCAGCCGCCCATGGCCATGGTAGGCATCGTTGAGGAAACCGCCGAT‑GTACTTGCTGCCGTCGGCCAGTTGCAGGCTGCCTTGGCCCGACA‑GACGCCAGTCGGCGAACTCGCCCTTGTAGT‑GGCTGCCATCGCTGCCCAGCAACTCCCCTTCGCCGATCAGCGAGCCGTCCTTGAATTCACCGATCCACACGTCGCCATCGGCATTTTCGTAGCGGCCTCGGCCTTCCAGGCGGTTGTCCTTGAATTCGCCGATGTACTGCTCGCCATCAACG  >  NC_002947/639690‑639986

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