Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,799,301	A→G	8.2%	A330A (GCT→GCC) ‡	ansB ←	glutaminase‑asparaginase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,799,301	0	A	G	8.2%	168.8 / 5.2	61	A330A (GCT→GCC) ‡	ansB	glutaminase‑asparaginase
Reads supporting (aligned to +/- strand):  ref base A (27/29);  new base G (2/3);  total (29/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 = 6.10e-01

GGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACCGT  >  NC_002947/2799154‑2799449                                                                                                                                                       |                                                                                                                                                     ggCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGa                                                                                                                                                      <  2:354963/148‑1 (MQ=255)         ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             <  3:347424/149‑1 (MQ=255)         ggCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             <  2:307609/149‑1 (MQ=255)          gCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcct                                                                                                                                            <  8:100153/149‑2 (MQ=255)           cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAcc                                                                                                                                             >  3:350471/1‑147 (MQ=255)           cccTCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAc                                                                                                                                              >  7:165619/1‑146 (MQ=255)           cccTCAGGGCCGGTGCGCCGGCGCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGAc                                                                                                                                              <  8:165619/146‑1 (MQ=255)              tCAGGGCCGGTGCGCCGGCTCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGCGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        >  1:127270/1‑149 (MQ=255)              tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTCCCATCGCCAGGATCCGCGCTTGTTGCGGGTTCCTGTCGTGAGCCCCGCCCCGGTc                                                                                                                                        >  4:310220/1‑149 (MQ=255)              tCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        >  1:218856/1‑149 (MQ=255)                agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctgaacc                                                                                                                                       <  1:147810/148‑7 (MQ=255)                agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGGCTCACGAcctgaacc                                                                                                                                       <  3:184753/148‑7 (MQ=255)                agagCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTc                                                                                                                                        <  6:45834/147‑1 (MQ=255)                  agCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTc                                                                                                                                    <  4:251285/149‑1 (MQ=255)                       gcgcACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCATACTCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  <  1:126033/149‑1 (MQ=255)                       gcgcACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTc                                                                                                                               <  8:6304/149‑1 (MQ=255)                            ccGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCg                                                                                                                          >  6:393763/1‑149 (MQ=255)                                  cTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGtt                                                                                                                    <  3:334203/149‑1 (MQ=255)                                            aGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACg                                                                                                          <  2:287445/149‑1 (MQ=255)                                              cACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCa                                                                                                        <  2:59035/149‑1 (MQ=255)                                                 tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTa                                                                                                     >  4:274304/1‑149 (MQ=255)                                                 tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTa                                                                                                     >  1:65935/1‑149 (MQ=255)                                                 tCAGTACTCCCAGAAGATCCGCTGCAGCTCTTAGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCAACGCCAGGATACTCCCGTATTGCGGGATCACGACGTGAGCCACGCCCCAGTCGTTCTTGTAATCCGGATCTACGCCAATACGCAGTa                                                                                                     >  5:13162/1‑149 (MQ=255)                                                   aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTTAGCCACGACCCAGTCGTTCTTGTCATACGGCTGTTCGGCATTACGCa                                                                                                        <  6:173918/144‑1 (MQ=255)                                                   aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTAc                                                                                                    >  2:251122/1‑148 (MQ=255)                                                   aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACg                                                                                                   <  2:218856/149‑1 (MQ=255)                                                   aGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTAATCCGGCTGTTCGGCATTACGCa                                                                                                        >  5:173918/1‑144 (MQ=255)                                                       cTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAac                                                                                               >  6:352143/1‑149 (MQ=255)                                                           cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  <  3:225884/110‑1 (MQ=255)                                                           cAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCtt                                                                                                                                  >  4:225884/1‑110 (MQ=255)                                                              aagaTCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGACACGACCCAGTCGTTATTGTCGTCCGGCTGTTCGGCATTACGCAGTACGAAACCAcccgg                                                                                        >  8:149133/1‑147 (MQ=255)                                                                      cTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACAt                                                                                <  3:171228/149‑1 (MQ=255)                                                                          aGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGcga                                                                            <  5:86852/149‑1 (MQ=255)                                                                                 tGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCg                                                                      >  3:176315/1‑148 (MQ=255)                                                                                   cTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGga                                                                    <  5:352143/148‑1 (MQ=255)                                                                                   cTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGACATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGgat                                                                   >  4:202949/1‑149 (MQ=255)                                                                                      tCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGgatgat                                                                <  4:125905/149‑1 (MQ=255)                                                                                        cTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATAAGGCACTACGAACCCACCCGGGGTGACATGCGCCGACCGGATGctct                                                              >  4:98628/1‑149 (MQ=255)                                                                                           ggTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGc                                                                                                                        >  6:176665/1‑88 (MQ=255)                                                                                           ggTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGc                                                                                                                        <  5:176665/88‑1 (MQ=255)                                                                                             gcTTGGGCA‑‑‑TCGCTACCATCGCCAGTATCAGAGCTTTTTGAGGTTTCAGGTCGTGAGCCACGACCCAGTAGGTCTTGTAATCCTGCTTTTCGGCATTACTAAGTACGAAACCACACTTTTTGACATGCGACGAGCTGATGATCTGAACg                                                         <  6:77500/148‑1 (MQ=255)                                                                                                     a‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCtt                                                 <  2:127270/149‑1 (MQ=255)                                                                                                     a‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCtt                                                 <  4:176936/149‑1 (MQ=255)                                                                                                          cGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTgc                                               >  2:153853/1‑149 (MQ=255)                                                                                                           gCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTgc                                               <  4:350471/148‑1 (MQ=255)                                                                                                               ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATc                                                               >  8:54455/1‑128 (MQ=255)                                                                                                               ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATc                                                               <  7:54455/128‑1 (MQ=255)                                                                                                               ccATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCa                                            <  6:161755/147‑1 (MQ=255)                                                                                                                cATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGtt                                         <  2:65935/149‑1 (MQ=255)                                                                                                                 aTCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTc                                        >  5:370269/1‑149 (MQ=255)                                                                                                                           tCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGcc                              <  3:310220/149‑1 (MQ=255)                                                                                                                             cgcgCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCgg                                                                     <  8:82996/108‑1 (MQ=255)                                                                                                                             cgcgCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCgg                                                                     >  7:82996/1‑108 (MQ=255)                                                                                                                              gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  3:216295/1‑148 (MQ=255)                                                                                                                              gcgcTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  5:34645/1‑148 (MQ=255)                                                                                                                                             tCAGGTCGTGAGCCACGACCCAGTCGTTATTGTCATCCGGCTGTTCGGCATTACTCAGTACGACAACACCCTGGTTGACATGCGACGCGCGGATGATCTGAACGCCGTTCTTGGGCAGTTCTTGCAGCGCCGGGACAAGCCGcgccga             >  4:10237/1‑148 (MQ=255)                                                                                                                                              tgggTCGTGGCTCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATa           >  1:305971/3‑149 (MQ=255)                                                                                                                                               aGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAc          >  1:374120/1‑149 (MQ=255)                                                                                                                                               aGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAc          <  5:393763/149‑1 (MQ=255)                                                                                                                                                ggTCGTGGCTCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGGATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATAcc         >  1:286200/1‑149 (MQ=255)                                                                                                                                                    gTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            <  6:280702/126‑1 (MQ=255)                                                                                                                                                    gTGAGCCACGACCCAGTCGTTCTTGTCATACGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCgg                            >  5:280702/1‑126 (MQ=255)                                                                                                                                                       aGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACACTGGTTGACATGCGAAGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACCATACCGAACCGt  >  5:38445/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                     GGCAACGAGGCCCTCAGAGCCGGCGCACCGGCCCTGAGGGGGAAGCACTCAGTACTCCCAGAAGATCCGCTGCAGCTCTTTGCTGTCCTGGGTCTTGGTCA‑‑‑TCGCTACCATCGCCAGGATCCGCGCTTTTTGCGGGTTCAGGTCGTGAGCCACGACCCAGTCGTTCTTGTCATCCGGCTGTTCGGCATTACGCAGTACGAAACCACCCTGGTTGACATGCGACGAGCGGATGATCTGCACGCCGTTCTTGCGCAGTTCTTGCAGCGCCGGGACTACCCGCGACGATACCGAACCGT  >  NC_002947/2799154‑2799449

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