Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	446,976	0	A	C	41.9%	90.4 / 20.9	62	W244G (TGG→GGG)	cyoE	protoheme IX farnesyltransferase
Reads supporting (aligned to +/- strand):  ref base A (27/9);  new base C (0/26);  total (27/35)
Fisher's exact test for biased strand distribution p-value = 3.45e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATAC  >  NC_000913/446828‑447123                                                                                                                                                      |                                                                                                                                                     gccagcagcGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAaccacc                                                                                                                                                    <  1:300942/151‑1 (MQ=255)               cGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAAccgcc                                                                                                                                     >  1:63139/1‑151 (MQ=255)                    ccGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCgcggc                                                                                                                                >  1:458976/1‑151 (MQ=255)                          ccATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGacc                                                                                                                            <  1:994608/149‑1 (MQ=255)                           cATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGaccacc                                                                                                                         >  2:801491/1‑151 (MQ=255)                              aaaaTCACGGGAGATCATGACCGGGGGGGCAGTGAGGGCGGTGAGGGAGAGCCCGAACAGCTGGCGCGCCCAGATTCTGTCAACCGCACCTTCATAG‑CCGCGGAGAGCCATGCCAAACCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGa                                                                                                                      <  2:84585/151‑1 (MQ=255)                              aaaaTCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGaccacc                                                                                                                         >  2:349963/1‑148 (MQ=255)                              aaaaTAAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACGGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTACCCACCAG‑ACTCTAACCGCCGCGGCGACCCCCAg                                                                                                                       <  2:967363/150‑1 (MQ=255)                                aaTCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGata                                                                                                                    <  2:7316/151‑1 (MQ=255)                                      cGGACATCATCCCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAC‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTtat                                                                                                               <  1:861572/150‑1 (MQ=255)                                      cGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCAAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAA‑ACGGTAACCGCCGCGGCGACCACCAGATATTtata                                                                                                              <  1:828208/151‑1 (MQ=255)                                       ggACATCATCACCGAGGGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATAc                                                                                                             <  2:313970/151‑1 (MQ=255)                                        gACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATAcc                                                                                                            >  2:566767/1‑151 (MQ=255)                                          catcatCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATAccc                                                                                                           >  2:535264/1‑150 (MQ=255)                                                caCCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAC‑CCGCGCAGAGCCATACCTACCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTaa                                                                                                    <  2:656332/151‑1 (MQ=255)                                                    gagaGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCg                                                                                                 >  2:379538/1‑150 (MQ=255)                                                             gTGATGGCGATGATAGAGAAGCCAAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATCCCTACCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAagagaga                                                                                       <  1:765821/151‑1 (MQ=255)                                                                  ggCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATc                                                                                  >  2:891058/1‑151 (MQ=255)                                                                  ggCGATGAAAGAGAAGCCGAACAGCTTGCGCCCCCAGATTCTGTCATCAGCAACTTTATAC‑CCGCGCAGAGCCATCCCAAACCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAagagagag                                                                                      <  1:252/147‑1 (MQ=255)                                                                       tgagAGAGAAGCGAAACAGCTGGCGCGCCCATATTCGGTCATCACAAACTTAAAAC‑CCGCGAAGAACCATCCCTACCCCCCAG‑ACGCAAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCa                                                                                 <  2:850007/147‑1 (MQ=255)                                                                         atAGAGAACCCAAACAGCTTGCGCCCCCAGATTCTGTCATCAGCAACTTTATAC‑CCGCGCAGAGCCATCCCAACCCCCCAG‑CCGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTgg                                                                           <  1:629011/151‑1 (MQ=255)                                                                          tAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGc                                                                          >  2:483323/1‑151 (MQ=255)                                                                           agagAAGCCAAACAGCGTGCGCGCCCAGTTTCTGTCTACAGCACCTTTATAA‑CCGCGCAGAGCCATCCCGAACCCCCAG‑ACGCAAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCa                                                                         <  2:1036669/151‑1 (MQ=255)                                                                               aacccGAACAGCTGGCGCCCCCAGATTCTGTCATCAGCAACTTAAAAC‑CCGCGAAGAGCCATCCCAACCCCCCAG‑ACGCAACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACg                                                                      <  2:463765/147‑1 (MQ=255)                                                                               aaGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCACCTTTATAA‑CCGCGCAGAGCCATACCAACCCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAAc                                                                       <  1:72959/149‑1 (MQ=255)                                                                               aaGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACg                                                                      >  2:458740/1‑150 (MQ=255)                                                                               aaGCCAAAAAGCTGGCGCCCCCAGATTCTGTCATCACAAACTTTATAC‑CCGCGCAGAGCCATCCCAACCCCCCAG‑ACGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAAc                                                                       <  2:1024761/149‑1 (MQ=255)                                                                                    gAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCaa                                                                  >  2:253408/1‑149 (MQ=255)                                                                                     aaCAGCTTGCGCCCCCAGATTCTGTAATCAGAACCTTTAAAC‑CCGCGCAGGGCCATCCCAACCCACCAG‑CCGAAAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAg                                                                <  1:891383/150‑1 (MQ=255)                                                                                        agtttGCGCGCCCAGATTCTGTCACCAGCACCTTTAAAC‑CCGCGAAGAGCCATCCCTACCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCg                                                             <  2:520490/147‑1 (MQ=255)                                                                                        aGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCg                                                             >  1:706306/1‑150 (MQ=255)                                                                                         gCTTGCGCCCCAATATTCGGTCACCAGCAACTTAAAAC‑CCGCGC‑GAACCCAGCCCAACCCCCAGCCCCCAAACCGCCGCGGCGACCACCAGATATTTAACCCCGGCGTATCCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCga                                                            <  2:473881/150‑1 (MQ=255)                                                                                            tGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatgat                                                        >  1:889444/1‑151 (MQ=255)                                                                                             gcgcgcCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatgat                                                        >  2:848465/1‑150 (MQ=255)                                                                                                 gcCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATAc                                                    >  1:340014/1‑150 (MQ=255)                                                                                                 gcCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATAc                                                    >  2:747788/1‑150 (MQ=255)                                                                                                 gcCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATAc                                                    >  1:638698/1‑150 (MQ=255)                                                                                                  cccAGTTTCTGCCACCAG‑AACTTTAAAACCCGCGCAGAGCCTTCCCTACCCACCAG‑ACGCAAACCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCgatg                                                          <  2:594024/143‑1 (MQ=255)                                                                                                   ccAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGc                                                 >  1:616233/1‑151 (MQ=255)                                                                                                     aGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCgt                                               >  2:877698/1‑151 (MQ=255)                                                                                                        ttCTGTCATCAGCAACTTTATAA‑CCGCGCAGACCCATACCTACCCCCCAG‑CCGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCgt                                               <  1:287586/148‑1 (MQ=255)                                                                                                          cTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCgtgatg                                           >  2:180554/1‑150 (MQ=255)                                                                                                                  aGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAttcttc                                  >  2:451252/1‑151 (MQ=255)                                                                                                                    cAACTTTATAC‑CCGCGAAGGGCCATCCCTACCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAtt                                      <  1:360898/145‑1 (MQ=255)                                                                                                                    cAACTTAATAC‑CCGCGCAGAGCAATACCAAACCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCgtgatgtgat                                       <  2:300154/144‑1 (MQ=255)                                                                                                                     aaCTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGcc                               >  2:434274/1‑151 (MQ=255)                                                                                                                     aaCTTAATAC‑CCGCCCAGACCAATCCCTACCCCCCAG‑ACCCAAACCGCCGCGGCGACCACCAGATATTTATCCCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGAttcttc                                  <  1:79215/148‑1 (MQ=255)                                                                                                                      aCTTTATAA‑CCGCGAAGAGCAATCCCTACCCCCCGG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGcc                               <  2:821794/150‑1 (MQ=255)                                                                                                                        tttATAC‑CCGCGCAGACCCATCCCTACCCCCCAG‑CCGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAc                             <  1:820945/150‑1 (MQ=255)                                                                                                                         ttATAC‑CCGCGCAGAGCCATCCCAACCCCCCAG‑CCGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCAc                             <  1:474403/149‑1 (MQ=255)                                                                                                                           ataA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGa                          >  1:41565/1‑150 (MQ=255)                                                                                                                           ataA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGa                          >  2:701389/1‑150 (MQ=255)                                                                                                                             aa‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAAt                       >  1:677963/1‑151 (MQ=255)                                                                                                                             aa‑CCGCGCAGACCCATCCCTACCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGa                          <  1:591839/148‑1 (MQ=255)                                                                                                                              accgccgagggcccATCCCAACCCCCCAG‑CCGCAACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAAt                       <  1:954143/139‑1 (MQ=255)                                                                                                                                ccGCGCAGACCCATCCCAACCCACCGG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGc                     <  1:60360/151‑1 (MQ=255)                                                                                                                                 cgcgCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGcc                    >  2:596362/1‑151 (MQ=255)                                                                                                                                   cgcAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCtt                  >  1:176557/1‑151 (MQ=255)                                                                                                                                       gacccATCCCTAACCCCCAA‑CCGCAACCCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTAc              <  1:919643/148‑1 (MQ=255)                                                                                                                                         cccATACCAACCCCCCAG‑CCGCAACCCGCCGCGGCGACCCCCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCa            <  2:259854/150‑1 (MQ=255)                                                                                                                                          cctacccaaaCCCCCGG‑ACTCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCAc           <  2:147062/143‑1 (MQ=255)                                                                                                                                              cccaaaCCCCCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGc       <  1:371171/147‑1 (MQ=255)                                                                                                                                               ccTACCCCCCAG‑ACGCTACCCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGc       <  1:567033/150‑1 (MQ=255)                                                                                                                                                   accaccAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATAc  >  2:1032605/1‑151 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GCCAGCAGCGTATGCGAGTCCGGTACCATAAAATCAACGGACATCATCACCGAGAGGGCAGTGATGGCGATGATAGAGAAGCCGAACAGCTTGCGCGCCCAGATTCTGTCATCAGCAACTTTATAA‑CCGCGCAGAGCCATACCTAACCACCAG‑ACGCTAACCGCCGCGGCGACCACCAGATATTTATACCCAGCGTAACCGCCAAGAGAGAGCATCAGCGTGGCAACGGCAAAGGCGATGATATACAGCGTGATGTGATTCTTCGCCACCGAAATGCCTTTTACCACTGGCAATAC  >  NC_000913/446828‑447123

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