Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,432,990	G→C	6.4%	intergenic (‑48/+22)	folC ← / ← accD	bifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,432,990	0	G	C	6.4%	222.2 / 6.7	62	intergenic (‑48/+22)	folC/accD	bifunctional folylpolyglutamate synthase/ dihydrofolate synthase/acetyl‑CoA carboxylase, beta (carboxyltransferase) subunit
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (2/2);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GCTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCCGGA  >  NC_000913/2432843‑2433138                                                                                                                                                    |                                                                                                                                                     gcTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTCGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGt                                                                                                                                                     >  1:50898/1‑149 (MQ=255)   tCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTgg                                                                                                                                                   >  3:224375/1‑149 (MQ=255)     aaGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                                 <  4:323318/149‑1 (MQ=255)     aaGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgcc                                                                                                                                                 <  4:40611/149‑1 (MQ=255)       ggCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccct                                                                                                                                               >  1:35503/1‑149 (MQ=255)         ccGAGATCGATAGTTTTACTGTGCAGGTTTTCCATATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTAGGGTGAGAATAGCATAAGGGCAGGGCCACTGgctctgc                                                                                                                                             <  1:479770/149‑1 (MQ=255)           gagaTCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGgccctgccc                                                                                                                                           <  1:444142/149‑1 (MQ=255)              aTCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCtta                                                                                                                                        >  2:225991/1‑149 (MQ=255)               tCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttat                                                                                                                                       <  2:89865/149‑1 (MQ=255)                  aTAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcag                                                                                                                                    >  1:519835/1‑149 (MQ=255)                   tAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagt                                                                                                                                   <  2:498107/149‑1 (MQ=255)                      ttttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGgctctgcc                                                                                                                                            <  3:255544/137‑1 (MQ=255)                      ttttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGGGCCACTGgctctgcc                                                                                                                                            >  4:255544/1‑137 (MQ=255)                      ttttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttat                                                                                                                                >  2:429853/1‑149 (MQ=255)                        ttACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatcagttatca                                                                                                                              <  2:101369/149‑1 (MQ=255)                            tgtgCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGcc                                                                                                                          >  3:276006/1‑149 (MQ=255)                                  ggTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGt                                                                                                                    <  1:115140/149‑1 (MQ=255)                                  ggTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGt                                                                                                                    <  1:498565/149‑1 (MQ=255)                                    ttttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTc                                                                                                                  >  3:340911/1‑149 (MQ=255)                                    ttttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTc                                                                                                                  <  2:39471/149‑1 (MQ=255)                                     tttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTcc                                                                                                                 <  3:73696/149‑1 (MQ=255)                                     tttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTcc                                                                                                                 <  3:82010/149‑1 (MQ=255)                                      ttCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCt                                                                                                                <  4:527175/149‑1 (MQ=255)                                        ccAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGa                                                                                                              <  3:268870/149‑1 (MQ=255)                                         cAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGAt                                                                                                             >  1:353266/1‑149 (MQ=255)                                           gATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatca                                                                                                                                     <  1:368343/123‑1 (MQ=255)                                           gATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCttatca                                                                                                                                     >  2:368343/1‑123 (MQ=255)                                           gATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATcc                                                                                                           <  3:107676/149‑1 (MQ=255)                                                 aaaGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTAcc                                                                                                     >  2:79866/1‑149 (MQ=255)                                                   aGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgg                                                                                                   >  2:541846/1‑149 (MQ=255)                                                     ccACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCgggg                                                                                                 >  3:336576/1‑149 (MQ=255)                                                              cAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTa                                                                                        >  1:514053/1‑149 (MQ=255)                                                              cAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTa                                                                                        >  1:51683/1‑149 (MQ=255)                                                                     gACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTc                                                                                 <  4:486256/149‑1 (MQ=255)                                                                        gcgGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCAcg                                                                              >  4:550066/1‑149 (MQ=255)                                                                           gCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACgcgg                                                                           <  2:520964/149‑1 (MQ=255)                                                                            cTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCt                                                                                   <  1:2953/140‑1 (MQ=255)                                                                            cTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCt                                                                                   >  2:2953/1‑140 (MQ=255)                                                                               gaggagTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCt                                                                       >  2:134965/1‑149 (MQ=255)                                                                                  gagTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCa                                                                    >  3:127087/1‑149 (MQ=255)                                                                                            gATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGgcg                                                          <  2:239518/149‑1 (MQ=255)                                                                                                    tGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGa                                                  <  1:225991/149‑1 (MQ=255)                                                                                                            gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCAcgcg                                                                            >  4:530237/1‑115 (MQ=255)                                                                                                            gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCAcgcg                                                                            <  3:530237/115‑1 (MQ=255)                                                                                                            gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCaa                                          <  2:56038/149‑1 (MQ=255)                                                                                                            gCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCaa                                          <  1:370426/149‑1 (MQ=255)                                                                                                             cTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAAc                                         >  4:328655/1‑149 (MQ=255)                                                                                                                aTTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTc                                      >  1:391209/1‑149 (MQ=255)                                                                                                                            gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTc                          <  1:84474/149‑1 (MQ=255)                                                                                                                            gagaATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTc                          <  4:118207/149‑1 (MQ=255)                                                                                                                               aaTAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGcc                       <  2:136965/149‑1 (MQ=255)                                                                                                                                     aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTc                 <  4:126489/149‑1 (MQ=255)                                                                                                                                     aaaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTc                 >  2:557489/1‑149 (MQ=255)                                                                                                                                      aaaGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCa                >  4:491206/1‑149 (MQ=255)                                                                                                                                        aGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAgg              <  2:353266/149‑1 (MQ=255)                                                                                                                                         gggCAGGGCCACTGGCTCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgc             >  2:138975/1‑149 (MQ=255)                                                                                                                                         gggCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgc             <  3:14883/149‑1 (MQ=255)                                                                                                                                         gggCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgc             <  2:519835/149‑1 (MQ=255)                                                                                                                                          ggCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgcg            >  2:344035/1‑149 (MQ=255)                                                                                                                                           gCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgcgc           >  3:430028/1‑149 (MQ=255)                                                                                                                                           gCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGgcgc           <  1:368228/149‑1 (MQ=255)                                                                                                                                                gCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTc      >  4:300876/1‑149 (MQ=255)                                                                                                                                                    gTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCcgga  <  2:514053/149‑1 (MQ=255)                                                                                                                                                    gTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCcgga  <  2:51683/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GCTCAAGGCCGAGATCGATAGTTTTACTGTGCAGGTTTTCCAGATAAGAAAGCCACGAAGCCAGAGGCGACGCGGCTTGAGGAGTGCGTTTGATAATCATGGTATCCGCTGATTCGTTACGGTGAGAATAGCAAAAGGGCAGAGCCAGTGGCCCTGCCCTTATCAGTTATCAGGCCTCAGGTTCCTGATCCGGTACCGGGGGTACCACTACGCCTTCACGCGGCGCTTCAGGATTCGGCGCTGGCAGATTCATCAACTTCGCCAGAATGCTCGCCAGTTTCAGGCGCATTTCCGGA  >  NC_000913/2432843‑2433138

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