Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,801,372	C→T	7.1%	C400Y (TGT→TAT) ‡	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,372	0	C	T	7.1%	220.3 / 6.1	70	C400Y (TGT→TAT) ‡	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (30/35);  new base T (3/2);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

TTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAC  >  NC_000913/1801228‑1801510                                                                                                                                                 |                                                                                                                                           ttttCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATcc                                                                                                                                        >  1:5429/1‑149 (MQ=255)    tcttcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATGTATCCGtt                                                                                                                                     <  3:498650/149‑1 (MQ=255)      ttcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTaa                                                                                                                                   <  1:550654/149‑1 (MQ=255)       tcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAAc                                                                                                                                  >  1:82470/1‑149 (MQ=255)          aGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTc                                                                                                                               <  2:566331/149‑1 (MQ=255)           gcgcAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCa                                                                                                                              >  2:423928/1‑149 (MQ=255)             gcAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATc                                                                                                                            <  3:529630/149‑1 (MQ=255)               aaCCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATcgc                                                                                                                          >  1:561743/1‑149 (MQ=255)                aCCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATcgcg                                                                                                                         <  2:561743/149‑1 (MQ=255)                   gCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAt                                                                                                                      <  4:381809/149‑1 (MQ=255)                      aGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTg                                                                                                                   >  3:177080/1‑149 (MQ=255)                        gTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGtt                                                                                                                 <  4:185122/149‑1 (MQ=255)                           cGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGCGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGttctt                                                                                                              <  4:652153/149‑1 (MQ=255)                            gCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGttcttc                                                                                                             >  1:133013/1‑149 (MQ=255)                            gCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGttcttc                                                                                                             <  1:442579/149‑1 (MQ=255)                              cTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCCGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAg                                                                                                           >  2:284356/1‑149 (MQ=255)                                  gCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACTATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACa                                                                                                       <  3:265107/149‑1 (MQ=255)                                  gCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACa                                                                                                       >  2:484392/1‑149 (MQ=255)                                  gCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACa                                                                                                       <  2:11983/149‑1 (MQ=255)                                     cGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACTAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACagaa                                                                                                    <  1:360410/149‑1 (MQ=255)                                     cGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACagaa                                                                                                    <  1:45054/149‑1 (MQ=255)                                        tCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGat                                                                                                 <  1:528949/149‑1 (MQ=255)                                         cccACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGata                                                                                                >  1:505232/1‑149 (MQ=255)                                         cccACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGata                                                                                                <  2:293492/149‑1 (MQ=255)                                          ccACATTTCGTCGATGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGatat                                                                                               >  3:862/1‑149 (MQ=255)                                           cacaTTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATg                                                                                              >  2:164091/1‑149 (MQ=255)                                                    tcgCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATc                                                                                     <  2:133013/149‑1 (MQ=255)                                                      gCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCt                                                                                   <  3:713786/149‑1 (MQ=255)                                                       cTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATGTATCCGttt                                                                                                                                    >  4:170174/1‑98 (MQ=255)                                                       cTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATGTATCCGttt                                                                                                                                    <  3:170174/99‑2 (MQ=255)                                                        tGCCAATACGTTTTTCAGGACGCGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTgg                                                                                 <  2:221380/149‑1 (MQ=255)                                                        tGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTgg                                                                                 <  3:135711/149‑1 (MQ=255)                                                            aaTACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTaa                                                                             >  1:194454/1‑149 (MQ=255)                                                             aTACGTTTTTCAGGACGAGGGGAGATTTTGACAACGAACTTCTCGGAGCCAAAAGTGCCGTACAGATCATAGACTAAACGGATACATCCGTTATCTTCATCGCGAATTTGTTCTACAGTACAGAAGATATGCGCGTCACCCTGGGTaaa                                                                            >  4:649952/1‑149 (MQ=255)                                                                 gTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCca                                                                        <  1:284356/149‑1 (MQ=255)                                                                  tttttCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCcac                                                                       <  2:5429/149‑1 (MQ=255)                                                                        aGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCacgcacgc                                                                 >  4:663220/1‑149 (MQ=255)                                                                        aGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCacgcacgc                                                                 >  1:412303/1‑149 (MQ=255)                                                                             gAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATc                                                            >  2:736650/1‑149 (MQ=255)                                                                                   agagTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCTAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCa                                                      <  2:320169/149‑1 (MQ=255)                                                                                             gacgaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATcc                                                                         >  3:544098/1‑120 (MQ=255)                                                                                             gacgaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATcc                                                                         <  4:544098/120‑1 (MQ=255)                                                                                             gacgaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAAc                                            >  1:656320/1‑149 (MQ=255)                                                                                                gaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATc                                                                                                                            <  1:646929/66‑1 (MQ=255)                                                                                                gaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATc                                                                                                                            >  2:646929/1‑66 (MQ=255)                                                                                                 aTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGa                                        <  4:862/149‑1 (MQ=255)                                                                                                 aTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGa                                        <  1:423928/149‑1 (MQ=255)                                                                                                    ttctCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACgg                                     >  1:583551/1‑149 (MQ=255)                                                                                                       tcGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTc                                  >  2:146414/1‑149 (MQ=255)                                                                                                            gCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCATCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTAc                             >  2:176464/1‑149 (MQ=255)                                                                                                               aaaaGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAttt                                                                                                                    >  1:433170/1‑59 (MQ=255)                                                                                                               aaaaGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAAttt                                                                                                                    <  2:433170/59‑1 (MQ=255)                                                                                                                aaaGTTCTCTGCATATCATAGACGAAACGGATACATCCGATAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCGTCCTCGGCAAATCCACCCACGTGCATCATGCCATGCACCGAACCTGCCGGCGCGTTattgtg                         >  2:240190/1‑144 (MQ=255)                                                                                                                  aGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGTAGCGAACCTGACGGCTCGTTACGGTGGc                       <  3:114722/149‑1 (MQ=255)                                                                                                                   gTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCa                      <  2:265028/149‑1 (MQ=255)                                                                                                                        gTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTAc                 <  2:505232/149‑1 (MQ=255)                                                                                                                           cATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCaa              >  1:419883/1‑149 (MQ=255)                                                                                                                            atatCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCaaa             <  2:588822/149‑1 (MQ=255)                                                                                                                            atatCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCaaa             <  4:432875/149‑1 (MQ=255)                                                                                                                                cATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCg         >  2:648025/1‑149 (MQ=255)                                                                                                                                 aTAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACAGTGGCAGCTACCAAACTCgg        >  2:91173/1‑149 (MQ=255)                                                                                                                                  tAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGc       <  3:221965/149‑1 (MQ=255)                                                                                                                                    gACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAACTCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCa     <  1:176464/149‑1 (MQ=255)                                                                                                                                    gACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCa     <  3:233622/149‑1 (MQ=255)                                                                                                                                      gttaaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATa   >  2:472299/4‑149 (MQ=255)                                                                                                                                       ttaaCGGATGTATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAc  >  4:598698/3‑149 (MQ=255)                                                                                                                                       tAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAc  >  1:58091/1‑149 (MQ=255)                                                                                                                                       tAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAc  >  4:721703/1‑149 (MQ=255)                                                                                                                                          aCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGt                          >  4:398304/1‑122 (MQ=255)                                                                                                                                          aCGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGt                          <  3:398304/122‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                           TTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAC  >  NC_000913/1801228‑1801510

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