Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,800,755	C→T	55.4%	V606I (GTT→ATT)	thrS ←	threonyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,800,755	0	C	T	55.4%	27.1 / 112.3	94	V606I (GTT→ATT)	thrS	threonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (24/18);  new base T (26/26);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 5.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTTCTCA  >  NC_000913/1800656‑1800853                                                                                                    |                                                                                                   gTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAAc                                                                                                    <  1:3903344/100‑1 (MQ=255)  tttAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACt                                                                                                   >  1:3955772/1‑100 (MQ=255)  tttAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACt                                                                                                   <  1:797374/100‑1 (MQ=255)  tttAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACt                                                                                                   <  1:892037/100‑1 (MQ=255)    tAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACttt                                                                                                 >  1:1970704/1‑100 (MQ=255)     aaGACTGCGGCTGCGAATCTCTTGTTGCGGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTg                                                                                                <  2:896235/100‑1 (MQ=255)           gCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCCTTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTgct                                                                                          >  2:324018/1‑98 (MQ=255)           gCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGAt                                                                                          <  1:1678617/100‑1 (MQ=255)            cGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGagt                                                                                         <  1:3838266/100‑3 (MQ=255)             ggCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTc                                                                                        <  2:227516/100‑1 (MQ=255)              gCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTcc                                                                                       >  1:3262198/1‑100 (MQ=255)              gCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTcc                                                                                       <  1:747447/100‑1 (MQ=255)              gCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTcc                                                                                       <  1:2932113/100‑1 (MQ=255)              gCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTcc                                                                                       <  2:10688/100‑1 (MQ=255)                tGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCAc                                                                                     >  1:789441/1‑100 (MQ=255)                tGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCAc                                                                                     >  1:2936276/1‑100 (MQ=255)                       ttctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTa                                                                              <  2:974334/99‑1 (MQ=255)                       ctctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTa                                                                              >  2:2758115/1‑100 (MQ=255)                       ctctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTa                                                                              <  2:2076953/100‑1 (MQ=255)                        tctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTaa                                                                             >  1:974340/1‑99 (MQ=255)                         ctTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATc                                                                            <  1:1383487/100‑1 (MQ=255)                           tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcac                                                                          >  2:3581619/1‑100 (MQ=255)                           tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATcac                                                                          <  1:83939/100‑1 (MQ=255)                           tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATcac                                                                          >  1:3508428/1‑100 (MQ=255)                           tgttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATcac                                                                          <  1:1328640/100‑1 (MQ=255)                            gttgCAGCTTCTTGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcacc                                                                         >  1:2726505/1‑100 (MQ=255)                            gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcacc                                                                         >  1:3110405/1‑100 (MQ=255)                            gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcacc                                                                         >  2:3987470/1‑100 (MQ=255)                            gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcacc                                                                         <  2:441417/100‑1 (MQ=255)                            gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTTCCTGATTCCACCTCTTTATcacc                                                                         >  2:1346333/1‑100 (MQ=255)                            gttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATcacc                                                                         >  2:493407/1‑100 (MQ=255)                             ttgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcacc                                                                         >  1:2647460/1‑99 (MQ=255)                              tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcaccac                                                                       <  1:1204524/100‑1 (MQ=255)                              tgCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATcaccac                                                                       >  1:2422952/1‑100 (MQ=255)                               gCAGCTTCTCGATCACTTCATTTCCGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACcaca                                                                      <  2:272663/100‑1 (MQ=255)                               gCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACcaca                                                                      <  1:2345558/100‑1 (MQ=255)                               gCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACcaca                                                                      <  1:2416311/100‑1 (MQ=255)                                  gattCTCGCTCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGGGCGAACGGCAACTTTGCCTTATTCCACCTCTTGATCACCACAGAc                                                                   >  2:2830160/3‑100 (MQ=255)                                  gCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGAc                                                                   >  2:52397/1‑100 (MQ=255)                                  gCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGAc                                                                   >  2:2150836/1‑100 (MQ=255)                                   cTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGAcc                                                                  <  2:844713/100‑1 (MQ=255)                                   cTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGAcc                                                                  <  1:3901548/100‑1 (MQ=255)                                     tctcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCCCCACAGACCAg                                                                >  1:956447/1‑100 (MQ=255)                                     tctcGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAg                                                                <  1:1690938/100‑1 (MQ=255)                                        cGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCat                                                             >  1:2658687/1‑100 (MQ=255)                                            cACTTCAGTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCatata                                                          <  2:325233/99‑1 (MQ=255)                                             aCTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATg                                                        <  2:1809732/100‑1 (MQ=255)                                              cTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATgg                                                       >  1:1230718/1‑100 (MQ=255)                                              cTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATgg                                                       >  2:1418694/1‑100 (MQ=255)                                               ttCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATggg                                                      <  1:92273/100‑1 (MQ=255)                                                 cATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgac                                                    >  2:223613/1‑100 (MQ=255)                                                 cATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgac                                                    >  2:3274546/1‑100 (MQ=255)                                                 cATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgac                                                    >  1:3247776/1‑100 (MQ=255)                                                   tttACGTCCATGCTTCCCAGGTCTTTACCACGGCGTGTGCGAACGGCCCTTTTGCCTGATTCCACCTCTTTATCACCACAGACAAGCAGATATGGgacgc                                                  >  2:1828966/1‑100 (MQ=255)                                                     tACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcga                                                <  2:1204064/100‑1 (MQ=255)                                                      aCGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgac                                               <  1:929264/100‑1 (MQ=255)                                                       cgcccATGCTTCCTAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCATGATTCCACCTCTTTATCACCACAGCCCAGCATATATGGgacgcgan                                               <  2:3027386/96‑2 (MQ=255)                                                       cgatcttGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacg                                              <  1:1233159/94‑1 (MQ=255)                                                       cGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacg                                              >  1:1278415/1‑100 (MQ=255)                                                        gTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgccgc                                             >  2:3329862/1‑100 (MQ=255)                                                        gTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacgc                                             >  2:419960/1‑100 (MQ=255)                                                         tCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCCACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCa                                            >  2:1316889/1‑100 (MQ=255)                                                          ccATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCaa                                           >  2:3313658/1‑100 (MQ=255)                                                           cATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCaaa                                          <  1:659499/100‑1 (MQ=255)                                                           cATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCaaa                                          >  1:1760281/1‑100 (MQ=255)                                                            aTGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAg                                         <  2:1036270/100‑1 (MQ=255)                                                             tGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgt                                        >  1:3459584/1‑100 (MQ=255)                                                             tGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgt                                        <  2:3226792/100‑1 (MQ=255)                                                             tGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGgacgcgacgcgatcg                                        >  2:1233151/1‑95 (MQ=255)                                                               cTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAgtgt                                      >  1:3855157/1‑100 (MQ=255)                                                                 tCCCAGTTCTTTACCACGGCGGGTGCGAACGGAAACTTTGCCTGATTCCACCTCTTTATCACCGCAGACCAGCATATATGGGACGCGACGCAAAGTGTGc                                    <  2:749029/100‑1 (MQ=255)                                                                  cccAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCt                                   >  1:2737637/1‑100 (MQ=255)                                                                   ccAGGTCTTTACCACGGCGGGTGCGAACGGCAATATTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTc                                  >  2:1867419/1‑100 (MQ=255)                                                                   ccAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGCCGCAAAgtgt                                      >  2:1358875/1‑96 (MQ=255)                                                                    cAGGTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcg                                 >  1:307935/1‑100 (MQ=255)                                                                      ggTCTTTACCACGGCGGGTGCGAACGGCGATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                               <  2:3914623/100‑1 (MQ=255)                                                                      ggTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                               <  1:71304/100‑1 (MQ=255)                                                                      ggTCTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                               <  2:961704/100‑1 (MQ=255)                                                                      ggTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGCCGCAAAGTGTGCTcgcg                               >  2:913082/1‑100 (MQ=255)                                                                      ggTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTcgcg                               <  1:1080514/100‑1 (MQ=255)                                                                         cTTTACCACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAt                            >  1:1240226/1‑100 (MQ=255)                                                                         cTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAt                            >  1:2415510/1‑100 (MQ=255)                                                                           ttACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGAttt                          <  1:1205955/100‑1 (MQ=255)                                                                              ccACGGCGGGTGCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTaa                       >  2:267276/1‑100 (MQ=255)                                                                                     gggTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAt                >  2:1093012/1‑100 (MQ=255)                                                                                        tGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtctn             <  1:882215/100‑2 (MQ=255)                                                                                         gCGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttc            <  1:520925/100‑1 (MQ=255)                                                                                          cGAACGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAAtcttct           <  2:3855132/100‑1 (MQ=255)                                                                                            aaCGGCAATTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCa         <  1:2984071/100‑1 (MQ=255)                                                                                               ggCAACTTTGCCTGATTCCACCTCTTTATCACCCCAAACCAGCATATATGGGACGCGACGCAAAGTGAGCACGCGAATTTTAAAGCCAATCTTCTCAttt      >  2:3080881/1‑100 (MQ=255)                                                                                               ggCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCAttt      >  2:2646520/1‑100 (MQ=255)                                                                                                gCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTtc     >  2:2045018/1‑100 (MQ=255)                                                                                                gCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTtc     <  1:3164779/100‑1 (MQ=255)                                                                                                gCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATAGATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTtc     >  1:3680887/1‑100 (MQ=255)                                                                                                   aCTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTTCTCa  <  1:315186/100‑1 (MQ=255)                                                                                                    |                                                                                                   GTTTAAGACTGCGGCTGCGAATCTCTTGTTGCAGCTTCTCGATCACTTCATTTACGTCCATGCTTCCCAGGTCTTTACCACGGCGGGTGCGAACGGCAACTTTGCCTGATTCCACCTCTTTATCACCACAGACCAGCATATATGGGACGCGACGCAAAGTGTGCTCGCGGATTTTAAAGCCAATCTTCTCATTTCTCA  >  NC_000913/1800656‑1800853

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