Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	992,868	A→T	7.4%	intergenic (+23/‑129)	der → / → PP_0858	GTPase Der/methionine/glutamine aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	992,868	0	A	T	7.4%	192.6 / 4.1	66	intergenic (+23/‑129)	der/PP_0858	GTPase Der/methionine/glutamine aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base T (2/3);  total (36/30)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCGCAGCAAACTG  >  NC_002947/992722‑993014                                                                                                                                                   |                                                                                                                                                   cAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGCCAAACGCTAGGCAGCGCCTTGTGGGGGCGGGCAAGa                                                                                                                                                 <  2:1504526/150‑1 (MQ=255)        ggcgAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGCAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCCTGTGAGGGCGGTCAAGTCCGcccc                                                                                                                                          <  2:1179117/150‑1 (MQ=255)         gcgAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGccccc                                                                                                                                         <  1:399169/150‑1 (MQ=255)         gcgAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCGAGACCGccccc                                                                                                                                         >  2:561414/1‑150 (MQ=255)                cccGTACGAGAGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCcacatggc                                                                                                                                  <  2:1402759/150‑5 (MQ=255)                    tACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCt                                                                                                                              <  1:487838/150‑1 (MQ=255)                        aGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGa                                                                                                                          <  2:960559/150‑1 (MQ=255)                             aagaagAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCag                                                                                                                     >  1:424751/1‑150 (MQ=255)                              agaagaACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                                    >  1:868698/1‑150 (MQ=255)                               gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCGATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                                    >  2:1097953/1‑149 (MQ=255)                               gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCGAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   <  1:561414/150‑1 (MQ=255)                               gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   >  1:607091/1‑150 (MQ=255)                               gaagaaCACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCaga                                                                                                                    <  1:1097953/149‑1 (MQ=255)                                        ccTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATccccc                                                                                                          >  2:1284018/1‑150 (MQ=255)                                          ttACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGccccc                                                                                                                                         >  1:717544/1‑117 (MQ=255)                                          ttACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGccccc                                                                                                                                         <  2:717544/117‑1 (MQ=255)                                          ttACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGTAATCCCCCTa                                                                                                        >  2:1537025/1‑150 (MQ=255)                                            aCCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTAtt                                                                                                      >  1:597333/1‑150 (MQ=255)                                               gATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGtt                                                                                                   >  2:377270/1‑150 (MQ=255)                                                aTCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGttt                                                                                                  >  2:285825/1‑150 (MQ=255)                                                  cGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCa                                                                                                <  2:597333/150‑1 (MQ=255)                                                    ccAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAAc                                                                                              >  2:1117559/1‑150 (MQ=255)                                                    ccAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGAGCCAGAGCAATCCCCCTATTGTTTCAAc                                                                                              >  2:1440997/1‑150 (MQ=255)                                                      aGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCcacatg                                                                                                                                    >  2:999656/1‑108 (MQ=39)                                                      aGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGTCTTGCCCGCCCCcacatg                                                                                                                                    <  1:999656/110‑3 (MQ=38)                                                        gTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCttt                                                                                          <  1:1284018/150‑1 (MQ=255)                                                             caaGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGTAATCCCCCTATTGTTTCAACCTTTTTCCt                                                                                     <  1:1537025/150‑1 (MQ=255)                                                              aaGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATcc                                                                                                                       <  2:219458/115‑1 (MQ=255)                                                              aaGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATcc                                                                                                                       >  1:219458/1‑115 (MQ=255)                                                                    cgccgcCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCt                                                                              <  1:525414/150‑1 (MQ=255)                                                                        gccTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGAt                                                                          >  1:478938/1‑150 (MQ=255)                                                                          cTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   >  2:243403/1‑107 (MQ=255)                                                                          cTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   <  1:243403/107‑1 (MQ=255)                                                                           tGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGATAAACGCTAAGCAGCGCCATGTGGGGTCGGGCCAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGAcccc                                                                       >  1:250531/1‑149 (MQ=255)                                                                               gTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   <  2:1039317/102‑1 (MQ=255)                                                                               gTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCagag                                                                                                                   >  1:1039317/1‑102 (MQ=255)                                                                               gTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTAt                                                                   >  1:787484/1‑150 (MQ=255)                                                                                      cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAAc                                                                                              >  1:1555299/1‑116 (MQ=255)                                                                                      cacaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAAc                                                                                              <  2:1555299/116‑1 (MQ=255)                                                                                       acaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATccc                                                                                                            >  1:1340489/1‑101 (MQ=255)                                                                                       acaAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATccc                                                                                                            <  2:1340489/101‑1 (MQ=255)                                                                                        caAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACt                                                          >  2:1242999/1‑150 (MQ=255)                                                                                         aagaagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTc                                                         >  1:939939/1‑150 (MQ=255)                                                                                           gaagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCa                                                       <  1:1242999/150‑1 (MQ=255)                                                                                            aagGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAc                                                      >  2:1613464/1‑150 (MQ=255)                                                                                             agGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAcc                                                     <  1:1613464/150‑1 (MQ=255)                                                                                                 cGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGAcc                                                                                 <  2:598344/118‑1 (MQ=255)                                                                                                 cGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGAcc                                                                                 >  1:598344/1‑118 (MQ=255)                                                                                                       gaagCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAccc                                                    >  2:1548743/1‑141 (MQ=255)                                                                                                       gaagCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCAccc                                                    <  1:1548743/141‑1 (MQ=255)                                                                                                                      aCGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAag                            >  1:1090412/1‑150 (MQ=255)                                                                                                                        gCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAagag                          >  1:781888/1‑150 (MQ=255)                                                                                                                         cTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAgg                         <  1:377270/150‑1 (MQ=255)                                                                                                                          tAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGat                                          >  1:1081255/1‑132 (MQ=255)                                                                                                                          tAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGat                                          <  2:1081255/132‑1 (MQ=255)                                                                                                                            agcagcGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCACCTATTGt                                                                                                    <  1:908186/72‑1 (MQ=255)                                                                                                                            agcagcGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCACCTATTGt                                                                                                    >  2:908186/1‑72 (MQ=255)                                                                                                                                    cATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCa                                  >  2:238553/1‑130 (MQ=255)                                                                                                                                    cATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCa                                  <  1:238553/130‑1 (MQ=255)                                                                                                                                      tgtgGGGGCGGTCTTGCCCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCg            >  1:118211/1‑150 (MQ=255)                                                                                                                                        tgGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCgca          <  1:285825/150‑1 (MQ=255)                                                                                                                                         gggggCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCgcag         <  2:787484/150‑1 (MQ=255)                                                                                                                                          ggggCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCgcagc        >  2:856741/1‑150 (MQ=255)                                                                                                                                                ggcAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTACCGTGGAATAGGGTTGCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACTCCAGTAAAGAGGGCTCCATGATCCGCAGCAAACTg  >  2:558352/1‑150 (MQ=255)                                                                                                                                                  cAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCt                                                                              >  1:95171/1‑72 (MQ=255)                                                                                                                                                  cAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCt                                                                              <  2:95171/72‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CAAGGGCGGCGAGAACCCGTACGAGGGTAAGAAGAACACCCTTACCGATCGCCAGGTCAACAAGAAACGCCGCCTGATGTCGCACCACAAGAAGGCCGAGAAGAAGCGTCGCGACAAACGCTAAGCAGCGCCATGTGGGGGCGGGCAAGACCGCCCCCACAGACGCCCTTTGATCCAGAGCAATCCCCCTATTGTTTCAACCTTTTTCCTGACCGCTCGATCCGCTATGCTCGGACTCCACCCCGTGCCGGATATACCCCAGGAAAGAGGGCTCCATGATCCGCAGCAAACTG  >  NC_002947/992722‑993014

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