Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,931,602	C→A	5.2%	P241Q (CCG→CAG)	purT →	phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,931,602	0	C	A	5.2%	281.3 / 2.2	77	P241Q (CCG→CAG)	purT	phosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base C (24/49);  new base A (2/2);  total (26/51)
Fisher's exact test for biased strand distribution p-value = 6.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

AGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACAT  >  NC_000913/1931465‑1931735                                                                                                                                          |                                                                                                                                      aGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGc                                                                                                                                     <  1:930487/140‑1 (MQ=255) aGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGc                                                                                                                                     <  1:930285/140‑1 (MQ=255)   gCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCACATGAGCCCGCtt                                                                                                                                   <  1:737144/140‑1 (MQ=255)   gCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCtt                                                                                                                                   <  1:1611556/140‑1 (MQ=255)   gCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCtt                                                                                                                                   <  1:827188/140‑1 (MQ=255)       cgtTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGccc                                                                                                                               >  1:640068/1‑140 (MQ=255)       cgtTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGccc                                                                                                                               >  1:1677871/1‑140 (MQ=255)            aGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGaa                                                                                                                          >  1:154006/1‑140 (MQ=255)             gattGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAAc                                                                                                                         <  2:928450/138‑1 (MQ=255)             gTTTGACTTCGAAATTACCCTGCTCACCGTCAGCGCGGTGGATGGCGTCCATTTCTGGGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAAc                                                                                                                         <  2:1270695/140‑1 (MQ=255)             gTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAAc                                                                                                                         <  2:928474/140‑1 (MQ=255)             gTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAAc                                                                                                                         <  2:855985/140‑1 (MQ=255)              tttGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACg                                                                                                                        >  1:1627146/1‑140 (MQ=255)                tGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTg                                                                                                                      >  1:1439233/1‑140 (MQ=255)                tGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTg                                                                                                                      >  1:1556821/1‑140 (MQ=255)                tGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTg                                                                                                                      >  1:1557282/1‑140 (MQ=255)                  aCTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTgcg                                                                                                                    >  2:1556613/1‑140 (MQ=255)                          aTTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGAt                                                                                                            <  2:528067/140‑1 (MQ=255)                                      aCCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt                                                                                                <  2:489250/140‑1 (MQ=255)                                      aCCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt                                                                                                <  2:1215010/140‑1 (MQ=255)                                      aCCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt                                                                                                <  2:489192/140‑1 (MQ=255)                                      aCCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAgt                                                                                                <  1:702989/140‑1 (MQ=255)                                                cgGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCa                                                                                      <  2:499530/140‑1 (MQ=255)                                                  gTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACt                                                                                    >  2:1328058/1‑140 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  2:641939/140‑1 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  2:217331/140‑1 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  2:837245/140‑1 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  1:1461010/140‑1 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  2:49325/140‑1 (MQ=255)                                                     gATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTggg                                                                                 <  1:789496/140‑1 (MQ=255)                                                               aTTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATggg                                                                       <  1:1371334/140‑1 (MQ=255)                                                               aTTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATggg                                                                       <  1:1370278/140‑1 (MQ=255)                                                               aTTTCTGTGCACAAGTAGGTCATCGCCAGGAAGTTGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCACGCTTGCCCTTGATCGTGCGCAGGTGATTGCCCGTAAAGTGGTGATGGCACTGGGCGGTTATggg                                                                       <  1:1233162/140‑1 (MQ=255)                                                                      tgCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTg                                                                >  1:502924/1‑140 (MQ=255)                                                                      tgCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTg                                                                >  1:503044/1‑140 (MQ=255)                                                                        cACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGgt                                                              >  2:1035565/1‑140 (MQ=255)                                                                        cACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGGAAAGTGGTGCTGGCACTGGGCGGGTAGGGGGGGGTTGgt                                                              >  1:1636230/1‑140 (MQ=255)                                                                            aGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCg                                                          <  2:768764/140‑1 (MQ=255)                                                                             gTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTCAAGTGGTGCTGGCACTGGGCGGTTCTGGGTTGTTTGGTGTCGa                                                         >  1:1104787/1‑140 (MQ=255)                                                                              tAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTTTGGGTTGTTTGGTGTCGAg                                                        <  1:1076756/140‑1 (MQ=255)                                                                              tAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAg                                                        <  1:651354/140‑1 (MQ=255)                                                                              tAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAg                                                        <  1:17464/140‑1 (MQ=255)                                                                                ggTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCAGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCg                                                                                                      <  2:1017048/92‑1 (MQ=255)                                                                                ggTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCAGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCg                                                                                                      >  1:1017048/1‑92 (MQ=255)                                                                                ggTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCAGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCg                                                                                                      <  2:1016320/92‑1 (MQ=255)                                                                                ggTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCAGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCg                                                                                                      >  1:1016320/1‑92 (MQ=255)                                                                                 gTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTa                                                     <  2:129375/140‑1 (MQ=255)                                                                                 gTCATCGACAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTa                                                     <  2:746990/140‑1 (MQ=255)                                                                                  tCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTAt                                                    <  1:674592/140‑1 (MQ=255)                                                                                  tCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTAt                                                    <  1:674246/140‑1 (MQ=255)                                                                                     tCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTCTTTg                                                 >  1:1685674/1‑140 (MQ=255)                                                                                        ccAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCt                                              >  2:268980/1‑140 (MQ=255)                                                                                         cAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCtg                                             <  2:1071897/140‑1 (MQ=255)                                                                                                aTGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGGTGTTTGGTGTCGAGCTATTTGTCTGTGGtgat                                      >  1:787683/1‑140 (MQ=255)                                                                                                     gACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGt                                 <  2:1439233/140‑1 (MQ=255)                                                                                                        tACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGAt                              <  2:116755/140‑1 (MQ=255)                                                                                                        tACCGAGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGAt                              <  2:856852/140‑1 (MQ=255)                                                                                                            gTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTc                          <  1:214374/140‑1 (MQ=255)                                                                                                                  ccTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAg                    >  2:169018/1‑140 (MQ=255)                                                                                                                  ccTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAg                    <  1:1294524/140‑1 (MQ=255)                                                                                                                  ccTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAg                    <  1:1294123/140‑1 (MQ=255)                                                                                                                    tGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGt                  >  2:1305787/1‑140 (MQ=255)                                                                                                                    tGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGt                  >  2:1305342/1‑140 (MQ=255)                                                                                                                    tGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGt                  >  1:1070791/1‑140 (MQ=255)                                                                                                                    tGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGt                  >  2:803123/1‑140 (MQ=255)                                                                                                                    tGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTCATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGt                  <  1:342183/140‑1 (MQ=255)                                                                                                                          ccACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTcccc            <  1:1064433/140‑1 (MQ=255)                                                                                                                          ccACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGGGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTcccc            <  1:1063418/140‑1 (MQ=255)                                                                                                                           cacaGCAAATGAGCCCGCGTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTccccg           >  2:1531351/1‑139 (MQ=255)                                                                                                                            acaGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTc          >  1:30774/1‑140 (MQ=255)                                                                                                                              aGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGGGGTGATTTTCAGTGAGGTCTCCCCTCGt        <  1:1380654/140‑1 (MQ=255)                                                                                                                              aGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGGGGTGATTTTCAGTGAGGTCTCCCCTCGt        <  1:1380536/140‑1 (MQ=255)                                                                                                                              aGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGt        <  1:1109820/140‑1 (MQ=255)                                                                                                                              aGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGt        <  1:1475119/140‑1 (MQ=255)                                                                                                                                   aTGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCcaca   <  1:1135236/140‑1 (MQ=255)                                                                                                                                   aTGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCcaca   <  1:711505/140‑1 (MQ=255)                                                                                                                                    tGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGCCTCCCCTCGTCCACAt  <  2:1424004/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      AGGCGTCGTTAAGTTTGACTTCGAAATTACCCTGCTAACCGTCAGCGCGGTGGATGGCGTCCATTTCTGTGCACCAGTAGGTCATCGCCAGGAAGATGGCGACTACCGTGAATCCTGGCAACCACAGCAAATGAGCCCGCTTGCCCTTGAACGTGCGCAGGAGATTGCCCGTAAAGTGGTGCTGGCACTGGGCGGTTATGGGTTGTTTGGTGTCGAGCTATTTGTCTGTGGTGATGAGGTGATTTTCAGTGAGGTCTCCCCTCGTCCACAT  >  NC_000913/1931465‑1931735

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