Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,898,162	T→G	Y47* (TAT→TAG)	cbrC →	UPF0167 family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,898,162	0	T	G	100.0%	227.0 / NA	62	Y47* (TAT→TAG)	cbrC	UPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base G (36/26);  total (36/26)

TTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGC  >  NC_000913/3898017‑3898299                                                                                                                                                  |                                                                                                                                          tttCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAGCAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTg                                                                                                                                         <  2:618620/148‑1 (MQ=255) tttCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGt                                                                                                                                      >  1:363812/1‑151 (MQ=255)   tCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAAACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGtt                                                                                                                                     <  1:265/150‑1 (MQ=255)     tATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGtt                                                                                                                                     <  1:783309/148‑1 (MQ=255)      aTGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTtgat                                                                                                                                  <  2:254517/150‑1 (MQ=255)             aaaaTATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTg                                                                                                                           >  1:483791/1‑150 (MQ=255)                 tatCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACAt                                                                                                                      >  2:256815/1‑151 (MQ=255)                      ggCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTg                                                                                                                 >  2:917323/1‑151 (MQ=255)                              ccccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCgtgg                                                                                                          >  1:375359/1‑150 (MQ=255)                               cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCgtggtg                                                                                                        >  1:443373/1‑151 (MQ=255)                               cccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCgtggt                                                                                                         >  1:775046/1‑150 (MQ=255)                                ccAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGgtgt                                                                                                       <  2:623017/151‑1 (MQ=255)                                    ttCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTg                                                                                                   >  1:30711/1‑151 (MQ=255)                                       aaaTATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCgg                                                                                                >  2:370048/1‑151 (MQ=255)                                        aaTATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGa                                                                                               <  2:299274/151‑1 (MQ=255)                                          tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACg                                                                                             >  2:100926/1‑151 (MQ=255)                                          tatCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACg                                                                                             >  2:913090/1‑151 (MQ=255)                                             catcCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGtt                                                                                          <  2:660253/151‑1 (MQ=255)                                                 ccAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTctgc                                                                                      >  1:895835/1‑151 (MQ=255)                                                  cAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTctgct                                                                                     <  1:735865/151‑1 (MQ=255)                                                   aaGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTctgctg                                                                                    >  2:835059/1‑151 (MQ=255)                                                        aCTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGaa                                                                               <  2:560976/151‑1 (MQ=255)                                                           ggAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGaaaa                                                                             <  1:660594/150‑1 (MQ=255)                                                             aaaCAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGtt                                                                          >  1:754664/1‑151 (MQ=255)                                                               aCAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTg                                                                        <  2:194971/151‑1 (MQ=255)                                                               aCAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTg                                                                        >  1:199469/1‑151 (MQ=255)                                                                      cATTTGAACAGGATAAAACCGTAGAGTTCGATTGCTGTGAACAACAGTCGTCAGTTTATTACTGGGGTCCCTTTTAGTGCGTTGCTGTAGTTGAACATCTCTGTCCGTGGTGTATTGCGTACGGTTCTGTTGCTGAAAAGTTTGCAg                                                                     >  1:42490/1‑147 (MQ=255)                                                                       aTTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGt                                                                >  2:8189/1‑151 (MQ=255)                                                                        tttGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGtt                                                               >  1:315069/1‑151 (MQ=255)                                                                         ttGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGttt                                                              >  2:194736/1‑151 (MQ=255)                                                                          tGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGtttt                                                             >  2:729470/1‑151 (MQ=255)                                                                           gAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTc                                                            <  1:831842/151‑1 (MQ=255)                                                                           gAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTc                                                            <  1:728495/151‑1 (MQ=255)                                                                            aaCAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCa                                                           <  2:104528/151‑1 (MQ=255)                                                                             aCAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAg                                                          <  1:223270/151‑1 (MQ=255)                                                                               aGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAgg                                                         >  2:553843/1‑150 (MQ=255)                                                                               aGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGga                                                        <  1:926008/151‑1 (MQ=255)                                                                                  aTAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGgatga                                                     <  1:594713/151‑1 (MQ=255)                                                                                     aaaCCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGc                                                  >  2:386679/1‑151 (MQ=255)                                                                                          gTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCa                                             <  1:605644/151‑1 (MQ=255)                                                                                            agagTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCAt                                            <  1:179355/150‑1 (MQ=255)                                                                                            agagTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATa                                           <  2:746745/151‑1 (MQ=255)                                                                                             gagTGCGATTGCTGTGAACAACAGACGTGAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCGGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCAt                                            <  1:417944/149‑1 (MQ=255)                                                                                                tGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGa                                         >  1:496076/1‑149 (MQ=255)                                                                                                   gATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGgtg                                    >  2:87736/1‑151 (MQ=255)                                                                                                       gCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg                                 >  2:169290/1‑150 (MQ=255)                                                                                                       gCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTg                                 >  1:338760/1‑150 (MQ=255)                                                                                                               caacaGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGt                        >  2:637622/1‑151 (MQ=255)                                                                                                                 acaGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTat                      <  2:477811/151‑1 (MQ=255)                                                                                                                 acaGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTat                      <  2:455535/151‑1 (MQ=255)                                                                                                                 acaGACGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTACAGGTAGTTTTCAGGATGATGTCAGCATAGAAGGTGTTGAATTTGAGTat                      <  2:425333/151‑1 (MQ=255)                                                                                                                     aCGTCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTatgatg                  <  2:922358/151‑1 (MQ=255)                                                                                                                      cGTTAGTTTATTTCTAGGGTCCCTTGTAGTGCGTTGCTGAAGTTGAAAATCTCTTTCCGGGGGGTATTGCGGACGGTTCTGCTGCTGACAAGTTTGCAGGTAGTTTTCAGAATGATGCCGGTATAGAAGGTGTTGAGTTTGAGGatga                    >  1:92006/1‑148 (MQ=255)                                                                                                                        tCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTAtgatga                 >  1:233159/1‑149 (MQ=255)                                                                                                                        tCAGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTAtgatga                 >  2:866767/1‑149 (MQ=255)                                                                                                                          aGTTTATTACTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAag               >  1:289677/1‑149 (MQ=255)                                                                                                                             ttattaCTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCTGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGAc          >  1:429304/1‑151 (MQ=255)                                                                                                                              tattaCTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACg         >  1:616601/1‑151 (MQ=255)                                                                                                                               attaCTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACg         >  2:609124/1‑150 (MQ=255)                                                                                                                                  aCTCGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAAtt     >  2:69584/1‑151 (MQ=255)                                                                                                                                     cGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAAttt    <  1:379025/149‑1 (MQ=255)                                                                                                                                     cGGGTCCCTTTTAGTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGc  >  2:305770/1‑151 (MQ=255)                                                                                                                                                  |                                                                                                                                          TTTCTATGACTCAAAATATCAGGCCGTTACCCCAATTCAAATATCATCCCAAGCCACTGGAAACAGGCGCATTTGAACAGGATAAAACCGTAGAGTGCGATTGCTGTGAACAACAGACGTCAGTTTATTACTCGGGTCCCTTTTATTGCGTTGATGAAGTTGAACATCTCTGTCCGTGGTGTATTGCGGACGGTTCTGCTGCTGAAAAGTTTGCAGGTAGTTTTCAGGATGATGCCAGCATAGAAGGTGTTGAATTTGAGTATGATGAAGAGGACGAATTTGC  >  NC_000913/3898017‑3898299

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