Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,801,371	A→G	7.1%	C400C (TGT→TGC) ‡	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,801,371	0	A	G	7.1%	218.1 / 6.2	70	C400C (TGT→TGC) ‡	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (30/35);  new base G (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 = 8.69e-01

TGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAC  >  NC_000913/1801223‑1801510                                                                                                                                                     |                                                                                                                                            tgttgtTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATa                                                                                                                                             >  4:378902/1‑149 (MQ=255)      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                                                                                                                            >  2:646929/1‑66 (MQ=255)                                                                                                     gaTCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATc                                                                                                                            <  1:646929/66‑1 (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)                                                                                                                                                     |                                                                                                                                            TGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCCCACATTTCGTCGCTGCCAATACGTTTTTCAGGACGAGTGGAGAGTTTGACGACGATCTTCTCGAAGCCAAAAGTGCTGTACATATCATAGACTAAACGGATACATCCGTTAACTTCATCGCGAATTTGTTCTTCAGTACAGAAGATATGCGCGTCATCCTGGGTAAATCCACGCACGCGCATCAGGCCATGCAGCGAACCTGACGGCTCGTTACGGTGGCAGCTACCAAACTCGGCCATAC  >  NC_000913/1801223‑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: ‑