Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,246,073	C→T	100%	G435D (GGC→GAC)	treA ←	periplasmic trehalase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,246,073	0	C	T	96.1%	228.3 / ‑0.1	76	G435D (GGC→GAC)	treA	periplasmic trehalase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (31/42);  total (32/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

CCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCC  >  NC_000913/1245974‑1246168                                                                                                    |                                                                                                ccGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGgc                                                                                                 >  2:170817/1‑100 (MQ=255)      gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  2:130237/1‑100 (MQ=255)      gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  2:2172360/1‑100 (MQ=255)      gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  1:1508959/1‑100 (MQ=255)      gTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGtcgccg                                                                                            >  2:315535/1‑100 (MQ=255)        tttGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCggg                                                                                          <  1:2573697/100‑1 (MQ=255)         ttGTAATCCTTCTGTCGCGACCCACTGTATCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGt                                                                                         <  1:144405/100‑1 (MQ=255)            tAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgc                                                                                      <  1:154105/100‑1 (MQ=255)            tAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgc                                                                                      >  2:324582/1‑100 (MQ=255)             aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTgca                                                                                     <  2:1478617/100‑1 (MQ=255)             aaTCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCAATTGTTGCCCACTTTTCACCGCCGTGGTGTTCAGCTCCCCGGGtt                                                                                        >  1:1125779/1‑97 (MQ=255)               tCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagc                                                                                   <  2:94274/100‑1 (MQ=255)                ccTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTgcagca                                                                                  <  1:2379305/100‑1 (MQ=255)                 cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                 >  2:207314/1‑100 (MQ=255)                 cTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGcagcag                                                                                 >  2:1999897/1‑100 (MQ=255)                    cTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGAtg                                                                              <  1:2058918/100‑1 (MQ=255)                      gTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGAtgtg                                                                            >  1:1090356/1‑100 (MQ=255)                          cgACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGtttt                                                                        <  2:1428599/100‑1 (MQ=255)                           gACCCACGTTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCTGGGTTGCAGCAGATGTGtttt                                                                        <  1:1264816/99‑1 (MQ=255)                            aCCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcg                                                                      >  1:137373/1‑100 (MQ=255)                              tcACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                                    <  1:970470/99‑1 (MQ=255)                              ccACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtc                                                                    <  2:796361/100‑1 (MQ=255)                               cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGCCGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtcg                                                                   >  2:85872/1‑100 (MQ=255)                               cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTtcgtca                                                                   >  2:970472/1‑99 (MQ=255)                               cACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGACCACTTTTCACCAACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTGtcgtcg                                                                   >  2:582024/1‑100 (MQ=255)                                 cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgcc                                                                 <  1:2519774/100‑1 (MQ=255)                                 cTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgcc                                                                 <  1:1320136/100‑1 (MQ=255)                                     aaCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcg                                                             <  2:2494301/100‑1 (MQ=255)                                      aCGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTcgtcgccgtcgc                                                            >  1:1765560/1‑100 (MQ=255)                                           gCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGATGCAGCAGATGTGTTTTCGTCGCCGTCGCCActt                                                       >  2:2642478/1‑100 (MQ=255)                                            cccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCAtttt                                                      >  2:1802141/1‑100 (MQ=255)                                             ccAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttg                                                     <  2:2712195/100‑1 (MQ=255)                                              cAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgt                                                    >  1:499215/1‑100 (MQ=255)                                               acccATTTGCCGCATCCCATTCTTGCCCACTTTTCACCGACGTGGTGTTCCGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgtt                                                   <  1:2082194/98‑1 (MQ=255)                                                gCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTttgttg                                                  >  1:218218/1‑100 (MQ=255)                                                   aTTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGgcg                                               >  1:250582/1‑100 (MQ=255)                                                       ggCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                           <  1:2505334/100‑1 (MQ=255)                                                       ggCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                           >  2:1231574/1‑100 (MQ=255)                                                       ggCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGATGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGAt                                           <  1:1608586/100‑1 (MQ=255)                                                           cATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCttt                                       <  2:2678179/100‑1 (MQ=255)                                                           cATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCttt                                       <  1:2013254/100‑1 (MQ=255)                                                            aTCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTc                                      <  2:2901207/100‑1 (MQ=255)                                                              cccATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCg                                     <  2:1117066/99‑1 (MQ=255)                                                                cATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTg                                  >  1:872507/1‑100 (MQ=255)                                                                cATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTg                                  <  2:1950751/100‑1 (MQ=255)                                                                  ttgttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGcc                                <  1:1999890/100‑1 (MQ=255)                                                                   tgttgCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCACTGCCg                               <  1:2638777/100‑1 (MQ=255)                                                                       gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCAtt                           <  1:506613/100‑1 (MQ=255)                                                                       gCCCACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTAGCCATTTTGTTGGCGCGTTCTTTCGCTGCTGTAtt                           >  1:931594/1‑100 (MQ=255)                                                                        cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          <  1:1437321/100‑1 (MQ=255)                                                                        cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          <  2:1508963/100‑1 (MQ=255)                                                                        cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          >  2:1231752/1‑100 (MQ=255)                                                                        cccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTg                          <  2:321463/100‑1 (MQ=255)                                                                         ccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGa                         >  1:2330251/1‑100 (MQ=255)                                                                         ccACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGa                         >  2:1705063/1‑100 (MQ=255)                                                                          cACTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGAc                        >  1:2012960/1‑100 (MQ=255)                                                                           aCTTTTCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACg                       <  2:1449345/100‑1 (MQ=255)                                                                               ttCACCGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACa                   <  1:1845512/100‑1 (MQ=255)                                                                                 cACCGACGTGGAGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGc                 >  1:1059602/1‑100 (MQ=255)                                                                                   ccGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg               <  2:2067878/100‑1 (MQ=255)                                                                                   ccGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg               >  2:2601297/1‑100 (MQ=255)                                                                                   ccGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg               >  2:1492869/1‑100 (MQ=255)                                                                                    cGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCggg              <  1:921004/100‑1 (MQ=255)                                                                                    cGACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCgg               <  1:753328/99‑1 (MQ=255)                                                                                     gACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGa             <  1:317949/100‑1 (MQ=255)                                                                                     gACGTGGTGTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGa             <  2:45897/100‑1 (MQ=255)                                                                                      aCGTGGTGGTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGaa            <  1:2756335/100‑1 (MQ=255)                                                                                        gtggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACa          <  2:2504342/100‑1 (MQ=255)                                                                                         tggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTGGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAg         >  1:953147/1‑100 (MQ=255)                                                                                         tggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAg         >  1:2070859/1‑100 (MQ=255)                                                                                          ggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAgg        <  2:2148031/100‑1 (MQ=255)                                                                                          ggtgTTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAg         <  1:118069/99‑1 (MQ=255)                                                                                           gtgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAggg       <  2:1608971/100‑1 (MQ=255)                                                                                            tgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggc      <  2:2079888/100‑1 (MQ=255)                                                                                            tgtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggc      >  2:2644698/1‑100 (MQ=255)                                                                                             gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     >  2:54578/1‑100 (MQ=255)                                                                                             gtTCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcg     >  1:884498/1‑100 (MQ=255)                                                                                              ttCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcgg    <  2:1658806/100‑1 (MQ=255)                                                                                               tCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcggc   >  1:542423/1‑100 (MQ=255)                                                                                               tCAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGggcggc   <  1:2270644/100‑1 (MQ=255)                                                                                                cAGGTCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGcc  >  1:1605577/1‑100 (MQ=255)                                                                                                    |                                                                                                CCGTAGTTTTGTAATCCTTCTGTCGCGACCCACTGTAACGGTGCCCAGCCATTTGGCGCATCCCATTGTTGCCCACTTTTCACCGACGTGGTGTTCAGGCCGCCGGGTTGCAGCAGATGTGTTTTCGTCGCCGTCGCCATTTTGTTGGCGCGATCTTTCGCTGCCGCATTGACGTACAGCGGGAACAGGGCGGCC  >  NC_000913/1245974‑1246168

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