Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,439,056	T→C	6.1%	S577G (AGC→GGC)	pfo ←	pyruvate‑flavodoxin oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,439,056	0	T	C	6.1%	199.0 / 3.0	66	S577G (AGC→GGC)	pfo	pyruvate‑flavodoxin oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (32/29);  new base C (2/2);  total (34/32)
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.51e-01

GCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCT  >  NC_000913/1438913‑1439204                                                                                                                                                |                                                                                                                                                     gcTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGc                                                                                                                                                 >  1:26659/1‑149 (MQ=255) gcTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGc                                                                                                                                                 <  4:88756/149‑1 (MQ=255)    gtgCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAg                                                                                                                                              <  2:9323/149‑1 (MQ=255)     tgCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAgg                                                                                                                                             >  5:144061/1‑149 (MQ=255)                  gTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTc                                                                                                                                <  2:328452/149‑1 (MQ=255)                      cAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGa                                                                                                                             <  5:115984/148‑1 (MQ=255)                       aaCGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATg                                                                                                                           >  8:345614/1‑149 (MQ=255)                        aCGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATgg                                                                                                                          <  6:144061/149‑1 (MQ=255)                            aaCTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaa                                                                                                                       >  6:282295/1‑148 (MQ=255)                             aCTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaa                                                                                                                     <  5:282295/149‑1 (MQ=255)                              cttcttCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATTGaaaaaa                                                                                                                    <  2:234505/149‑1 (MQ=255)                              cttcttCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaaa                                                                                                                    <  2:151775/149‑1 (MQ=255)                              cttcttCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGaaaaaa                                                                                                                    <  6:221865/149‑1 (MQ=255)                               ttcttcTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAg                                                                                                                   >  5:101478/1‑149 (MQ=255)                                     tACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCt                                                                                                             >  2:121673/1‑149 (MQ=255)                                          aTTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCAt                                                                                                         >  4:141773/1‑148 (MQ=255)                                           ttCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGa                                                                                                       <  3:141773/149‑1 (MQ=255)                                              tcgcgCCAGCGCCAGAGCATTCCAGTTGCGTTCCACCAGATCCTGGCCTTTGATACTTTATCTTTTGGCTATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCATGCTGTATTTGCGt                                                                                                                                 <  1:233243/119‑1 (MQ=255)                                              aCGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCATTCAGAATTTGCGt                                                                                                                                 >  2:233243/1‑120 (MQ=255)                                              aCGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACgg                                                                                                    >  8:81573/1‑149 (MQ=255)                                              aCGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACgg                                                                                                    >  7:312033/1‑149 (MQ=255)                                                    cagcgccagAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTaa                                                                                              <  8:312033/149‑1 (MQ=255)                                                            gagCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGgcc                                                                                      >  7:157769/1‑149 (MQ=255)                                                               ccTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTCTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGgccgc                                                                                    >  8:277186/1‑148 (MQ=255)                                                                 tGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAAGACGGGCCGCCAg                                                                                 >  4:318565/1‑149 (MQ=255)                                                                 tGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCCATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAACCCATCTGCATGACGGTATTAATACGGGCCGCCAg                                                                                 >  8:303692/1‑149 (MQ=255)                                                                 tGCCAGTTGCGTTCCACCAGATCCTGGCCTTAGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCCATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAg                                                                                 >  8:308970/1‑149 (MQ=255)                                                                            ttCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTtgc                                                                      >  3:310406/1‑149 (MQ=255)                                                                               caccaGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTaa                                                                                              <  8:133084/122‑1 (MQ=255)                                                                               caccaGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTaa                                                                                              >  7:133084/1‑122 (MQ=255)                                                                               caccaGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTATGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTAAGATGGAAAAAAGCCATCTGCATGACGGTATGAATACGGGCCGCTAGGCCACATTCGCggg                                                                   >  6:20652/1‑149 (MQ=255)                                                                                accaGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGc                                                                  >  2:265758/1‑149 (MQ=255)                                                                                 ccaGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGc                                                                  >  3:275638/1‑148 (MQ=255)                                                                                    gATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGa                                                                                                                             <  4:4661/86‑1 (MQ=255)                                                                                    gATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGa                                                                                                                             >  3:4661/1‑86 (MQ=255)                                                                                      tCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACg                                                                                          <  6:290172/119‑1 (MQ=255)                                                                                        cTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGa                                                                >  8:134001/1‑143 (MQ=255)                                                                                        cTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGa                                                                <  7:134001/143‑1 (MQ=255)                                                                                            ccTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTcgccgc                                                      <  6:303973/149‑1 (MQ=255)                                                                                                tGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCCATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGt                                                    <  7:308970/147‑1 (MQ=255)                                                                                                tGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCCATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGt                                                    <  7:303692/147‑1 (MQ=255)                                                                                                   tACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAAt                                                >  1:97849/1‑148 (MQ=255)                                                                                                     cTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATca                                              >  7:165656/1‑148 (MQ=255)                                                                                                     cTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATca                                              <  4:251523/148‑1 (MQ=255)                                                                                                       gTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGc                                                                                                           <  4:113953/85‑1 (MQ=255)                                                                                                       gTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGc                                                                                                           >  3:113953/1‑85 (MQ=255)                                                                                                        tAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACAt                                           <  1:121673/148‑1 (MQ=255)                                                                                                        tAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACAt                                           <  7:277186/148‑1 (MQ=255)                                                                                                         aaCTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATcaca                                            >  6:29308/1‑146 (MQ=255)                                                                                                                     tCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCtttt                             <  6:101478/149‑1 (MQ=255)                                                                                                                      cGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCttttt                            >  3:211223/1‑149 (MQ=255)                                                                                                                         acacTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGAAGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGAGCTTTTTTCt                         >  4:187057/4‑149 (MQ=255)                                                                                                                          cccTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACAt                                           <  5:54608/130‑1 (MQ=255)                                                                                                                          cccTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACAt                                           >  6:54608/1‑130 (MQ=255)                                                                                                                           ccTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTTGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTgg                       <  4:310406/149‑1 (MQ=255)                                                                                                                           ccTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATcaca                                            >  7:54338/1‑128 (MQ=255)                                                                                                                           ccTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATcaca                                            <  8:54338/128‑1 (MQ=255)                                                                                                                             tGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCtttttt                           <  8:143941/143‑1 (MQ=255)                                                                                                                             tGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCtttttt                           >  7:143941/1‑143 (MQ=255)                                                                                                                             tGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGtt                     >  5:101663/1‑149 (MQ=255)                                                                                                                               cAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGttt                    <  8:165656/148‑1 (MQ=255)                                                                                                                                aaTTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTaa                  >  5:233336/1‑149 (MQ=255)                                                                                                                                 aTTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAgcgc                                   <  2:281830/131‑1 (MQ=255)                                                                                                                                 aTTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAgcgc                                   >  1:281830/1‑131 (MQ=255)                                                                                                                                    attgcaGGGTGCGATTGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACg              <  8:116032/144‑1 (MQ=255)                                                                                                                                         aGGGCGCTATCGCCATGCAGAATTTGCGTCAGATGGAATAAAGCCATCTGGATGACGGTATTAATACTGGACGGTAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTg         <  5:20652/149‑1 (MQ=255)                                                                                                                                                tATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAActtct  >  5:287379/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     GCTGTGCGGATTTACCGGTTGCAACGGAACTTCTTCTACGGATTCACGCGCCAGCGCCAGAGCCTGCCAGTTGCGTTCCACCAGATCCTGGCCTTTGCTACTGTAACTTTTGGCAATCGCACCCTGCAATTCTGCGAGGGCGCTATCGCCAGGCAGAATTTGCGTCAGATGGAAAAAAGCCATCTGCATGACGGTATTAATACGGGCCGCCAGGCCACATTCGCGGGCGATTTTCGCCGCGTTAATCACATAGAAGCGCGCTTTTTTCTGGTTTAACACGGCCTGAACTTCT  >  NC_000913/1438913‑1439204

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