Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	646,747	0	T	G	40.0%	109.1 / 12.6	61	T255P (ACC→CCC)	citG	2‑(5''‑triphosphoribosyl)‑3'‑dephosphocoenzyme‑A synthase
Reads supporting (aligned to +/- strand):  ref base T (8/28);  new base G (24/0);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 3.80e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAA  >  NC_000913/646600‑646894                                                                                                                                                     |                                                                                                                                                       atatTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑g                                                                                                                                                    <  1:166792/151‑1 (MQ=255)     ttCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAAt                                                                                                                                                 <  1:958023/150‑1 (MQ=255)      tCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATg                                                                                                                                                <  2:605291/150‑1 (MQ=255)       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑c                                                                                                                                              <  1:282341/150‑1 (MQ=255)       ccAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGAGGTTT‑Gaa                                                                                                                                                  >  2:267015/1‑147 (MQ=255)         aaaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑cccc                                                                                                                                           <  1:664637/151‑1 (MQ=255)          aaTTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑ccccc                                                                                                                                          <  1:426892/151‑1 (MQ=255)               aTCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAAG‑CCCCCctttt                                                                                                                                     >  1:415326/1‑151 (MQ=255)                  aaGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCctttt                                                                                                                                     >  2:981580/1‑148 (MQ=255)                   aGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCttttt                                                                                                                                     >  2:270791/1‑147 (MQ=255)                    gTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTTCaa                                                                                                                                 >  2:358433/1‑150 (MQ=255)                     tGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATG‑CCCCCctttt                                                                                                                                     >  2:923506/1‑145 (MQ=255)                      gCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCaata                                                                                                                               <  1:497081/150‑1 (MQ=255)                        ttAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAATGACCCCCTTTTTG‑aataa                                                                                                                              >  1:393328/1‑149 (MQ=255)                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCaataat                                                                                                                             <  1:7194/149‑1 (MQ=255)                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATg                                                                                                                            <  2:575365/150‑1 (MQ=255)                         tAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑CAATG‑CCCCCctttttttataa                                                                                                                              >  1:119636/1‑141 (MQ=255)                          aaaTAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGtt                                                                                                                          <  2:315576/151‑1 (MQ=255)                            ataatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGGCGGAAGG‑CCCCCTTTTTG‑aataaat                                                                                                                            >  1:35207/1‑145 (MQ=255)                             taatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATG‑CCCCCCTTTTTCcataaatttt                                                                                                                         >  2:241337/1‑144 (MQ=255)                              aatTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTgc                                                                                                                       <  2:410951/150‑1 (MQ=255)                                ttAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTGTAATAATAttttggcc                                                                                                                    >  1:611121/1‑146 (MQ=255)                                  aaaTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAATG‑CCCCCCTTTTTCAATAAAGTTTGCGcccc                                                                                                                  >  1:489396/1‑149 (MQ=255)                                    aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTTCCATTATGttttccgctccc                                                                                                                >  2:895428/1‑142 (MQ=255)                                    aTCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAATAtttttgcc                                                                                                                    >  2:290612/1‑142 (MQ=255)                                         tgCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GCATG‑CCCCCCTTTTTCAATAAAGTTTTTGCCTCCCGCCGt                                                                                                           >  1:1009111/1‑151 (MQ=255)                                            tAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATT‑CCCCCCTTTTTCTATAAAGTTTGCGCCCcgc                                                                                                                >  1:468553/1‑143 (MQ=255)                                             aaaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGc                                                                                                        <  2:494623/150‑1 (MQ=255)                                              aaaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCa                                                                                                      <  2:741512/151‑1 (MQ=255)                                               aaaCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAg                                                                                                     <  2:256330/151‑1 (MQ=255)                                                 aCCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAgcg                                                                                                   <  1:852923/151‑1 (MQ=255)                                                  ccAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GGATG‑TCCCCCTTTTTCAAAAATGTTTTCGGCCCGCGCTGTAgac                                                                                                       >  2:646457/1‑144 (MQ=255)                                                   cAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCa                                                                                                 <  1:62262/151‑1 (MQ=255)                                                          ggATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTTCAATAATGTTTTCGCCCCCCGCCGTAGCCCGCGCCGGccccc                                                                                          >  1:1100198/1‑151 (MQ=255)                                                                gTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCAGAGATAATCGAGATCGGC‑GGGGGGTC‑GGAAG‑CCCCCCTTTTTTAATAAAGTTTTCGCCTCGCGGCGTTGCCAGCGCCGGCCCCCCCCGcc                                                                                    >  2:21217/1‑151 (MQ=255)                                                                 tAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTT‑GAAAG‑GCCCCCTTTTGCAATAAATTTTTCGCCTCCCGCCGTTGGCAGCGCCGGcccccccccccc                                                                                   >  1:1007977/1‑145 (MQ=255)                                                                   ggTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAATGTTTTGCCCCCGCGCTGTAGCCAGCGCAGGgccccccccccccc                                                                                 >  2:233324/1‑143 (MQ=255)                                                                              ccgccTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑TAATG‑GCCCCCTTTTTCAAAAAAATTTTTGGCCCGCGCCGGAGCCCCCGCAGGCCCCCCCCCCCACCCCCCGGTGCa                                                                       >  1:46797/1‑150 (MQ=255)                                                                                   tGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGCGGCTA‑CTGTT‑ACCCCTATTTTTCATAATGGTTAGGCCTCCCCCGCTAGCCAGCGGTTGGCCCCCTCCCCTCcccgcgccgggacag                                                                   >  2:214298/1‑140 (MQ=255)                                                                                         gagaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATc                                                           <  2:598881/150‑1 (MQ=255)                                                                                           gaTTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGcc                                                        <  2:379817/151‑1 (MQ=255)                                                                                                        acacTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GAATG‑CCCCCCTTTTTCAATAAAGTTTTCGCCCCGCGCCGTAGCCAGCGCAGGCCCCCCCCCCCCCCCCGCGAAGCAACGTT‑GGGATCGCCGTTTATCGcc                                              >  2:927300/1‑148 (MQ=255)                                                                                                                tGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGCGGGTC‑GAAAT‑TCCCCTTTTTTCTATAAAGTTTTGGGCCCGCCCTGTTGGCCACGCAGGCCCCCCCCGCCCCCCCGCGCGGCAACGCT‑TGTATCGCCCTTTGTTGCCCTCCGTAg                                      >  2:824615/1‑148 (MQ=255)                                                                                                                   cGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCaa                                <  2:959948/151‑1 (MQ=255)                                                                                                                     aaCTGCCGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTTCAATAAAGTTTTCGCCCCGCCCTGGAGCCCGCGCAGGCCCCCCCCCCCCCCCCGCCGTGGCACCGT‑GGGATCGCCGTTGATCGcccaccgta                                       >  2:360080/1‑135 (MQ=255)                                                                                                                         gCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTAt                           <  2:475480/150‑1 (MQ=255)                                                                                                                          ccGGAGATAATCGAGATCGGCGGGGGGTTG‑AAAGG‑CCCCCTTTTTGCAA‑AAAGATTGCTCCCCGCCGTGTAACCCGCGCCGGGCCCCCCCCCCCCCCCCCGAAGGAACGTTGGGTTCCGCCG‑TGCTCGCCATCAGGCAGAGCaaagga                            >  1:373862/1‑144 (MQ=255)                                                                                                                           cGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGa                        <  1:746096/151‑1 (MQ=255)                                                                                                                           cGGAGATAATCGAGATCGGC‑GGGGGGTC‑GGATG‑CCCCCCTTTTGTAATAATGTTTTCGCCCCGCGCCGTAGCCACCGCAGGCCCCCCCCCCCCCCCCCCGACGCCACGCT‑GGGGTCGCCGCTGATCGCCCCCCGTAGGGGCAAggga                            >  1:857013/1‑145 (MQ=255)                                                                                                                             gagaTAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGc                      <  1:200596/151‑1 (MQ=255)                                                                                                                                       agaTCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt            <  1:214550/151‑1 (MQ=255)                                                                                                                                        gaTCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt            <  2:814758/150‑1 (MQ=255)                                                                                                                                          tCGGC‑GGGGGGTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGtt            >  1:962341/1‑148 (MQ=255)                                                                                                                                           cggc‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAgg        <  1:165541/151‑1 (MQ=255)                                                                                                                                             gc‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt      <  2:664945/151‑1 (MQ=255)                                                                                                                                             gc‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGAt      <  1:46316/151‑1 (MQ=255)                                                                                                                                              c‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc     >  2:682/1‑151 (MQ=255)                                                                                                                                              c‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc     >  2:419333/1‑151 (MQ=255)                                                                                                                                              c‑ggGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATc     >  2:117490/1‑151 (MQ=255)                                                                                                                                                gggggTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCt    >  2:461225/1‑151 (MQ=255)                                                                                                                                                  gggTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTaa  <  2:661264/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       ATATTTCCAAATTTATCAAGTGCTTAAATAATTAAATCTGTGCTAAAAACCAGGTAAGGATCAGTAGGTCAGCACTGCCGCCTGGACTGAGATTTCGTTCGATACACTCCCTGTCGAACTGCCGGAGATAATCGAGATCGGC‑GGGGGTTC‑GAATG‑CCCCCTTTTTGCAATAATGTTTGCGCCTCGCGCTGTAGCCAGCGCAGGCCCCCCTCGCCACCGCGCGATGCAACGTT‑GGTATCGCCGTTGATCGCCATCAGTAGGAGCAAGGTATCGAGCAATGCCAGTTCAGGATCTAA  >  NC_000913/646600‑646894

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