Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,859,828	0	T	G	51.5%	‑1.6 / 114.6	66	I23M (ATT→ATG)	pphB	serine/threonine‑specific protein phosphatase 2
Reads supporting (aligned to +/- strand):  ref base T (14/18);  new base G (14/20);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.85e-01

CGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTAAACCAACCCTGGTTTACG  >  NC_000913/2859682‑2859975                                                                                                                                                   |                                                                                                                                                    cGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATg                                                                                                                                                 <  2:2953/151‑1 (MQ=255)    gAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTg                                                                                                                                              <  1:925658/151‑1 (MQ=255)        atTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAat                                                                                                                                           >  2:752474/1‑150 (MQ=255)           aaCGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATc                                                                                                                                        <  2:762869/150‑1 (MQ=255)           aaCGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATc                                                                                                                                        <  2:557392/150‑1 (MQ=255)              gACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGtt                                                                                                                                    >  1:107249/1‑151 (MQ=255)                 aaCTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGttatt                                                                                                                                 <  2:549675/151‑1 (MQ=255)                     gACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACaa                                                                                                                             <  2:451347/151‑1 (MQ=255)                        gAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGAGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGACCATCACTATAGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATcc                                                                                                                          >  1:449427/1‑151 (MQ=255)                          ataaaataaaAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCg                                                                                                                        <  2:153233/151‑1 (MQ=255)                            aaaataaaaACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGc                                                                                                                       <  2:322863/150‑1 (MQ=255)                            aaaataaaaACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCt                                                                                                                      <  1:767527/151‑1 (MQ=255)                            aaaataaaaACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGc                                                                                                                       <  1:714233/150‑1 (MQ=255)                                 aaaaaCACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACa                                                                                                                  <  1:668022/150‑1 (MQ=255)                                  aaaaCACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATc                                                                                                                >  2:89315/1‑151 (MQ=255)                                  aaaaCACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATc                                                                                                                <  1:408965/151‑1 (MQ=255)                                    aaCACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCa                                                                                                               >  2:151839/1‑150 (MQ=255)                                            gTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCttt                                                                                                      >  1:465037/1‑151 (MQ=255)                                              ataataGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCttttt                                                                                                    >  1:488806/1‑151 (MQ=255)                                              ataataGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCttttt                                                                                                    >  1:222534/1‑151 (MQ=255)                                               taataGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCAATATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCACGATTACATCAACTCTAtttttt                                                                                                   >  2:454241/1‑151 (MQ=255)                                                aataGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTc                                                                                                  >  1:660418/1‑151 (MQ=255)                                                 ataGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTcc                                                                                                 >  2:846645/1‑151 (MQ=255)                                                   aGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTGCATCAACTCTCTTTTTTcccc                                                                                               >  1:669553/1‑151 (MQ=255)                                                   aGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTccc                                                                                                <  1:423384/150‑1 (MQ=255)                                                      aaaGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCaaa                                                                                            <  2:506468/151‑1 (MQ=255)                                                          cTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCg                                                                                        >  2:8529/1‑151 (MQ=255)                                                             aTTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGac                                                                                      <  1:93824/150‑1 (MQ=255)                                                                 tACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGacttact                                                                                 >  2:708993/1‑151 (MQ=255)                                                                  aCAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGacttactt                                                                                <  2:369377/151‑1 (MQ=255)                                                                  aCAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGacttact                                                                                 <  1:547659/150‑1 (MQ=255)                                                                    aGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTAt                                                                              >  2:865309/1‑151 (MQ=255)                                                                    aGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTAt                                                                              <  2:149696/151‑1 (MQ=255)                                                                    aGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTAt                                                                              <  1:232007/151‑1 (MQ=255)                                                                       gTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTc                                                                           <  1:763970/151‑1 (MQ=255)                                                                         aaaaCCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCt                                                                          <  1:722384/150‑1 (MQ=255)                                                                             ccATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCgg                                                                     <  1:789186/151‑1 (MQ=255)                                                                              cATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGc                                                                    >  2:387389/1‑151 (MQ=255)                                                                                 gCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGa                                                                  <  1:517997/150‑1 (MQ=255)                                                                                  ccATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGat                                                                 >  2:56151/1‑150 (MQ=255)                                                                                   cATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGata                                                                >  1:701051/1‑150 (MQ=255)                                                                                        acacGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATAtt                                                           <  2:280565/150‑1 (MQ=255)                                                                                          acGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGAt                                                        <  2:778556/151‑1 (MQ=255)                                                                                          acGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTTTCGTCGATAATATTGAt                                                        >  1:54137/1‑151 (MQ=255)                                                                                          acGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGAt                                                        <  1:646623/151‑1 (MQ=255)                                                                                           cGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATc                                                       <  1:208420/151‑1 (MQ=255)                                                                                           cGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATc                                                       >  1:425319/1‑151 (MQ=255)                                                                                           cGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGAt                                                        >  2:195882/1‑150 (MQ=255)                                                                                               aTCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTgg                                                   <  1:291447/151‑1 (MQ=255)                                                                                                    aaaaTCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCgg                                              <  2:592115/151‑1 (MQ=255)                                                                                                         cAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGt                                         <  1:873026/151‑1 (MQ=255)                                                                                                           aTGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCt                                       >  2:801814/1‑151 (MQ=255)                                                                                                             gCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCtt                                      <  1:752455/150‑1 (MQ=255)                                                                                                             gCCCATCACTATCGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTg                                     <  2:640047/151‑1 (MQ=255)                                                                                                                        tCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTgcg                           <  1:925978/150‑1 (MQ=255)                                                                                                                         cGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGcc                         >  2:552559/1‑151 (MQ=255)                                                                                                                         cGCCATATATGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGcc                         >  1:706250/1‑151 (MQ=255)                                                                                                                          gCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCt                        <  2:94758/151‑1 (MQ=255)                                                                                                                           ccATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTg                       >  2:6859/1‑151 (MQ=255)                                                                                                                               atatGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTa                    <  2:137793/150‑1 (MQ=255)                                                                                                                               atatGGGTCGTTGGTGATATGCATGGGGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTa                    >  2:578788/1‑150 (MQ=255)                                                                                                                                tatGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTaa                   >  1:505213/1‑150 (MQ=255)                                                                                                                                 atGGGTCGTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTaaa                  <  1:460725/150‑1 (MQ=255)                                                                                                                                     gTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTAaaccaac             <  1:89313/151‑1 (MQ=255)                                                                                                                                        gTTGGTGATATGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTAAACCAACCCt          <  2:840387/151‑1 (MQ=255)                                                                                                                                                tatGCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTAAACCAACCCTGGTTTACg  >  1:476344/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CGGGAATATTAACGACAACTGACGAATAAAATAAAAACACCCTGTATAATAGGAAAGCTTATTTTACAGGGTAAAACCATGCCATCTACACGCTATCAAAAAATCAATGCCCATCACTATCGCCATATATGGGTCGTTGGTGATATTCATGGTGAATATCAGTTATTACAATCCCGCTTACATCAACTCTCTTTTTTCCCCAAAATCGACTTACTTATTTCTGTCGGCGATAATATTGATCGTGGACCGGAGAGTCTTGACGTCCTGCGCCTGCTAAACCAACCCTGGTTTACG  >  NC_000913/2859682‑2859975

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