Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,984,441	0	A	C	57.6%	46.5 / 87.5	60	G384G (GGT→GGG)	yqeF	short chain acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (21/4);  new base C (0/34);  total (21/39)
Fisher's exact test for biased strand distribution p-value = 2.44e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.04e-10

GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591                                                                                                                                                        |                                                                                                                                                       gAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCacac                                                                                                                                                       >  2:56540/1‑151 (MQ=255)     atGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTg                                                                                                                                                   >  2:204975/1‑151 (MQ=255)       gTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGc                                                                                                                                                 >  1:1084120/1‑151 (MQ=255)         cATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGccc                                                                                                                                               <  2:680437/151‑1 (MQ=255)          aTGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGcccc                                                                                                                                              <  1:534313/151‑1 (MQ=255)                      aaaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACaa                                                                                                                                   >  1:8227/1‑150 (MQ=255)                       aaTGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACa                                                                                                                                    >  1:540806/1‑148 (MQ=255)                            gggAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTg                                                                                                                            >  2:360970/1‑151 (MQ=255)                                 cccAAAAAATTTATGTCGTCTTGCCAGCCATCATTATTAAACAGATTAATGCCGGCAAATTTTTACGGTT‑AAGGAGGAAATGAAGAGCTTTTCGTCACGTTAAATGG‑AAAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAg                                                                                                                        <  1:614444/150‑1 (MQ=255)                                 cccAAAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATAAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGTTTTTCGTCACTTTCAATGG‑TAAATGCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  2:752920/149‑1 (MQ=255)                                    aaaaaaTGTTTGTTGACTTGCCACCCATCATTATTCAACAAATAAATGACGGCAGATATTTATCGGA‑AAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGG‑TCAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGt                                                                                                                       <  1:924787/146‑1 (MQ=255)                                    aaaaaaTGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCATATTTTTATCGGT‑AAGGAGAAAATGCAGATCTTTTCGTCCCGTTAAGTGG‑TCAACCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCa                                                                                                                         <  1:171320/144‑1 (MQ=255)                                      aaagTGTATGGCGCTTTCCCAGCCATCATTATTCAACAGATTAATGACGGCATTTATTTATCGGT‑AAGGAGAAAATGCAGACTTTTTCGCCACGTTAAATGG‑TAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCttt                                                                                                                  <  1:141523/151‑1 (MQ=255)                                           ttatgtCGTCTTGCCACCCATTATTTTTAAACAAATAAAGGACGGCAAATATTTACCGGT‑AAGGAGAAAATGCAGAGTTTTTCGTCCCGTTCAAGGG‑TCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               <  1:1253602/148‑1 (MQ=255)                                            tatgtCGACTTGCCACCCATCATTATTCAAAAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTTTTCGTCACGTTAAAGGG‑CCAATCCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCg                                                                                                               <  1:807341/148‑1 (MQ=255)                                                 cgattggccagccaTCATTATTCAACAGATAAATGACGGCAGATATTTATCGGT‑AAGGAGAAAATGCAGATTTTTTCGTCCCTTTCAATGG‑TAAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTa                                                                                                        <  2:344188/144‑1 (MQ=255)                                                   aCTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGt                                                                                                     >  2:332898/1‑151 (MQ=255)                                                     tggccgcccatttttattCAAAAGATAAAGGACGGAGAAAATTTCTGGGT‑AGGGAGAAACGGCAGAGCTTTTCGCCACGTT‑AATGGATCAATGCCCCCCCCGGCCCGCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCAtt                                                                                                         <  1:497542/132‑1 (MQ=255)                                                        ccagccaTCATTATTCAACAGATAAATGACGGCAGATTTTTATCGGT‑AGGGAGGAAATGAAGAGCTATTCGTCACGTTAAATGG‑TAAATGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGttt                                                                                                   <  2:1057009/148‑1 (MQ=255)                                                        ccagccaTAATTATTCAACAGATTAATGACGGCAGATTTTTATGGGT‑AAGGAGAAAATGCAGAGTTTTTCGTCCCGTTCAATGG‑TAAATGCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCa                                                                                                 <  2:1050401/150‑1 (MQ=255)                                                              atcatTATACAACAGATTAATGACGGCAGATTTTAATGGGT‑AAGGAGAAAATGAAGAGTATTTCGTCGCTTTCAATGG‑CAAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCAtt                                                                                            <  2:206664/149‑1 (MQ=255)                                                                catTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCAt                                                                                        >  1:837695/1‑151 (MQ=255)                                                                  ttattaAACAGATTAAGGACGGCAGATTTTTATCGGT‑AAGGAGGAAATGCAGATCTATTCGTCACGTTAAATGG‑TAAATGCCCCGCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       <  2:819652/150‑1 (MQ=255)                                                                  ttattaAAAAAATAAAGGACGGAAGATTTTTTGCGGT‑AAGGAGAAAAGGCAGTGCTTTTCGCACCGTTAAATGG‑TCACCCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTc                                                                                          <  1:23143/147‑1 (MQ=255)                                                                  ttattCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       >  2:1084888/1‑150 (MQ=255)                                                                  ttataCAACAGATAAATGACGGCAGATATTTATCGGT‑AAGGAGAAAATGCAGAGCTTTTCGTCACGTTCAATGG‑TCAATGCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATg                                                                                       <  1:1309574/150‑1 (MQ=255)                                                                     tacaaaaaataaaTGAGGGAAGATATTTA‑CGGTAAAGGAGAAATGGAAGGCCTTTTGGTCGCGTTCAATGG‑TAATTCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGaa                                                                                     <  1:347840/138‑1 (MQ=255)                                                                        aaaaGATTAATGACGGCAGATATTTAGCGGT‑AAGGAGGAAAGGCAGAGCTTTTCGTCACGTTAAATGG‑TCAATCCCACCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                 <  2:455508/147‑1 (MQ=255)                                                                        aaCAGATTAATGCCGGAATATTTTTATGGGT‑AGGGAGGAAAGGAAGAGCTATTCGTCACTTTAAATGG‑TCACTCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCa                                                                                  <  2:459334/149‑1 (MQ=255)                                                                        aaCAGATTAATGACGGAAAATATTTATCGGT‑AGGGAGGAAATGCAGAGCTATTCATCACGTTCATTGC‑CCAACCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCag                                                                                 <  1:9243/150‑1 (MQ=255)                                                                           aGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAgaga                                                                              >  1:302692/1‑150 (MQ=255)                                                                             aTTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGTTTTTCGTCACGTTAAATGG‑TCAACCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaaa                                                                            <  1:915492/150‑1 (MQ=255)                                                                             aTTAATGACGGAAGATATTTAGGGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCCCACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGaa                                                                             <  1:310916/149‑1 (MQ=255)                                                                                attggcgggatatatTTATGGGG‑AGGGAGGAAATGCGATTTTTTTGGTCAGGTACAGTGT‑TCAGTGCCACCCCTTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                         <  1:873899/139‑1 (MQ=255)                                                                                aaTGACGGAAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCGCGTTAAAGGG‑TCAATGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAg                                                                        <  2:1237561/151‑1 (MQ=255)                                                                                  tGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCa                                                                         >  1:562215/1‑148 (MQ=255)                                                                                      ggaaaatatTTATCGGA‑AGGGAGAAAATGCAGAGTTTTTCGTCCCTTTAAATGG‑TCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCg                                                                  <  1:703193/146‑1 (MQ=255)                                                                                      ggCAGATATTTATCGGT‑AAGGAGGAAAGGCAGAGTTTTTCGTCACGTTCAAGGG‑TCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGAt                                                                     <  2:198570/148‑1 (MQ=255)                                                                                         aaatatTTATGGGA‑AAGGAGAAAAGGCAGAGTTTTTCGTCACTTTAAAGGG‑TCAATGCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGc                                                                <  2:7988/148‑1 (MQ=255)                                                                                              tttATCGGT‑AAGGAGAAAATGCAGACTTTTTCGCCACTTTCAATGT‑TCAATCCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCg                                                           <  2:875448/150‑1 (MQ=255)                                                                                                tATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGa                                                         >  1:292741/1‑150 (MQ=255)                                                                                                    ggT‑AAGGAGAAAATGCAGACCTTTTCGTCACGTTCAATGG‑TCAAGGCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGaa                                                        <  2:1269718/147‑1 (MQ=255)                                                                                                               aaaTGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          >  2:61278/1‑151 (MQ=255)                                                                                                               aaaGGCAGAGTTTTTCGCCACTTTAAATGG‑TCAATCCCACCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCg                                              <  1:700872/147‑1 (MQ=255)                                                                                                                aaTGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTg                                          >  2:412936/1‑150 (MQ=255)                                                                                                                    cAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:1083452/1‑150 (MQ=255)                                                                                                                    cAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  2:100968/1‑150 (MQ=255)                                                                                                                    cAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:1165927/1‑150 (MQ=255)                                                                                                                     agagCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:1116791/1‑149 (MQ=255)                                                                                                                     agagCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      >  1:981111/1‑149 (MQ=255)                                                                                                                      gagTTATTCGTCACGTTAAATGG‑TCAATCCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGa                                     <  2:464442/149‑1 (MQ=255)                                                                                                                       agttATTCGTCACGTTAAATGG‑TCAATGCCCCCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCg                                      <  1:610182/144‑1 (MQ=255)                                                                                                                            ttCGTCGCGTTAAATGG‑TCAACCCCACCCCCGGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGAt                              <  2:1165125/150‑1 (MQ=255)                                                                                                                            ttCGTCACTTTAAATGG‑TCAATCCCACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGa                               <  2:694123/149‑1 (MQ=255)                                                                                                                                     ttaaaTGG‑TAAATCCCCCCCCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGccacc                     <  2:296009/147‑1 (MQ=255)                                                                                                                                        aaTGG‑TAAATGCCCCCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcatt                  <  2:1188132/150‑1 (MQ=255)                                                                                                                                           gg‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac               >  1:841675/1‑150 (MQ=255)                                                                                                                                           gg‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgac               >  2:28009/1‑150 (MQ=255)                                                                                                                                                 aTGCCCCCCCCGGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACcattgacatt            <  1:1117141/148‑1 (MQ=255)                                                                                                                                                   gCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGAccc       >  1:149579/1‑151 (MQ=255)                                                                                                                                                    ccACCCCCTGCCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCattgacattgac         <  2:1088276/148‑1 (MQ=255)                                                                                                                                                        acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:1244195/1‑151 (MQ=255)                                                                                                                                                        acCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGc  >  2:50682/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       GAATATGTCATGAGAATGGAAAAATGAGGGAGCCCAGAAAATGTATGTCGACTTGCCAGCCATCATTATTCAACAGATTAATGACGGCAGATATTTATCGGT‑AAGGAGGAAATGCAGAGCTATTCGTCACGTTCAATGG‑TCAATGCCACACCCTGGCCCCCGCCGATACAAAGCGTTGCCAGTCCTTTGCGGGCATTACGTTTCACCATTTCATGAACCAGAGAAACCAGGATTCGGCAACCGGAAGCGCCTATCGGGTGACCGAGTGCGATCGCGCCACCATTGACATTGACCCGACGC  >  NC_000913/2984292‑2984591

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