Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,434,882	0	T	G	41.7%	93.8 / 22.2	60	Y252S (TAC→TCC)	truA	tRNA pseudouridine(38‑40) synthase
Reads supporting (aligned to +/- strand):  ref base T (3/32);  new base G (25/0);  total (28/32)
Fisher's exact test for biased strand distribution p-value = 6.31e-14
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434735‑2435031                                                                                                                                                    |                                                                                                                                                       ggTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAAt                                                                                                                                                     <  2:139465/151‑1 (MQ=255)        atTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATTATACCGGGCAGGGGAATACAcc                                                                                                                                                >  2:1168569/1‑149 (MQ=255)             ccTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACCGTCAAGGGAATCTACcgcg                                                                                                                                             >  1:736775/1‑147 (MQ=255)               tATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATAATACCCGTCAGGGGAATCCACcccgcc                                                                                                                                           >  2:905369/1‑145 (MQ=255)                aTGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGa                                                                                                                                       <  2:1073568/150‑1 (MQ=255)                  gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATa                                                                                                                                     >  1:1085006/1‑150 (MQ=255)                  gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATa                                                                                                                                     <  1:901017/150‑1 (MQ=255)                  gCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATAc                                                                                                                                    <  1:93689/151‑1 (MQ=255)                   ccTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGGTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTTGGATATCATAACCGACAAGGGAATCCACCCCGACAAGaaa                                                                                                                                     >  1:1032074/1‑147 (MQ=255)                       gTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGGTTTGTTAGATGATAACAGGGGGGGTAAACCAACCCGCCCCGAAACCGccc                                                                                                                               >  1:960067/1‑151 (MQ=255)                         aCTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCt                                                                                                                              <  2:443874/150‑1 (MQ=255)                                    ttCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTg                                                                                                                    <  2:447749/149‑1 (MQ=255)                                      cAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATAACGGGCAGGGGAATCCACCGCGACCAGAAAAAGCCCTTCCCCTTTTTcc                                                                                                                  >  1:569667/1‑149 (MQ=255)                                       aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGTAAATCATAACGGGCAGGGGAATACACCGCGTCCCGATAAAACCCCTCCCCCTTTGcccgt                                                                                                               >  1:1113446/1‑148 (MQ=255)                                       aGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGGCAGGGGAATCCACCGCCACCCGATAAAGCCCTTCCCCTTTTTcc                                                                                                                  >  2:146213/1‑148 (MQ=255)                                        ggCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATAACGGTCAGGGGAAACCACCGCGACCAGATACAGCCCTTCCCCTTTTTCCGtt                                                                                                               >  1:397295/1‑150 (MQ=255)                                         gCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCACGACCAGATACAGCCCTTCCGCTTTTGCCGTtg‑c                                                                                                            <  1:382359/151‑1 (MQ=255)                                            gTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTtg‑ct                                                                                                           <  1:333852/149‑1 (MQ=255)                                                 tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCGGTTAGGGGAATCCACCGCGCCCAGATACAACCCTTCCGCTTTTGTCGTTG‑CTgccgcc                                                                                                     >  1:956180/1‑150 (MQ=255)                                                 tACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGCTTTTTAAGATCATAACGGGCCGGGGAAACCACCCCGGCCCGATAAAGCCCCTCCGCCTTTTCCGCTG‑CTgccgccc                                                                                                    >  2:239961/1‑150 (MQ=255)                                                  aCGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCa                                                                                                    <  2:526810/150‑1 (MQ=255)                                                        gCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCt                                                                                             <  1:371155/151‑1 (MQ=255)                                                           aCAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGt                                                                                           <  2:403827/150‑1 (MQ=255)                                                            cAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGACCCATTGGCGGTTTTTGAAGATCATAACGGGCCGGGGAATCCACCCCCACCAGATACAGCCCCTCCCCTTTTTCCCTTG‑CTGCCGCCCGCGGTCTGTcc                                                                                         >  2:201898/1‑151 (MQ=255)                                                               ccTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCttt                                                                                      <  2:821499/151‑1 (MQ=255)                                                                cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCAAACCCGGCAGGGGAATCCACCGCGAACCGATACAACCCCTCCGCCTTTTCCGGTG‑CTGCCCCCAGCGTTCTGTCCttt                                                                                      >  2:649881/1‑150 (MQ=255)                                                                cTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGGAGATCATACCCGGCAGGGGAAAACACCGCGGCCCGAAACAGCCCTTCCGCTTTTTTCGTTG‑CTGCCCCCCGCGTTTCGTCCttt                                                                                      >  2:934624/1‑150 (MQ=255)                                                                cTGAAATGCTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCCGGCCGGGGAATCCACCCCGACCCGATACAGCCCTTCCGCTTTTTCCGTTT‑CTGCCCCCCGCGGTCTGTCCttt                                                                                      >  2:351130/1‑150 (MQ=255)                                                                    aaTGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGGCAGGGGAATCCACCGCGGCCAGATACAGCCCTTCCCCTTTTTCCGTTG‑CTGCCGCCAGCGGTCTGTCCttt                                                                                      >  1:820202/1‑146 (MQ=255)                                                                             aGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATGAAACCGGTCAGGGGAATACACCGGGGCCAGATATAGCCCCTCCGACTTTTTCGCTG‑CAGGCGCCAGGGGTCTGTCCTTTTCAGCCAGCAActc                                                                        >  1:1068232/1‑151 (MQ=255)                                                                                attattCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGGCAGGGGAAACCACCCCGCCCCGATAAAACCCTTCCGCTTTTTCCGGTG‑CTGCCCCCCGCGGTCTTTCCTTTTCCCCCCGCaa                                                                           >  2:329446/1‑145 (MQ=255)                                                                                        cGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAg                                                              <  2:696927/150‑1 (MQ=255)                                                                                         gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCCCACAGGGGAATACACCCCCACCAGATACAGCCCCTCCGCTTTTTCCGCTGTCTCCCCCCA‑CGCTCTGTCCTTTTCCGCCCACAACACTGCTAttc                                                                >  2:92056/1‑145 (MQ=255)                                                                                         gTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGct                                                            <  1:94129/151‑1 (MQ=255)                                                                                            aGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATAACGGGCAGGGGAATCCACCCCGACCAGATACAACCCCTCCCCTTTTTCCGTTG‑CTGCCCCCAGCGGTCTGTCCTTTTCCGCCAGCAAATCTGCTATCCAGCtctc                                                         >  1:975876/1‑151 (MQ=255)                                                                                             gTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTcc                                                        <  2:1080533/151‑1 (MQ=255)                                                                                                  ccAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTg                                                   <  2:262181/151‑1 (MQ=255)                                                                                                   cAGAAATAGCGGGCCCATTGGCGGTTTTTGGAGATCATACCGGTCAGGGGAATCCACCCCGACCAGATACAGCCCCTCCGCTTTTTCCGTTT‑CTGCCCCCAGCGGTCTGTCCCTTTCCGCCAGCAACTCTTCTTTCCAGCTCTCCCGCTgg                                                  >  1:311363/1‑151 (MQ=255)                                                                                                    aGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGt                                                 <  2:719212/151‑1 (MQ=255)                                                                                                     gAAATAGCGGGCCCATTGGCGGTTTTTGAAGATTATACCGGTCAGGGGAAACCACCGCGACCCGATAAAGCCCTTCCCCTTTTTCCGTTT‑CTGCCGCCCGCGGTCTTTCCTTTTCCGCCCGCCAATCTGCTTTCCCGCTCTCCCGCTgggt                                                >  2:747779/1‑149 (MQ=255)                                                                                                     gAAATAGCGGGCCCATTGGCGGTTTTGGGAGAACATAACGGTCAGGGGAAATCACCCCGACCAGATACAGCCCCTCCGCCTTTGCCGTTG‑CTGCCGCCAGCCGTCTGTCCCTTGCCGCCCGCAAACCCGCTATTCAGCCCTCCCGGTggg                                                 >  2:968096/1‑149 (MQ=255)                                                                                                     gAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCCGGCAGGGGAATCCACCCCGACCAGATAAAACCCCTCCCCTTTTGCCGTTT‑TTGCCGCCCGCGCTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCCGCCgg                                                  >  2:305163/1‑149 (MQ=255)                                                                                                       aaTAGCGGGCCCATTGGCGGTTTTTGAAGATCATACCGGTCAGGGGAATCCACCGCGACCAGATAAAGCCCTTCCGCTTTTTCCGTTG‑GTGCCGCCCGCGGTCTGTCCTTTGCCGCCCGCAACTCCGCTATCCAGCTCTCCGGCTGGgtgt                                              >  1:251718/1‑151 (MQ=255)                                                                                                          aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTggg                                           <  1:152527/151‑1 (MQ=255)                                                                                                          aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                            <  1:194459/150‑1 (MQ=255)                                                                                                          aGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTgg                                            <  2:724953/150‑1 (MQ=255)                                                                                                                          gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATc                             >  1:558407/1‑149 (MQ=255)                                                                                                                          gTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAg                           <  2:817613/151‑1 (MQ=255)                                                                                                                            tttGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAgg                          <  2:1084954/150‑1 (MQ=255)                                                                                                                                     tCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCATCTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACa                 <  2:188233/150‑1 (MQ=255)                                                                                                                                      cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               >  2:232993/1‑151 (MQ=255)                                                                                                                                      cATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAat               <  2:148392/151‑1 (MQ=255)                                                                                                                                         aCCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAatat             <  1:190658/150‑1 (MQ=255)                                                                                                                                          ccGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTc           <  2:673050/151‑1 (MQ=255)                                                                                                                                           cGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTcc          <  1:878834/151‑1 (MQ=255)                                                                                                                                            ggTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCt         <  2:763881/151‑1 (MQ=255)                                                                                                                                             gTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTg        <  2:230493/151‑1 (MQ=255)                                                                                                                                                aGGGTAATCCACCGCGACCAGATACAGCCCGTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGAcc     <  1:673977/151‑1 (MQ=255)                                                                                                                                                 gggTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCa    <  2:968027/151‑1 (MQ=255)                                                                                                                                                   gTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCata  <  2:439359/151‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                       GGTCCATATTTGCCTATGCCTTGTACTCGTCATCTTTCAGGCTGTAACTACGTTGGCTACAGCCTGAAATGTTCCGAGCATTATTCTCGTTAGTCCGCCAGAAATAGCGGGCCCATTGGCGGTTTTGGAAGATCATACCGGTCAGGGTAATCCACCGCGACCAGATACAGCCCTTCCGCTTTTGCCGTTG‑CTGCCGCCAGCGTTCTGTCCTTTGCCGCCAGCAACTCTGCTATCCAGCTCTCCGGCTGGTTGTGGGCACCGACTTCCATCAGGCTGCCGACAATATTCCTGACCATA  >  NC_000913/2434735‑2435031

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