Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,518,553	A→T	10.9%	C10S (TGC→AGC)	PP_4850 ←	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,518,553	0	A	T	10.9%	176.0 / 11.1	64	C10S (TGC→AGC)	PP_4850	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (25/32);  new base T (4/3);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 6.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.89e-01

AACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGG  >  NC_002947/5518405‑5518698                                                                                                                                                      |                                                                                                                                                  aaCACCAGCAGACCGGCCTGCAGGCTTGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCa                                                                                                                                                   <  3:24458/149‑1 (MQ=255)   cacCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGt                                                                                                                                                 <  1:325347/149‑1 (MQ=255)    acCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTc                                                                                                                                                >  5:101953/1‑149 (MQ=255)      cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTg                                                                                                                                              <  6:101953/149‑1 (MQ=255)      cagcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGCCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGGAGTATACAGCCTCAACGCCACACAGCGCGAATGCCGCaaaat                                                                                                                                               >  1:119794/1‑145 (MQ=255)       agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGc                                                                                                                                             <  5:72283/149‑1 (MQ=255)       agcagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGc                                                                                                                                             >  8:28448/1‑149 (MQ=255)         cagACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCa                                                                                                                                           <  7:28448/149‑1 (MQ=255)                tccAGCAGGGTGGTGATCAACTTCAGGGTCAAAGCGAGCAGCAGGGCGGCGATGACGATGATCAGGTTGGCGAGGGTGATCTGTCTGTGGGTAAAGG‑TCATTCAGCCTCGACGGCAGACAGCTGTAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                                    <  5:361345/148‑1 (MQ=255)                  cagcagGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTggcg                                                                                                                                  <  1:45229/149‑1 (MQ=255)                          tGGGCAGCAACTTCAGGTGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                                    <  5:182186/139‑1 (MQ=255)                          tGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGCCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTgg                                                                                                                                    >  6:182186/1‑139 (MQ=255)                           gggCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCa                                                                                                                         >  4:297250/1‑149 (MQ=255)                              cagcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGAGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAg                                                                                                                      >  4:293446/1‑149 (MQ=255)                               agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                                     >  4:204075/1‑149 (MQ=255)                               agcaACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAgg                                                                                                                     <  4:82951/149‑1 (MQ=255)                                  aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                         <  2:272953/126‑4 (MQ=38)                                  aaCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACTGCGGCATTCCCGCTGTCTGccgtc                                                                                                                                         >  1:272953/1‑123 (MQ=38)                                   aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                            <  5:296231/138‑1 (MQ=255)                                   aCTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAg                                                                                                                            >  6:296231/1‑138 (MQ=255)                                       cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                                 <  5:93106/145‑1 (MQ=255)                                       cAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCa                                                                                                                 >  6:93106/1‑145 (MQ=255)                                           ggCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCAttt                                                                                                         >  5:285884/1‑149 (MQ=255)                                            gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                                        >  1:149683/1‑149 (MQ=255)                                            gCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttc                                                                                                        <  1:183823/149‑1 (MQ=255)                                               aaGCCAGCAGCAGCCCGGCCATTACGATGATCAGGCTGGCGAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGAGGGCAGAGAGCGGTAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                                     <  1:97693/149‑1 (MQ=255)                                               aaGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttc                                                                                                     <  3:11807/149‑1 (MQ=255)                                                 gCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGTCGGCAGACAGCGGGAAGGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttcttctt                                                                                                   <  7:157570/149‑1 (MQ=255)                                                   cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                                 <  2:51582/149‑1 (MQ=255)                                                   cagcagcagCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                                 <  1:8709/149‑1 (MQ=255)                                                     gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGTAGGTCATTTCTTCTTCAgg                                                                                               <  6:285884/149‑1 (MQ=255)                                                     gcagcagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAgg                                                                                               <  3:297250/149‑1 (MQ=255)                                                         cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                                      <  2:34218/138‑1 (MQ=255)                                                         cagcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATttctt                                                                                                      >  1:34218/1‑138 (MQ=255)                                                           gcCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCg                                                                                         >  7:146843/1‑149 (MQ=255)                                                            cccGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGc                                                                                        >  8:70816/1‑149 (MQ=255)                                                               ggCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCa                                                                                     <  2:149683/149‑1 (MQ=255)                                                                gCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCAt                                                                                    >  7:92053/1‑149 (MQ=255)                                                                          tgatCTGGCTGGCCAGGGTGTTCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCTGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTg                                                                          <  2:119794/149‑1 (MQ=255)                                                                           gatCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                         >  3:76914/1‑149 (MQ=255)                                                                            atCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgc                                                                         <  8:92053/148‑1 (MQ=255)                                                                             tCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTgcgc                                                                       <  8:146843/149‑1 (MQ=255)                                                                                        aGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCtt                                                            >  1:173368/1‑149 (MQ=255)                                                                                           gTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAg                                                         >  8:118308/1‑149 (MQ=255)                                                                                             gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                              >  8:210884/1‑126 (MQ=255)                                                                                             gATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGc                                                                              <  7:210884/126‑1 (MQ=255)                                                                                                cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                                 >  8:274989/1‑104 (MQ=255)                                                                                                cTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCa                                                                                                 <  7:274989/104‑1 (MQ=255)                                                                                                 tGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCttt                                                   >  3:20985/1‑149 (MQ=255)                                                                                                     ggggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtg                                               <  7:118308/149‑1 (MQ=255)                                                                                                       ggggTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggt                                             <  4:76914/149‑1 (MQ=255)                                                                                                          gTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGc                                          <  2:173368/149‑1 (MQ=255)                                                                                                               gg‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAg                                     >  8:88393/1‑149 (MQ=255)                                                                                                                       cAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGc                               <  4:20985/148‑1 (MQ=255)                                                                                                                        aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                           <  3:102302/135‑1 (MQ=255)                                                                                                                        aGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtggtggtgg                                           >  4:102302/1‑135 (MQ=255)                                                                                                                               ctggCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCa                      >  8:225392/3‑149 (MQ=255)                                                                                                                                  gCAGACTGCGGCATTCCCGCTGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGGGGTGGCGTCAGCATTGGAGGTATTCaata                   >  3:26475/1‑146 (MQ=255)                                                                                                                                    aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                                 >  2:126337/1‑117 (MQ=255)                                                                                                                                    aGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTtgg                                                 <  1:126337/117‑1 (MQ=255)                                                                                                                                     gACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCAtttttt                >  6:95453/1‑149 (MQ=255)                                                                                                                                        aGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCtt             <  3:204075/149‑1 (MQ=255)                                                                                                                                          cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc           >  7:16252/1‑149 (MQ=255)                                                                                                                                          cGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGc           >  8:288163/1‑149 (MQ=255)                                                                                                                                                  ccGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGttt    <  5:95453/148‑1 (MQ=255)                                                                                                                                                   cGCAGTCTGCCAGCCTTGGCGGCGCAGGCATACGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTgg  <  3:293446/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  AACACCAGCAGACCGGCCAGCAGGCTGGGCAGCAACTTCAGGGGCAAAGCCAGCAGCAGCCCGGCCATGACGATGATCAGGCTGGCCAGGGTGATCTGACGGGGGGTAAAGG‑TCATACAGCCTCGACGGCAGACAGCGGGAATGCCGCAGTCTGCCAGCCTTGGCGGCGCAGGCATAGGCGCAGGTCATTTCTTCTTCAGGCAATCGCTCATGAAGGCCTTGCGCTCATCCCCCTTCAGGGCTTTGGTGGTGGCGTCAGCATTGCAGGTCTTCATTTTTTCTTGCTGGGTTTGG  >  NC_002947/5518405‑5518698

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