Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,182,928	A→T	E562V (GAA→GTA)	rpoB →	RNA polymerase, beta subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,182,928	0	A	T	100.0%	213.5 / NA	68	E562V (GAA→GTA)	rpoB	RNA polymerase, beta subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base T (37/31);  total (37/31)

GTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTAT  >  NC_000913/4182782‑4183069                                                                                                                                                   |                                                                                                                                              gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGtaa                                                                                                                                             >  2:16839/1‑149 (MQ=255) gTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCCCCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGtaa                                                                                                                                             >  2:14534/1‑149 (MQ=255)  tttATGGACCAGAACAACCCGCTGTCTGAGGTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAACCGTAAc                                                                                                                                            >  1:809480/1‑149 (MQ=255)  tttATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAAc                                                                                                                                            <  1:596943/149‑1 (MQ=255)        gACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                      <  1:37391/149‑1 (MQ=255)        gACCAGAACAACCCGCAGTCCAAGTGCACGCACAAACGTCCGATCTCCGCCCTCGGCCCAGGCGGTCTGACCCGTGAATGTGCATGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGa                                                                                                                                      <  2:169216/149‑1 (MQ=255)         aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGCCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGac                                                                                                                                     >  2:711186/1‑148 (MQ=255)         aCCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGaa                                                                                                                                     >  1:618462/1‑149 (MQ=255)          ccAGAACAACCCGCTGTCCGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAg                                                                                                                                    <  1:787275/149‑1 (MQ=255)           cAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACCACGGTCGCGTATGTCCAATCGTAACCCCTGAAgg                                                                                                                                   >  2:676721/1‑149 (MQ=255)              aacaacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTcc                                                                                                                                <  1:329735/149‑1 (MQ=255)                caacCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGa                                                                                                                              >  1:816231/1‑149 (MQ=255)                   cccGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACa                                                                                                                           >  1:271244/1‑149 (MQ=255)                            tGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGa                                                                                                                  >  2:56070/1‑149 (MQ=255)                             gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                                                                 >  1:295798/1‑149 (MQ=255)                             gagaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGAt                                                                                                                 >  1:479519/1‑149 (MQ=255)                              agaTTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGCCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATc                                                                                                                >  1:119439/1‑149 (MQ=255)                                  tACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCCCATCAAct                                                                                                            <  2:808836/149‑1 (MQ=255)                                     gCACAAACGTCGTATCTCCGCACTCGGCACAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGGAGGTGCGCACATGGGTATGATCAActctc                                                                                                         >  1:206236/1‑149 (MQ=255)                                        caAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGt                                                                                                      <  2:68508/149‑1 (MQ=255)                                              tcgtATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACCTCGGTCTGATCAACTCTCTGTCCgtgt                                                                                                >  1:185327/1‑149 (MQ=255)                                                   tctcGGCACTCGGCCCAGGCGGTCTGCCCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGca                                                                                           <  2:370769/149‑1 (MQ=255)                                                     tcCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGcaca                                                                                         >  2:306118/1‑149 (MQ=255)                                                      ccccACTCTGCCCAGGCGGTCTGTCCCGTTAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTACCCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAg                                                                                        <  2:754781/146‑1 (MQ=255)                                                        gCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTCCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGAc                                                                                      <  1:303213/149‑1 (MQ=255)                                                            tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                  >  1:795563/1‑149 (MQ=255)                                                            tCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAAc                                                                                  >  2:679308/1‑149 (MQ=255)                                                              ggCCCAGGCGGTCTGACCCGTGAACGTGCTGGCTTCGAAGTTCGAGACGTACAACCGACTCGCTACGATCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGAGCAACTCTATGTCCGTGTAAGCACAGACTAACGa                                                                                >  1:499957/1‑149 (MQ=255)                                                                 ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGACCATCGGTCTGATCAACTGCCTTCCCGTGTACGCCCAGCCTAAAAAATa                                                                             >  1:4170/1‑149 (MQ=255)                                                                 ccAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATa                                                                             >  1:166941/1‑149 (MQ=255)                                                                      cGGTCTGCCCCGTGCACGTGCCGGCTTCGGAGTTCGAGACGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCTGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCt                                                                        <  2:526478/149‑1 (MQ=255)                                                                       gttcggACCCGTGCACGCGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                                                       <  2:809480/144‑1 (MQ=255)                                                                       ggTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCtt                                                                       <  1:726636/149‑1 (MQ=255)                                                                        gTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACCGACTCACGAATACGGCTTc                                                                      >  2:141913/1‑149 (MQ=255)                                                                         tCTGACCCGAGAACGTGAAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACCCTCTGTCCGTGTACGCACAGACTAACGAAGACGGCTTcc                                                                     >  2:304561/1‑149 (MQ=255)                                                                           tggcccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGGCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                   <  1:39237/146‑1 (MQ=255)                                                                           tGACTCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCtt                                                                   >  2:709963/1‑149 (MQ=255)                                                                               ccGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACACTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTgaga                                                               >  2:490330/1‑149 (MQ=255)                                                                                  tGAACGTGCAGGCTGCGAAGTCCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTc                                                            <  2:327400/149‑1 (MQ=255)                                                                                   gaccGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCACCTTTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTcc                                                           <  2:615960/146‑1 (MQ=255)                                                                                     aCGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAAGACGGCTTCCTTGAGACTCCGt                                                         >  1:672609/1‑149 (MQ=255)                                                                                          cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCATTCGTAACCCCTGAAGGTCCGAACATCGGt                                                                                                                      <  2:132006/83‑1 (MQ=255)                                                                                          cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGCAGGTCCGAACATCGGt                                                                                                                      >  1:132006/1‑83 (MQ=255)                                                                                          cAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGt                                                    >  2:784432/1‑149 (MQ=255)                                                                                             gCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTaaa                                                 >  2:182893/1‑149 (MQ=255)                                                                                              cctCGAAGTTCGAGACTTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAg                                                <  1:4291/147‑1 (MQ=255)                                                                                               ttCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATATGGCTTCCTTGAGACTCCGTATCGTAAAGt                                               <  1:501451/149‑1 (MQ=255)                                                                                                     gTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTCACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGa                                         >  1:448939/1‑149 (MQ=255)                                                                                                           gaCGTACACCCGACTCACGCCGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTCTACGCCCCGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGgtgt                                   <  1:564034/149‑1 (MQ=255)                                                                                                            aCGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCTGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGtgtt                                  <  1:553175/149‑1 (MQ=255)                                                                                                             cGTACCCCCGACTCACTACGGTCGCGTATGTCCAATCGTCACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCCGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGtgttg                                 <  2:58107/149‑1 (MQ=255)                                                                                                                acacCCGTCTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTaa                              <  1:291064/149‑1 (MQ=255)                                                                                                                   cccGACTCACTACGAACGCGTATGACGAATCGTAACACCTGATTGTCGGAAGAGCGGTCAGATGGACGCTCGGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTg                           >  2:602121/1‑149 (MQ=255)                                                                                                                    ccGGCTCAGTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACGCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                          <  1:606224/149‑1 (MQ=255)                                                                                                                    ccGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGGGACTCCGTATCGTCACGTGACCGCCGGTGTTGTAACTGa                          >  1:789347/1‑149 (MQ=255)                                                                                                                    ccGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                          >  1:786005/1‑149 (MQ=255)                                                                                                                    ccGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGCCTGCTCAACTCTCTTTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGa                          >  1:360986/1‑149 (MQ=255)                                                                                                                     cGACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGAc                         <  2:185327/149‑1 (MQ=255)                                                                                                                      gACTCACTACGGTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACg                        >  1:693127/1‑149 (MQ=255)                                                                                                                                ggTCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGCATACGGCTTCCTTGAGACTCCGTATCGTAAAGTAACCGACGGTGTTGTAACTGACGAAATTCCCTa              >  1:634551/1‑149 (MQ=255)                                                                                                                                ggTCGCGGATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTa              <  1:16839/149‑1 (MQ=255)                                                                                                                                  tCGCGTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTAcc            >  2:192652/1‑149 (MQ=255)                                                                                                                                     cgTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGt         <  1:657863/149‑1 (MQ=255)                                                                                                                                     cgTATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGt         <  1:711186/149‑1 (MQ=255)                                                                                                                                       tATGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCt       <  2:295798/149‑1 (MQ=255)                                                                                                                                       tATGTCCAAGCGTCACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCt       <  1:619157/149‑1 (MQ=255)                                                                                                                                        aTGTCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTg      <  1:494611/149‑1 (MQ=255)                                                                                                                                           tCCAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACCCCGCATCGTAAAGTGACCGACGGTGTTGTAACTGACAAAATGCCCTACCTGTCTgatg   >  1:249119/1‑146 (MQ=255)                                                                                                                                           tCCAATCGTAACCCCTGAAGGTCCGAACATCGATCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGCCGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTa   >  2:291518/1‑149 (MQ=255)                                                                                                                                            ccAATCGTAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTAt  >  1:342018/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                              GTTTATGGACCAGAACAACCCGCTGTCTGAGATTACGCACAAACGTCGTATCTCCGCACTCGGCCCAGGCGGTCTGACCCGTGAACGTGCAGGCTTCGAAGTTCGAGACGTACACCCGACTCACTACGGTCGCGTATGTCCAATCGAAACCCCTGAAGGTCCGAACATCGGTCTGATCAACTCTCTGTCCGTGTACGCACAGACTAACGAATACGGCTTCCTTGAGACTCCGTATCGTAAAGTGACCGACGGTGTTGTAACTGACGAAATTCACTACCTGTCTGCTAT  >  NC_000913/4182782‑4183069

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