Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	954,156	C→T	6.6%	R104H (CGC→CAC)	focA ←	formate channel

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	954,156	0	C	T	6.6%	189.5 / 2.0	60	R104H (CGC→CAC)	focA	formate channel
Reads supporting (aligned to +/- strand):  ref base C (27/29);  new base T (2/2);  total (29/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGTTTGTAGGACGTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTG  >  NC_000913/954010‑954303                                                                                                                                                   |                                                                                                                                                    cggTTTGTAGGACGTTTAGTCCCCATGGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGTCCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTTATgcgc                                                                                                                                                   <  4:229747/149‑1 (MQ=255)       gTAGGACGTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTc                                                                                                                                             >  8:69998/1‑149 (MQ=255)         aGGACGTTTAGTCCCCATTTACCATTTGAGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGTGCCCACTCGc                                                                                                                                           >  8:30943/1‑149 (MQ=255)         aGGACGTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGTGCCCACTCGc                                                                                                                                           >  8:25754/1‑149 (MQ=255)          ggACGTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGcc                                                                                                                                          >  4:163086/1‑149 (MQ=255)             agTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGGGATGCGCCCACTCGCCTTa                                                                                                                                       <  4:246227/148‑1 (MQ=255)               tttAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAg                                                                                                                                      >  7:331779/1‑148 (MQ=255)               tttAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGc                                                                                                                                     <  3:163086/149‑1 (MQ=255)                    tCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGcaacaa                                                                                                                                >  8:125655/1‑149 (MQ=255)                    tCCCCACTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGAGCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCTCCTTAGcaacaa                                                                                                                                <  1:349488/149‑1 (MQ=255)                      cccATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGcaacaac                                                                                                                               <  1:21655/148‑1 (MQ=255)                       ccATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGcaacaacaa                                                                                                                             >  7:226230/1‑149 (MQ=255)                        cATTGACCATTTGCGGTCATATATTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGAGGCGCCCACTCGCCTTAGCAACAACAAt                                                                                                                            <  2:102682/149‑1 (MQ=255)                          ttGACCATTTGCGGTCATATACTCGACGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCACAGGTGATGCGCCCACTCGCCTTAGAAACAACAATc                                                                                                                           >  4:94217/1‑148 (MQ=255)                             aCCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAaca                                                                                                                       <  6:22514/149‑1 (MQ=255)                               cATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACg                                                                                                                     >  3:312938/1‑149 (MQ=255)                                  ttGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTg                                                                                                                  >  3:331082/1‑149 (MQ=255)                                    gCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTgg                                                                                                                 <  8:331779/148‑1 (MQ=255)                                     cGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACTAACATCAGTGCGCCGAACATGTTGCCAATATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGaa                                                                                                               <  3:76864/149‑1 (MQ=255)                                     cGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGaa                                                                                                               <  8:227667/149‑1 (MQ=255)                                        tCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGt                                                                                                             >  6:277921/1‑148 (MQ=255)                                          atatACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTaa                                                                                                           <  5:325005/148‑1 (MQ=255)                                              aCTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAagag                                                                                                       <  2:16272/148‑1 (MQ=255)                                              aCTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAagag                                                                                                       >  1:16272/1‑148 (MQ=255)                                               cTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGAt                                                                                                     <  7:64232/149‑1 (MQ=255)                                                tCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTaaa                                                                                                          >  2:66454/1‑143 (MQ=255)                                                tCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACAAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACAGTGGAAGTaaa                                                                                                          >  3:123076/1‑143 (MQ=255)                                                  gCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATc                                                                                                    <  5:322120/147‑1 (MQ=255)                                                       aaaGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCg                                                                                             >  1:20645/1‑149 (MQ=255)                                                          gCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGTGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAg                                                                                          <  7:25754/149‑1 (MQ=255)                                                          gCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGTGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGAGCGGCTCCGCAg                                                                                          <  7:30943/149‑1 (MQ=255)                                                              caTTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGacaa                                                                                      <  4:331082/149‑1 (MQ=255)                                                                  aaaaGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAacac                                                                                   >  8:73389/1‑148 (MQ=255)                                                                        aCAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGaa                                                                             <  3:177438/148‑1 (MQ=255)                                                                            aCAGAAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCATTTTTTCGCCAACTGACCACATGTGATGCGCCCACTCGCATTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGc                                                                        >  2:43728/1‑149 (MQ=255)                                                                                  gTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCaga                                                                  <  8:226230/149‑1 (MQ=255)                                                                                  gTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCa                                                                    >  5:148173/1‑147 (MQ=255)                                                                                     cgcCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCagagag                                                               <  4:312938/149‑1 (MQ=255)                                                                                     cgcCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCACACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCagagag                                                               <  8:57510/149‑1 (MQ=255)                                                                                      gcCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAgagaga                                                              <  5:266786/149‑1 (MQ=255)                                                                                            cAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGTTCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCaaa                                                        <  7:69998/149‑1 (MQ=255)                                                                                                  gCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATgccg                                                   <  7:73389/148‑1 (MQ=255)                                                                                                  gCCAAAATAGACATTTAGCCAGTTTTTCGACAAATGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATgccgc                                                  >  2:289828/1‑149 (MQ=255)                                                                                                    cAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATgccgcc                                                 >  5:111004/1‑148 (MQ=255)                                                                                                    cAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATgccgcc                                                 >  7:134294/1‑148 (MQ=255)                                                                                                     aaaaTATACATATAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAAAAACAATCAACACGGTGGAAGTAAAGAGATCGGCGCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCaa                                               <  1:184631/149‑1 (MQ=255)                                                                                                     aaaaTAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCATTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCaa                                               <  2:192472/149‑1 (MQ=255)                                                                                                      aaaTAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAAc                                              >  8:211044/1‑149 (MQ=255)                                                                                                              aTTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCagagag                                                               >  5:106134/1‑124 (MQ=255)                                                                                                              aTTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCagagag                                                               <  6:106134/124‑1 (MQ=255)                                                                                                              aTTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTg                                      >  1:50951/1‑149 (MQ=255)                                                                                                                 tAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGcc                                    <  5:86869/148‑1 (MQ=255)                                                                                                                   gCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATg                                 >  5:39830/1‑149 (MQ=255)                                                                                                                    ccAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGc                                >  5:110398/1‑149 (MQ=255)                                                                                                                       gTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGCGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATTCCGCCAACCAGTTTTGCCATGCCg                              >  4:196053/1‑148 (MQ=255)                                                                                                                              gCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCAt                      <  1:66454/149‑1 (MQ=255)                                                                                                                                      aCCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATtgtgcct               <  1:323965/148‑1 (MQ=255)                                                                                                                                        cccAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATtgtgcctgtg            >  5:339118/1‑149 (MQ=255)                                                                                                                                           aGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATtgtgcctgtgc           >  2:292022/1‑147 (MQ=255)                                                                                                                                                aTGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTgg     <  7:211044/148‑1 (MQ=255)                                                                                                                                                aTGCGCCCACTCGCCTTAGCAACAACAATCAACACAGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCTCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTgg     <  4:123076/148‑1 (MQ=255)                                                                                                                                                  gcgcCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTg  >  7:139382/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGGTTTGTAGGACGTTTAGTCCCCATTGACCATTTGCGGTCATATACTCGCCGGAAAGCCACATTAAAAGTACAAACAGCAGTGCGCCGACCAGGTTGCCAAAATAGACATTTAGCCAGTTTTTCGCCAACTGACCCCAGGTGATGCGCCCACTCGCCTTAGCAACAACAATCAACACGGTGGAAGTAAAGAGATCGGCTCCGCAGACAACACAAAGAATCAGCCCCAGAGAGAAGCAAATGCCGCCAACCAGTTTTGCCATGCCGAAGGGCATTGTGCCTGTGCCAGTGGTTG  >  NC_000913/954010‑954303

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