Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP000730	1,996,160	G→A	11.2%	noncoding (55/75 nt)	USA300HOU_t0027 ←	tRNA‑Gly

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP000730	1,996,160	0	G	A	11.2%	167.6 / 8.1	67	noncoding (55/75 nt)	USA300HOU_t0027	tRNA‑Gly
Reads supporting (aligned to +/- strand):  ref base G (34/25);  new base A (6/2);  total (40/27)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

ACTACTCCCGCAAATATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCAATA  >  CP000730/1996026‑1996299                                                                                                                                       |                                                                                                                                            actactCCCGCAAATATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCta                                                                                                                                       >  1:174371/1‑141 (MQ=255)     ctgccGCAATTATTAATATATAGAGCGGAAGGTAGTATTTACACTTGTACCTTATTCCAGGAAGGCATGTATTCTAAGAGTTTCAACACGCCTGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATAcc                                                                                                                                   <  2:591708/138‑1 (MQ=18)                 ttAAATTATGGAGCGGAAGATAGGATTTACACCGATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATCGGATTTGCACCTATACCTCGTTCCGggaa                                                                                                                       >  1:188256/1‑141 (MQ=255)                   aaaaaatggggggaaaGATTTGATATACACCTATACCTCACACCAGGAAGGAATTTATTTTATGAGTAAAAATACTCCCCCATTATTATGAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGggaagg                                                                                                                     <  2:356214/127‑1 (MQ=18)                   aaaTTATGGAGCGGAAGATAGGCTTTACACATATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTccttcccgg                                                                                                                     <  2:522459/141‑8 (MQ=18)                   aaaTTATGGAGCGGAAGATAGGATTTACACCTATACCTCATTCCAGGAAGGAGTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGgggagg                                                                                                                     >  2:393686/1‑141 (MQ=255)                   aaaTAAAGGAGCGGAAGAAAGGAGTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGAAGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGggaagg                                                                                                                     <  2:35893/141‑1 (MQ=21)                             gCGGAAGATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCCTTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTGGTTCCGGGAAGGAACGTGTTCt                                                                                                           >  1:348762/1‑141 (MQ=255)                                aaaGATAGGATATACACCTATACCTCATTCCAGGAAGGAATGGATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTaaa                                                                                                        <  2:738421/140‑1 (MQ=255)                                 aaGATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTGCTaaaa                                                                                                       >  2:599524/1‑141 (MQ=255)                                   gATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCt                                                                                                           <  2:742771/135‑1 (MQ=255)                                   gATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCt                                                                                                           >  1:742771/1‑135 (MQ=255)                                         aTTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTACACCTCGTTCCGGGAAGGAACGTGATCTAAAAGTTGAact                                                                                               >  1:295171/1‑141 (MQ=255)                                            tACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTCTACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAactact                                                                                            >  2:699587/1‑141 (MQ=255)                                               acCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTccc                                                                                         >  1:679978/1‑141 (MQ=255)                                                ccTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACATATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTGGAACTACTCCCg                                                                                        >  1:345683/1‑141 (MQ=255)                                                 cTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGc                                                                                       >  1:464400/1‑141 (MQ=255)                                                     aCCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTATAAAAGTTGAACTACTCCCGCATaa                                                                                   >  1:487078/1‑141 (MQ=255)                                                       cTCGTTACAGGAATGGATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTGAATTCTGGAGTGGAGGGTAGGATTTGCACCTATCCGTCGTTCCGGGAAGGCACGGGGTCTAGAAGTTGAACTACTCCCGCATAAAc                                                                                 >  2:408862/1‑141 (MQ=21)                                                       cTCATTCCAGGAATGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATAAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCGAAAAGTTGAACTCCTTCCGCATAAAc                                                                                 >  1:246064/1‑141 (MQ=25)                                                       cTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTCCTCCCGCCTAAAc                                                                                 >  2:135567/1‑141 (MQ=255)                                                       cTCATTCCAGGAAGGAATGTATTCTAAGAGCTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAAc                                                                                 <  1:248284/141‑1 (MQ=255)                                                           ttCCAGGAAGTAATGTATTCTAAGAGTTAAAATACACCAACATTCTTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCCCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCaaatattaaa                                                                             <  1:340181/141‑10 (MQ=14)                                                            tCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCAAATAccagaa                                                                            >  2:274221/1‑137 (MQ=255)                                                               aggaaggaATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCAAAAACCCggtgag                                                                         >  2:466165/1‑137 (MQ=255)                                                                 gaaggaaTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAggcgg                                                                       <  1:260088/141‑1 (MQ=255)                                                                 gaaggaaTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAggcgg                                                                       <  1:517972/141‑1 (MQ=255)                                                                    ggaaTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTTTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCaa                                                                    <  1:542412/141‑1 (MQ=255)                                                                     gaaTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTTTAAAAGTTGAACTAGTCCCGACTAAACCTGGAGGCGGCaga                                                                   >  2:316546/1‑139 (MQ=255)                                                                     gaaTGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACATCGTTCCGGGAAGGAACGTATTCAAAAAGTTGAACTAATCCCGCATAGACCTTGATGCAGCAAc                                                                   >  2:223805/1‑141 (MQ=21)                                                                           attttcggagTCGACATACACCCGCATTATTATTGAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGAtt                                                             <  2:490073/134‑1 (MQ=21)                                                                             tCTAAGAATTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTg                                                           <  2:174371/141‑1 (MQ=255)                                                                              cTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGa                                                          <  1:466165/141‑1 (MQ=255)                                                                                    gTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGGAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCGACCGGATTTGAACCggg                                                    >  2:386209/1‑140 (MQ=255)                                                                                    gTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTGCTAAAAATTGAACTACTCCCGCATAAACCTGGACGCGGCAACCGGATTTGAACCgcg                                                    >  2:516955/1‑139 (MQ=255)                                                                                      tGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGa                                                  <  1:601966/141‑1 (MQ=255)                                                                                        aaaTACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATa                                                <  2:464400/141‑1 (MQ=255)                                                                                          aTACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAGCCGGAtttttcccggctataaa                                              >  2:541782/1‑127 (MQ=21)                                                                                          aTACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATaaa                                              >  2:410300/1‑141 (MQ=255)                                                                                               cccGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGtt                                                                                                    <  1:398775/82‑1 (MQ=255)                                                                                               cccGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGtt                                                                                                    <  1:626042/82‑1 (MQ=255)                                                                                               cccGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGtt                                                                                                    >  2:398775/1‑82 (MQ=255)                                                                                               cccGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGtt                                                                                                    >  2:626042/1‑82 (MQ=255)                                                                                                  gCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCCGTGATAAAGGTTTTGc                                      <  1:394766/141‑1 (MQ=255)                                                                                                    attattattAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGGGATAAAGGTTTTGCAg                                    >  1:362733/1‑141 (MQ=255)                                                                                                       attattaAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGACGGAACGGGTTCTAAGCGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGAATGTGCACCGGTGATAGAGGGTTTGAAGAcg                                 >  1:666541/1‑140 (MQ=21)                                                                                                       attattaAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGAcc                                 >  1:666502/1‑141 (MQ=255)                                                                                                         tattaAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACctc                               >  2:246050/1‑141 (MQ=255)                                                                                                         tattaAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACctc                               <  1:245479/141‑1 (MQ=255)                                                                                                         tattaAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACctc                               <  1:654896/141‑1 (MQ=255)                                                                                                              aaTTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCt                          <  1:274221/141‑1 (MQ=255)                                                                                                                       gCGGAAGACAGTGTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTg                 <  2:377021/141‑1 (MQ=255)                                                                                                                        cGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTgg                <  1:120068/141‑1 (MQ=255)                                                                                                                         ggAAGAGAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGc               <  2:679978/141‑1 (MQ=255)                                                                                                                           aaGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCCCTTGGCtt             >  1:173480/1‑140 (MQ=255)                                                                                                                           aaGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTa             >  1:157382/1‑141 (MQ=255)                                                                                                                           aaGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTa             >  1:173425/1‑141 (MQ=255)                                                                                                                             gATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATCAAGGTTTTGCAGACCTCTTCCTTACCACTTGGCTATg           >  1:226091/1‑141 (MQ=255)                                                                                                                                aGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATgccg        >  1:328374/1‑141 (MQ=255)                                                                                                                                aGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATgccg        <  1:330989/141‑1 (MQ=255)                                                                                                                                  gATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATgccgcc      <  1:749343/141‑1 (MQ=255)                                                                                                                                   aTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACGTCGGACTTACCACTTGGCGATGCCGCCa     >  2:290889/1‑141 (MQ=255)                                                                                                                                   aTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCa     >  1:602287/1‑141 (MQ=255)                                                                                                                                   aTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCa     >  1:144946/1‑141 (MQ=255)                                                                                                                                   aTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGACAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCa     >  1:728458/1‑141 (MQ=255)                                                                                                                                   aTTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTACCGCGTCAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACGTGTCTATGCCCCCa     >  2:283989/1‑141 (MQ=37)                                                                                                                                      tGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCAATa  >  1:137518/1‑141 (MQ=255)                                                                                                                                       |                                                                                                                                            ACTACTCCCGCAAATATTAAATTATGGAGCGGAAGATAGGATTTACACCTATACCTCATTCCAGGAAGGAATGTATTCTAAGAGTTGAAATACTCCCGCATTATTATTAAATTATGGAGCGGAAGATAGGATTTGCACCTATACCTCGTTCCGGGAAGGAACGTGTTCTAAAAGTTGAACTACTCCCGCATAAACCTGGAGGCGGCAACCGGATTTGAACCGGTGATAAAGGTTTTGCAGACCTCTGCCTTACCACTTGGCTATGCCGCCAATA  >  CP000730/1996026‑1996299

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