Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	41.9%	116.7 / 13.1	62	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (35/1);  new base C (0/26);  total (35/27)
Fisher's exact test for biased strand distribution p-value = 1.29e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591                                                                                                                                                        |                                                                                                                                                        gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCa‑cac                                                                                                                                                        <  1:1156934/151‑1 (MQ=255)   atatGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CAccc                                                                                                                                                      >  1:837192/1‑151 (MQ=255)       gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGTCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑c                                                                                                                                                 >  1:58144/1‑151 (MQ=255)           tGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑ccccc                                                                                                                                             >  1:1104210/1‑151 (MQ=255)               aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCccgccg                                                                                                                                         >  2:632204/1‑151 (MQ=255)               aaTGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCccgccg                                                                                                                                         >  1:503887/1‑151 (MQ=255)                    aaaaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATAc                                                                                                                                     >  2:428805/1‑150 (MQ=255)                            gggAGCCAAAAAATGGTATGCCGACTTGCCACCCATAATTATAAAAAAGATAAATGACGGCAGATATTTATCGGAAAGG‑AGAAAATGCAGGGTTTTTCTTCTCGTTCAATTACAAACCCCA‑CCCCCGGC‑CCCCCGCCGATACAAAGCGtt                                                                                                                             <  2:902081/150‑1 (MQ=255)                             ggcgCCCAGAAAATTTATGACGACTTCCAAGACATCATTATTAAAAAGATAAGTCACGCCAGATTTTTGTCGGTAAGG‑AGGAACGGAGCCTCTATTGCTCCCCATCAATGCTAAACGCCC‑CCCCCTGC‑CCCCCGCCGATACAAAGCGTTGc                                                                                                                           <  2:157090/148‑1 (MQ=255)                                 ccaaaaaaaTGTATGTCGACTGGCCACCCATCATTATACAACAGATAAAGGACGGCAGATATTTATGGGAAAGG‑AGAAAATGCAGGGCTTTTCGTCAC‑TTTAATGGGCAAACCCC‑CCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                       <  2:750792/146‑1 (MQ=255)                                   cAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                     >  2:902173/1‑151 (MQ=255)                                    aaaaaaTGTATGTCAACTTGCCACCCATCATTATTCAAAAGATTAAGGACGGAAATTTTTAATCGGAAGGG‑AGGAAATGAAGAGCTATTCGCAACGTTAAATGGTCAACCCCC‑CCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTcc                                                                                                                     <  2:674029/148‑1 (MQ=255)                                    aGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCt                                                                                                                    >  1:484288/1‑151 (MQ=255)                                      aaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  >  2:1150394/1‑151 (MQ=255)                                      aaaaTGTATGCCGAGTCGCCACCCATCATTATCAAAAAAATAAATGACGAAAAATATTTAGCGGAAAGG‑AGAAAAGGCAGAGCTATTGGTCCCGTTAATGGGAAAACCCCA‑CCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCtt                                                                                                                   <  2:495652/150‑1 (MQ=255)                                        aaTGTATGTCGACTTGCCCGCCATAATTATTCAAGAAAATAATGGCGGCAGTTATTCTGCGGTAAGA‑AGAAAAGGAAGGGCTATTGGTAACGTCCAATGCTAATCGCCA‑CACCCCGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGc                                                                                                                <  2:1012626/151‑1 (MQ=255)                                               gtCGACTGGCCAGCCATCATTATTAAAAAGATTAAGGACGGCAGATATTTACCGGAAAGG‑AGAAAAGGCAGAGTTATTCGTCACGTTAAATGGTAAATGCCA‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                            <  1:170864/148‑1 (MQ=255)                                                tCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                         >  2:921979/1‑150 (MQ=255)                                                 ctaCTTGCCAGCCATCATTATCAAACAGATTAAGGACGGCAGATATTTATCGGAAAGG‑AGAAAATGCAGACTTTTTCGTCACTTTCAATGGTCAATGCCC‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGc                                                                                                            <  1:832119/144‑1 (MQ=255)                                                 cGACTGGCCAGCAATCATTATCCAAAAAATAAAGGACGGCAGTTTTTTATCGGTAGGG‑AGAAAAGGAAGATTTTTTCTTCCCGTTAAAGGTTCAATCCCC‑CCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTAc                                                                                                       <  1:933342/151‑1 (MQ=255)                                                  gACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACg                                                                                                      >  2:510898/1‑151 (MQ=255)                                                         cacccatatttattaAACAGATTAATGCCGGCAGATTTTTACCGGTAAGG‑AGGAAATGCAGGGTTATTCGTCACTTTAAATGGTAATTGCCA‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  1:904848/140‑1 (MQ=255)                                                          agccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCAcc                                                                                               >  2:1175149/1‑150 (MQ=255)                                                           gccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAt                                                                                             >  1:620285/1‑151 (MQ=255)                                                            ccaTCATTATTCAACAGATTAATGACGGCAGATATTTATCGGCAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                            >  2:991434/1‑151 (MQ=255)                                                            ccaTCATTATTCAAAAGATTAATGCCGGCAGATATTTATCGGAAAGG‑AGGAAAAGCAGCCCTTTTCGTCACGTTCAAGGGTAAATGCCA‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  1:828005/146‑1 (MQ=255)                                                             catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                           >  2:197190/1‑151 (MQ=255)                                                             catcatTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAttt                                                                                           >  1:25969/1‑151 (MQ=255)                                                                  ttattaAAAAGATTAATGCCGGCAGATATTTATCGGTAAGG‑AGAAAATGCAGAGCTTTTCGTCACTTTAAATGGTAAATCCCC‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  1:233817/147‑1 (MQ=255)                                                                   tattCAACAGATTATTGACGGCAGATATTTACCGGTAAGG‑AGGACATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGA‑CCCCCTCCGACACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGCTTCACCATTTCATGa                                                                                      >  2:1191680/1‑150 (MQ=255)                                                                    attCAAAAGATAAATGACGGAAAATTTTTATC‑GTAAGGAAGAAAATGAAGAGTTTTTCTCCACTTTCATTGTCAAATCCCC‑CCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGa                                                                                      <  2:1195256/149‑1 (MQ=255)                                                                      tCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAAcc                                                                                   >  1:331715/1‑150 (MQ=255)                                                                                aagggggagagttttttGTCGATAAGG‑AGAAACTAAAGCCCTATTCGACACGTAAGAGGGCCAA‑GCCACCCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                         <  2:539762/136‑1 (MQ=255)                                                                                aaggACGGCATATTTTTTCGGGTAAGG‑AGAAAAGGCAGAGCTATTCGTCATGTTCAATGGTCAATCCCA‑CCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          <  2:349921/146‑1 (MQ=255)                                                                                aaTGACGGCAGATTTTTATCGGTAAGG‑AGAAAATGCAGAGCTTTTCGTCACGTTCAAGGGTAAACCCCA‑CCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAAcc                                                                          <  2:393552/149‑1 (MQ=255)                                                                                 aTGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAgg                                                                       >  2:447439/1‑151 (MQ=255)                                                                                      ggCAGATATTTTTCGGTAAGG‑AGAAAATGCAGAGCTTTTCTTCCCGTTAAAGGGTAAATGCCC‑CCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTc                                                                   <  1:272578/150‑1 (MQ=255)                                                                                       gaaGATATTTATCGGTAGGG‑AGAAAAGGCAGAGTTATTCGTCACGTTAAAGGGTAAAGGCCC‑CCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                                  <  1:677862/148‑1 (MQ=255)                                                                                        cAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCgg                                                                 >  1:196145/1‑150 (MQ=255)                                                                                        cAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  1:511613/1‑151 (MQ=255)                                                                                         aGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                >  2:395588/1‑150 (MQ=255)                                                                                              tttATCGGAAAGG‑AGAAAATGCAGAGCTTTTCGCCACGTAAAAGGATCAATCCCC‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                           <  2:290441/150‑1 (MQ=255)                                                                                                 aTCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                         >  2:861585/1‑149 (MQ=255)                                                                                                    ggTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcgc                                                    >  2:649841/1‑151 (MQ=255)                                                                                                    ggTAAGG‑AGAAAATGCAGAGCTTTTCGTCACGTTAAATGGTCAAGCCCC‑CCCCCTGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAgcg                                                     <  1:394760/150‑1 (MQ=255)                                                                                                     gTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGA‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGcc                                                   >  2:297125/1‑151 (MQ=255)                                                                                                        agg‑aggAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTAt                                                >  2:911884/1‑151 (MQ=255)                                                                                                         gg‑aggaAATGCAGAGCTTTTCGTCACTTTAAATGGTAAATGCCA‑CCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATc                                               <  1:179244/151‑1 (MQ=255)                                                                                                               aaaTGCAGATCTTTTCGTCCGGTTAAGGGGTCAACGCCC‑CCCCCGCC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCgg                                             <  2:725892/148‑1 (MQ=255)                                                                                                               aaaTGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          >  1:603750/1‑151 (MQ=255)                                                                                                               aaaGGCAGATTTATTCGTCACGTTAAATGGTAAATCCCC‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCggg                                            <  1:912608/149‑1 (MQ=255)                                                                                                                       agCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTg                                  >  2:1085685/1‑151 (MQ=255)                                                                                                                         cTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                >  2:1154182/1‑151 (MQ=255)                                                                                                                          tATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                >  2:177249/1‑150 (MQ=255)                                                                                                                            ttCGTCTCTTTCAATGGTAAATGCCC‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCg                                <  1:164038/148‑1 (MQ=255)                                                                                                                               gTCCCGTTCAATGGTCAATGCCA‑CCCCCGGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATc                             <  1:844307/148‑1 (MQ=255)                                                                                                                                    gTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  1:972313/1‑151 (MQ=255)                                                                                                                                    gTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     >  2:412925/1‑151 (MQ=255)                                                                                                                                        aaTGGTCACTGCCC‑CCCCCGGC‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     <  1:197194/147‑1 (MQ=255)                                                                                                                                        aaTGGTCAATGCCA‑CCCCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     <  1:1191216/147‑1 (MQ=255)                                                                                                                                             tCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt            >  2:831603/1‑151 (MQ=255)                                                                                                                                              cAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGAATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  1:688947/1‑151 (MQ=255)                                                                                                                                               aaTGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacattg           >  1:39836/1‑150 (MQ=255)                                                                                                                                                        acCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:615793/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                        GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGTAAGG‑AGGAAATGCAGAGCTATTCGTCACGTTCAATGGTCAATGCCA‑CACCCTGG‑CCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591

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