Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	955,586	C→G	6.2%	G75G (GGC→GGG)	alaC →	aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	955,586	0	C	G	6.2%	202.0 / 2.4	64	G75G (GGC→GGG)	alaC	aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base G (2/2);  total (29/35)
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

CTCAAGATGGCCGCCCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGC  >  NC_002947/955440‑955727                                                                                                                                                   |                                                                                                                                              cTCAAGATGGCCGCCCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATc                                                                                                                                            >  1:1212511/1‑150 (MQ=255)       aTGGCCGCCCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCTATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCgcgc                                                                                                                                      <  1:1275300/150‑1 (MQ=255)       aTGGCCGCCCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCgcgc                                                                                                                                      >  1:510097/1‑150 (MQ=255)          gccgccCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGCCGTATCCACGTCGCGCGGGATGCCgcgcgac                                                                                                                                   <  1:1198431/150‑4 (MQ=255)            cgcccgccGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACCACGGCCTGTAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTgcg                                                                                                                                 >  1:607616/1‑150 (MQ=255)               ccgccgccgTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCg                                                                                                                              >  2:912110/1‑150 (MQ=255)                ctcctccgTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCgg                                                                                                                             <  1:1089846/145‑1 (MQ=255)                cgccgccgTGGCGAGGATATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGACGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCgg                                                                                                                             <  2:567097/150‑1 (MQ=255)                   cgccgTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGcc                                                                                                                          <  1:912110/150‑1 (MQ=255)                        tGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTTGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCAtctc                                                                                                                     <  2:98665/150‑1 (MQ=255)                           cGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCaa                                                                                                                  <  1:299199/150‑1 (MQ=255)                               gACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTgg                                                                                                              >  1:77929/1‑150 (MQ=255)                                aCATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGt                                                                                                             <  1:219492/150‑1 (MQ=255)                                  atcatcGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTAc                                                                                                           >  1:361067/1‑150 (MQ=255)                                     atcGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAg                                                                                                        <  1:900969/150‑1 (MQ=255)                                     atcGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACAGGGACAAg                                                                                                        >  1:54845/1‑150 (MQ=255)                                              aGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGGCATCTCCAACTGGTACAAGGACCGCTAc                                                                                               >  1:644031/1‑150 (MQ=255)                                              aGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTAc                                                                                               <  1:862459/150‑1 (MQ=255)                                              aGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTAc                                                                                               <  1:49744/150‑1 (MQ=255)                                                 aTGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAg                                                                                            <  1:564800/150‑1 (MQ=255)                                                  tGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGACTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAgg                                                                                           <  1:129768/150‑1 (MQ=255)                                                       aaCCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGAc                                                                                      <  2:96603/150‑1 (MQ=255)                                                        aCCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACa                                                                                     <  2:510097/150‑1 (MQ=255)                                                           ccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCAtct                                                                                                                      >  2:685053/1‑114 (MQ=255)                                                           ccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCAtct                                                                                                                      <  1:685053/114‑1 (MQ=255)                                                           ccGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCg                                                                                  >  2:1101827/1‑150 (MQ=255)                                                                   gcCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGaa                                                                          >  2:249526/1‑150 (MQ=255)                                                                    ccACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTCCTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGaaa                                                                         <  1:774635/150‑1 (MQ=255)                                                                        cccGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACaa                                                                                                         <  1:420497/114‑1 (MQ=255)                                                                        cccGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACaa                                                                                                         >  2:420497/1‑114 (MQ=255)                                                                             cgcACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCa                                                                <  1:249526/150‑1 (MQ=255)                                                                               cacaTCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGATACGAAAATGAAATCGACCCGGAAAGCGAAGCCATc                                                              <  1:354147/150‑1 (MQ=255)                                                                                       cgAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCAt                                                      <  2:670828/150‑1 (MQ=255)                                                                                        gagaAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTTACCATc                                                     >  1:147244/1‑150 (MQ=255)                                                                                                    cAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCCTCGTCACCATCGGCTCGAAagag                                         >  1:245892/1‑150 (MQ=255)                                                                                                     aGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAgg                                        >  1:669070/1‑150 (MQ=255)                                                                                                     aGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAgg                                        <  2:607616/150‑1 (MQ=255)                                                                                                     aGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGACTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAgg                                        <  2:1159006/150‑1 (MQ=255)                                                                                                            ccAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGc                                 <  2:378344/150‑1 (MQ=255)                                                                                                            ccAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGc                                 >  1:1050665/1‑150 (MQ=255)                                                                                                            ccAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGc                                 >  1:378344/1‑150 (MQ=255)                                                                                                                  tGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGaaa                                                                         >  1:25768/1‑104 (MQ=255)                                                                                                                  tGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGaaa                                                                         <  2:25768/104‑1 (MQ=255)                                                                                                                     aGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGAt                        >  1:759698/1‑150 (MQ=255)                                                                                                                         acCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCCCCATCCGCGCCAAAGAGGGCCTGGCCCACCTGATGCTg                    >  2:1066324/1‑150 (MQ=255)                                                                                                                          cccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTgg                   <  1:1101827/150‑1 (MQ=255)                                                                                                                           ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCGTCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGc                  >  1:375612/1‑150 (MQ=255)                                                                                                                           ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAgg                                        >  2:942431/1‑128 (MQ=255)                                                                                                                           ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAgg                                        <  1:942431/128‑1 (MQ=255)                                                                                                                           ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCAc                             >  1:806843/1‑139 (MQ=255)                                                                                                                           ccACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCAc                             <  2:806843/139‑1 (MQ=255)                                                                                                                                gCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCAccc             <  2:361067/150‑1 (MQ=255)                                                                                                                                  tACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAACGAGGGCCTGGCCCACCTGATGCTGGCCACCCTc           <  2:77929/150‑1 (MQ=255)                                                                                                                                     tCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCAGGAAAGCGAAGCCATCGCCCCCATCGACCCAAAAGAGGGCCTTGCCCACCGGATGCTGGCCAACCtagcc        >  2:530948/1‑146 (MQ=255)                                                                                                                                      ccACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGTCCCACCTGATGCTGGCCACCCTCGAcc       >  1:965421/1‑150 (MQ=255)                                                                                                                                        aCGTCGCGCGGCATCCCGCGCCTGCGACGGGCCATCTCCAACTGGTACAATGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAg     >  2:979182/1‑150 (MQ=255)                                                                                                                                          gTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCCTCGGCACGAAAGAGGGCCTGGTCCACCTGATGCTGGCCACCCTCGACCAggg   <  2:147244/150‑1 (MQ=255)                                                                                                                                          gTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAggg   <  2:54845/150‑1 (MQ=255)                                                                                                                                          gTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAggg   <  2:469568/150‑1 (MQ=255)                                                                                                                                          caggcgcgGGATGCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTc           >  1:1266716/4‑142 (MQ=255)                                                                                                                                          caggcgcgGGATGCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTc           <  2:1266716/139‑1 (MQ=255)                                                                                                                                           aggcgcgGGATGCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAAGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGc  >  1:476489/3‑150 (MQ=255)                                                                                                                                                 ggCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTc                                                >  1:1230883/1‑98 (MQ=255)                                                                                                                                                 ggCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTc                                                <  2:1230883/98‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                              CTCAAGATGGCCGCCCGCCGCCGTGGCGAGGACATCATCGACCTGAGCATGGGCAACCCCGATGGCGCCACCCCGCCGCACATCGTCGAGAAGCTGGTACAGGTCGCCCAGCGTGAAGACACCCACGGCTACTCCACGTCGCGCGGCATCCCGCGCCTGCGCCGGGCCATCTCCAACTGGTACAAGGACCGCTACGAGGTCGACATCGACCCGGAAAGCGAAGCCATCGTCACCATCGGCTCGAAAGAGGGCCTGGCCCACCTGATGCTGGCCACCCTCGACCAGGGC  >  NC_002947/955440‑955727

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