Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,448,947	A→C	62.4%	T477P (ACG→CCG)	feaB →	phenylacetaldehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,448,947	0	A	C	62.4%	89.1 / 95.5	90	T477P (ACG→CCG)	feaB	phenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (13/21);  new base C (22/34);  total (35/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAA  >  NC_000913/1448848‑1449042                                                                                                    |                                                                                                ctGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGaa                                                                                                 >  2:10823/1‑100 (MQ=255)   ggAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCg                                                                                               >  2:2994880/1‑100 (MQ=255)         taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAGACAGCCATACCTTAATTGACGCTAACGTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGt                                                                                         <  1:2029726/100‑1 (MQ=255)         taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACATCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGt                                                                                         >  1:2088992/1‑100 (MQ=255)         taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATGCCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                                                                         <  2:2628644/100‑1 (MQ=255)         taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                                                                         <  2:237638/100‑1 (MQ=255)          aGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCg                                                                                          >  1:100601/1‑98 (MQ=255)           gCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGa                                                                                       >  1:184415/1‑100 (MQ=255)              aTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAttt                                                                                    <  2:428322/100‑1 (MQ=255)                 gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTgg                                                                                 <  2:26703/100‑1 (MQ=255)                 gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTgg                                                                                 <  1:3313405/100‑1 (MQ=255)                 gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTgg                                                                                 >  2:3213441/1‑100 (MQ=255)                  cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGc                                                                                <  2:2409704/100‑1 (MQ=255)                  cTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGc                                                                                >  2:615283/1‑100 (MQ=255)                   ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGcc                                                                               <  2:2486477/100‑1 (MQ=255)                   ttACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcc                                                                               <  2:361952/100‑1 (MQ=255)                    tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGccc                                                                              <  2:1592513/100‑1 (MQ=255)                       aggcaggGACGGTGTGGGTAAACAGCCATACCTTAACTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGa                                                                           >  2:2384181/1‑100 (MQ=255)                        ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAGCTTACCGTTTGGTGGGATGAGGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGa                                                                           >  1:3332559/1‑99 (MQ=255)                        ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                          >  2:1772430/1‑100 (MQ=255)                        ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                          >  2:2639677/1‑100 (MQ=255)                        ggcaggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGATTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGAc                                                                          >  1:3060994/1‑100 (MQ=255)                          caggGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActg                                                                        >  2:1055456/1‑100 (MQ=255)                             ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                     <  1:1711187/100‑1 (MQ=255)                             ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                     <  2:929626/100‑1 (MQ=255)                              gACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggctg                                                                    <  2:622391/100‑1 (MQ=255)                               aCGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGActggct                                                                     <  1:1815997/98‑1 (MQ=255)                                 ggTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                 <  1:2071349/100‑1 (MQ=255)                                  gtgtGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                <  1:1520631/100‑1 (MQ=255)                                  gtgtGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGCTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                <  1:2380149/100‑1 (MQ=255)                                      gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                            >  2:3258430/1‑100 (MQ=255)                                      gggTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg                                                            >  2:290299/1‑100 (MQ=255)                                         tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                         >  2:3291006/1‑100 (MQ=255)                                         tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                         <  2:579813/100‑1 (MQ=255)                                         tAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtg                                                         <  2:255880/100‑1 (MQ=255)                                          aaaCAGTCCTACCTTAGTTGACCCGGACTGACTGCTTGGTGTGCTGAATCAGTTAGGTCCGGTACGTGATTTTGGCCCCGACTGGCTGGACGGTAGgtgt                                                        <  1:457750/100‑1 (MQ=255)                                          aaaCAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                        <  2:3063929/100‑1 (MQ=255)                                             cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                     <  2:3096348/100‑1 (MQ=255)                                              aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaaa                                                    <  1:2758250/100‑1 (MQ=255)                                              aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaaa                                                    <  2:1324187/100‑1 (MQ=255)                                               gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                     <  1:2412257/98‑1 (MQ=255)                                               gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                   <  1:3290978/100‑1 (MQ=255)                                               gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                   <  2:1536578/100‑1 (MQ=255)                                               gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                   >  1:1983675/1‑100 (MQ=255)                                                 cATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTa                                                 >  2:3061835/1‑100 (MQ=255)                                                   tACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAg                                               >  1:240047/1‑100 (MQ=255)                                                     ccTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                             <  1:1327808/100‑1 (MQ=255)                                                     ccTTAAGTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTc                                             <  1:1478244/100‑1 (MQ=255)                                                      cTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCg                                            <  1:10823/100‑1 (MQ=255)                                                       ttAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCgg                                           <  2:2455584/100‑1 (MQ=255)                                                        tAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgt                                          >  1:1557099/1‑100 (MQ=255)                                                         aagtGACGGTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGTCTGGACGGTTGGTGTGAAACTAAg                                               <  1:791581/91‑1 (MQ=255)                                                          aTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgt                                        <  2:2749863/100‑1 (MQ=255)                                                          aTTGAAGATAACTTACCGTTTGGTGGAATGAAGCAGTCAGAACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                         <  1:32381/99‑1 (MQ=255)                                                           ttGACGCTAACTTACCGTTTGTTGTGAGGACGTGCTCAGTAACTGACCGTCATTTTGGACCTGACGGGCGGTTCGGGTGGTCTGAAAATAAGTc                                             >  1:2630910/1‑94 (MQ=255)                                                           ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                       <  1:1922400/100‑1 (MQ=255)                                                            tGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgt                                      >  2:2430205/1‑100 (MQ=255)                                                             gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtg                                     <  1:1989251/100‑1 (MQ=255)                                                              aCGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtgt                                    <  2:1557106/100‑1 (MQ=255)                                                               cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                   <  2:1063287/100‑1 (MQ=255)                                                               cGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                   >  2:940878/1‑100 (MQ=255)                                                                gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTAc                                  <  1:2639655/100‑1 (MQ=255)                                                                 cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACg                                 <  2:2810304/100‑1 (MQ=255)                                                                    aCTTCCCGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGCCGGTTGGTGTGCAACTAAATCGGTGCGAGTACGATa                              >  1:2582370/1‑100 (MQ=255)                                                                     cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                             <  2:1093691/100‑1 (MQ=255)                                                                     cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAt                             >  1:3339546/1‑100 (MQ=255)                                                                          cgctttgggggATGAAGCCGGCAGGAACGGGAAGTGAATTTGGCCTAGACTGGCTGGGCGGGTGGAGTGAAACTAAGTAGGTGTGTGTACGGTATTAAt                         >  1:1542522/9‑99 (MQ=255)                                                                          cGTTTGGTGGGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATc                        >  2:3131437/1‑100 (MQ=255)                                                                           gttttgggtgATGAAGCAGTCAGGAACGGNCCGTGTTTTTGGCCCNGACTGGCTGGACGGTTGGTCTGAAACTAANTCGGTGTGTGTACGGTNTTAATc                        <  1:3046842/90‑1 (MQ=255)                                                                             ttGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTgg                     <  2:1502175/100‑1 (MQ=255)                                                                               ggtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGCCGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                   >  1:2489866/1‑100 (MQ=255)                                                                               ggtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                   >  2:2938366/1‑100 (MQ=255)                                                                               ggtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCACCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                   >  1:279161/1‑100 (MQ=255)                                                                                gtggGATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTc                  <  2:2610432/100‑1 (MQ=255)                                                                                gtggGATGAAGCAGTCAGGACCAGGCCGAGATGTTGGCCACGACTAGCTGGACGGTGGGTGGGAAACGAAGTCGAGGTGTGAACGGTATTAATCTGGTTc                  >  1:2133561/1‑100 (MQ=255)                                                                                  gggATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGATACTAAGTCGGTTTGTGTACGGTATTAATCTGGTTCGc                >  1:848598/1‑100 (MQ=255)                                                                                   ggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt               <  1:3094235/100‑1 (MQ=255)                                                                                   ggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCt               <  1:1132911/100‑1 (MQ=255)                                                                                    gATGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTc              <  1:1046630/100‑1 (MQ=255)                                                                                     aTGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa             >  1:1693445/1‑100 (MQ=255)                                                                                       gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa           <  2:2171417/100‑1 (MQ=255)                                                                                       gAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa           <  2:784557/100‑1 (MQ=255)                                                                                        aaGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTTTACGGTATTACTCTGGTTCGCTCAt            >  1:3023490/1‑98 (MQ=255)                                                                                         aGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAg         <  1:3258402/100‑1 (MQ=255)                                                                                         aGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTGCGCTCATAAg         <  1:2384164/100‑1 (MQ=255)                                                                                             gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa     <  1:84770/100‑1 (MQ=255)                                                                                             gTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGGTGGTGTGAAACTAAGTAGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaa     >  1:712276/1‑100 (MQ=255)                                                                                               cAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa   <  2:1436127/100‑1 (MQ=255)                                                                                               cAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa   <  2:2913863/100‑1 (MQ=255)                                                                                               cAGGAACGGGCCGTGAGTTTGGCCCCGACTGGCTGGACGGTTTGTGTGAAACTAAGTCGGTGTGTGTACGGTAGTAAGCTGGTTCGCTCATAAGTaaaaa   <  2:2877193/100‑1 (MQ=255)                                                                                                 ggAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaaa  >  1:2391899/1‑99 (MQ=255)                                                                                                    |                                                                                                CTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAA  >  NC_000913/1448848‑1449042

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