Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,007,741	T→A	D117V (GAT→GTT)	yedD ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,007,741	0	T	A	98.9%	307.5 / ‑4.0	87	D117V (GAT→GTT)	yedD	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base A (37/49);  total (37/50)
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.68e-01

ATTCATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGATCCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGATCA  >  NC_000913/2007644‑2007836                                                                                                  |                                                                                                aTTCATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAAcc                                                                                               <  1:152639/100‑1 (MQ=255)  ttCATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACca                                                                                              >  1:379816/1‑100 (MQ=255)   tCATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACcac                                                                                             <  2:75481/100‑1 (MQ=255)    cATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACcaca                                                                                            >  2:2259636/1‑100 (MQ=255)     aTCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCacac                                                                                           <  1:1841461/100‑1 (MQ=255)      tCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACg                                                                                          <  2:1358463/100‑1 (MQ=255)       ccccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGc                                                                                         >  1:611186/1‑100 (MQ=255)       ccccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGc                                                                                         <  1:2337853/100‑1 (MQ=255)        cccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCt                                                                                        >  1:387829/1‑100 (MQ=255)        cccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCt                                                                                        >  1:866433/1‑100 (MQ=255)         ccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGATCCACACGCTc                                                                                       <  1:49430/100‑1 (MQ=255)         ccGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTc                                                                                       >  1:2307570/1‑100 (MQ=255)              ggcAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCg                                                                                  >  1:25788/1‑100 (MQ=255)               gcAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCTCACGCTCCATCGt                                                                                 <  1:1572648/100‑1 (MQ=255)                 aaGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCa                                                                               <  1:161349/100‑1 (MQ=255)                   gCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATa                                                                             <  1:1313941/100‑1 (MQ=255)                   gCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATa                                                                             <  1:113576/100‑1 (MQ=255)                     cGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATAcc                                                                           <  1:2080377/100‑1 (MQ=255)                        ggAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATc                                                                        >  1:1995264/1‑100 (MQ=255)                          agaTTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGt                                                                      >  1:2125820/1‑100 (MQ=255)                          agaTTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGt                                                                      <  2:460024/100‑1 (MQ=255)                            aTTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTAc                                                                    >  2:1876216/1‑100 (MQ=255)                            aTTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTAc                                                                    <  1:1757886/100‑1 (MQ=255)                             tttCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACt                                                                   >  2:2059575/1‑100 (MQ=255)                              ttCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTc                                                                  >  2:368665/1‑100 (MQ=255)                              ttCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTc                                                                  >  2:375476/1‑100 (MQ=255)                                   tACTGCAGTGGCGTCCGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATgg                                                             <  2:2241197/100‑1 (MQ=255)                                    aCGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGgt                                                            >  2:596749/1‑100 (MQ=255)                                     cGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGgtg                                                           <  1:1445018/100‑1 (MQ=255)                                      ggCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGgtgt                                                          <  1:1211156/100‑1 (MQ=255)                                      ggCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGgtgt                                                          >  1:1330541/1‑100 (MQ=255)                                       gCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGtt                                                         <  2:1995257/100‑1 (MQ=255)                                       gCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGtt                                                         <  2:137373/100‑1 (MQ=255)                                          gTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTAcc                                                      >  2:840095/1‑100 (MQ=255)                                            ggCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGt                                                    >  2:892202/1‑100 (MQ=255)                                              cGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCa                                                  <  1:732606/100‑1 (MQ=255)                                                   gTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACATGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTc                                             <  2:1299461/100‑1 (MQ=255)                                                   gTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTc                                             <  2:2125813/100‑1 (MQ=255)                                                    tAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCa                                            <  1:596750/100‑1 (MQ=255)                                                      aCGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCCTCGTACTCAATGGTGTTACCGTCACGCTCAAc                                          >  1:1021471/1‑100 (MQ=255)                                                       cGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAAcc                                         >  2:2278069/1‑100 (MQ=255)                                                       cGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAAcc                                         <  1:892204/100‑1 (MQ=255)                                                        ggCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCt                                        >  2:176481/1‑100 (MQ=255)                                                        ggCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCt                                        >  1:1000858/1‑100 (MQ=255)                                                          cgcTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCg                                      <  2:262533/100‑1 (MQ=255)                                                           gcTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGt                                     <  2:129775/100‑1 (MQ=255)                                                           gcTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGt                                     <  1:2043953/100‑1 (MQ=255)                                                             tttATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTAc                                   >  2:1218208/1‑100 (MQ=255)                                                             tttATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTAc                                   <  1:2259644/100‑1 (MQ=255)                                                             tttATCCAGCGCGGCGGCACACTCGGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTAc                                   <  1:5465/100‑1 (MQ=255)                                                              ttttCCAGCCCGGCGGCACACTCTGCCGTCGGGCGGACCACCCCCCCCCTCGTCATACCATCCTACCCAAATGTGTTACCCTCCCGCCCCACCCCGTAc                                   >  2:1752464/4‑99 (MQ=255)                                                              ttATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACa                                  >  2:2027746/1‑100 (MQ=255)                                                                aTCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGc                                <  1:29537/100‑1 (MQ=255)                                                                 tCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGc                                >  1:1413409/1‑99 (MQ=255)                                                                 tCCAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCt                               <  1:181780/100‑1 (MQ=255)                                                                  ccAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCt                               >  2:1434680/1‑99 (MQ=255)                                                                  ccAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTc                              >  1:2066633/1‑100 (MQ=255)                                                                   cAGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCa                             >  2:2338996/1‑100 (MQ=255)                                                                    aGCGCGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCAc                            >  2:248787/1‑100 (MQ=255)                                                                     gcgcGGCGGCACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACg                           >  1:2078291/1‑100 (MQ=255)                                                                         ggcggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGtttt                       <  2:138926/100‑1 (MQ=255)                                                                          gcggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTa                      >  1:33332/1‑100 (MQ=255)                                                                          gcggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTa                      <  2:379814/100‑1 (MQ=255)                                                                          gcggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTa                      <  2:1330538/100‑1 (MQ=255)                                                                           cggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTAc                     <  1:1993335/100‑1 (MQ=255)                                                                           cggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTAc                     <  1:2339005/100‑1 (MQ=255)                                                                           cggcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTAc                     <  2:2307561/100‑1 (MQ=255)                                                                             gcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACgg                   <  2:593226/100‑1 (MQ=255)                                                                             gcACACTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACgg                   <  1:840097/100‑1 (MQ=255)                                                                                cacTCTGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTaa                >  1:1159110/1‑100 (MQ=255)                                                                                  ctctGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACg              <  2:1785504/100‑1 (MQ=255)                                                                                  ctctGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACg              <  2:70922/100‑1 (MQ=255)                                                                                  ctctGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACg              <  2:1488163/100‑1 (MQ=255)                                                                                  ctctGCTGTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACg              <  1:521655/100‑1 (MQ=255)                                                                                    ctgctgTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGtt            >  2:1051967/1‑100 (MQ=255)                                                                                      gctgTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACGTCGTACAGCTCACGTTTTACGGTAACGTTAg          >  2:1570733/1‑100 (MQ=255)                                                                                        tgTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTg        <  1:375478/100‑1 (MQ=255)                                                                                        tgTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTg        <  2:1409392/100‑1 (MQ=255)                                                                                        tgTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTg        <  2:2087770/100‑1 (MQ=255)                                                                                        tgTCGCGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTg        <  1:368667/100‑1 (MQ=255)                                                                                         gTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTga       >  1:1862246/1‑100 (MQ=255)                                                                                         gTCGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTga       <  2:1001389/100‑1 (MQ=255)                                                                                           cGGGCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTgaga     >  1:1866067/1‑100 (MQ=255)                                                                                            gggCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGat    <  2:2066626/100‑1 (MQ=255)                                                                                             ggCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGatc   <  2:1159108/100‑1 (MQ=255)                                                                                              gCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGatca  >  1:1744534/1‑100 (MQ=255)                                                                                              gCGAACCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGatca  <  2:1489407/100‑1 (MQ=255)                                                                                                  |                                                                                                ATTCATCCCCGGCGGCAAGCCGGGGAGATTTCATTACGGCAGTGGCGTCGGTAACGGCGCTTTATCCAGCGCGGCGGCACACTCTGCTGTCGGGCGATCCACACGCTCCATCGTCATACCATCGTACTCAATGGTGTTACCGTCACGCTCAACCTCGTACAGCTCACGTTTTACGGTAACGTTAGTGAGATCA  >  NC_000913/2007644‑2007836

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