Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,407,898	A→C	6.2%	Q113P (CAG→CCG) ‡	alaA →	glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,407,898	0	A	C	6.2%	227.8 / 6.9	65	Q113P (CAG→CCG) ‡	alaA	glutamate‑pyruvate aminotransferase; glutamic‑pyruvic transaminase (GPT); alanine transaminase
Reads supporting (aligned to +/- strand):  ref base A (32/29);  new base C (2/2);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.83e-01

TCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTC  >  NC_000913/2407755‑2408044                                                                                                                                                |                                                                                                                                                    tCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCAt                                                                                                                                                <  3:189536/149‑1 (MQ=255)     gggTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑Ct                                                                                                                                           >  1:402642/1‑149 (MQ=255)               aTTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCggg                                                                                                                                 <  1:118760/149‑1 (MQ=255)                  ccAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGAc                                                                                                                              <  2:421842/149‑1 (MQ=255)                   cAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACg                                                                                                                             >  1:346490/1‑149 (MQ=255)                    aaaGGTCTTTACTCCGCGCGTAAATCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATTCCTGCATTG‑cc                                                                                                                                           <  3:543679/134‑2 (MQ=255)                    aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGATCTTATCGTTCAGGCAATGCCTGCATTG‑cc                                                                                                                                           >  4:543679/1‑133 (MQ=255)                    aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                                                                            <  2:231039/133‑1 (MQ=255)                    aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑c                                                                                                                                            >  1:231039/1‑133 (MQ=255)                    aaaGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                            >  1:231938/1‑149 (MQ=255)                        gTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATg                                                                                                                        >  4:240565/1‑149 (MQ=255)                        gTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATg                                                                                                                        >  4:110278/1‑149 (MQ=255)                        gTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATg                                                                                                                        >  1:533082/1‑149 (MQ=255)                          cTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGtt                                                                                                                      <  1:54322/149‑1 (MQ=255)                           tttACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTg                                                                                                                     >  4:548510/1‑149 (MQ=255)                                tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                            >  3:110765/1‑137 (MQ=255)                                tCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGa                                                                                                                            <  4:110765/137‑1 (MQ=255)                                           aGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGtt                                                                                                                      >  4:530575/1‑132 (MQ=255)                                           aGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGtt                                                                                                                      <  3:530575/132‑1 (MQ=255)                                           aGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGAtt                                                                                                     >  1:298423/1‑149 (MQ=255)                                            gCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAAc                                                                                                                                       >  3:260015/1‑114 (MQ=255)                                            gCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAAc                                                                                                                                       <  4:260015/114‑1 (MQ=255)                                            gCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTa                                                                                                    >  2:3356/1‑149 (MQ=255)                                             ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                                                   >  1:144904/1‑149 (MQ=255)                                             ccATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAc                                                                                                   >  4:310890/1‑149 (MQ=255)                                              catcatGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTAcc                                                                                                  <  2:402642/149‑1 (MQ=255)                                                    gcagcaCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCActct                                                                                            >  1:179776/1‑149 (MQ=255)                                                        caCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgataa                                                                                                                                  <  3:438530/107‑6 (MQ=37)                                                        caCTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCA‑GCAATGCCTGCATTGCCTGAAcgataa                                                                                                                                  >  4:438530/1‑102 (MQ=37)                                                            aCCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACcgcg                                                                                    <  1:120505/149‑1 (MQ=255)                                                                  cTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGtt                                                                              <  3:548510/149‑1 (MQ=255)                                                                    cGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTc                                                                            <  2:533082/149‑1 (MQ=255)                                                                     gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                                           <  3:110278/149‑1 (MQ=255)                                                                     gTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCg                                                                           <  3:240565/149‑1 (MQ=255)                                                                      tGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGc                                                                          <  2:144904/149‑1 (MQ=255)                                                                               tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                 >  2:379801/1‑149 (MQ=255)                                                                               tGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCg                                                                 >  1:141261/1‑149 (MQ=255)                                                                                       ccGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCg                                                         <  4:19745/149‑1 (MQ=255)                                                                                        cGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCgg                                                        >  3:278920/1‑149 (MQ=255)                                                                                         gTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGt                                                       <  4:145468/149‑1 (MQ=255)                                                                                          tGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTg                                                      >  2:156702/1‑149 (MQ=255)                                                                                            gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa                                                    <  1:3356/149‑1 (MQ=255)                                                                                            gAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCa                                                    >  1:444535/1‑149 (MQ=255)                                                                                                   tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCtt                                             >  1:13419/1‑149 (MQ=255)                                                                                                   tttACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCtt                                             >  4:530049/1‑149 (MQ=255)                                                                                                      aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGTGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCTGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                          <  4:34441/149‑1 (MQ=255)                                                                                                      aCATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGc                                          >  4:518615/1‑149 (MQ=255)                                                                                                         tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                        >  1:282400/1‑148 (MQ=255)                                                                                                         tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGa                                        <  2:282400/148‑1 (MQ=255)                                                                                                         tCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGAt                                       <  2:13419/149‑1 (MQ=255)                                                                                                            gCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGaa                                    >  1:378441/1‑149 (MQ=255)                                                                                                            gCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGaa                                    <  4:449331/149‑1 (MQ=255)                                                                                                                tGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCct                                <  3:310890/149‑1 (MQ=255)                                                                                                                tGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCct                                <  1:93480/149‑1 (MQ=255)                                                                                                                tGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCct                                <  2:179776/149‑1 (MQ=255)                                                                                                                       tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                         >  1:293425/1‑149 (MQ=255)                                                                                                                       tCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTg                         >  2:123012/1‑149 (MQ=255)                                                                                                                         ggAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGt                       >  1:219388/1‑149 (MQ=255)                                                                                                                            gCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTcc                    <  2:346490/149‑1 (MQ=255)                                                                                                                                      cAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCga          <  2:378441/149‑1 (MQ=255)                                                                                                                                        ggCAATGCAGGCATTG‑CTGAACAGCGGGGATGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGttt                                                                             >  1:130712/1‑80 (MQ=255)                                                                                                                                        ggCAATGCAGGCATTG‑CTGAACAGCGGGGATGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGttt                                                                             <  2:130712/80‑1 (MQ=255)                                                                                                                                        ggCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatg        >  1:619392/1‑149 (MQ=255)                                                                                                                                         gCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCgatga       <  1:153997/149‑1 (MQ=255)                                                                                                                                              gCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTc  >  2:579648/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                    TCAAGGGTATTGCGATTCCAAAGGTCTTTACTCCGCGCGTAAAGCCATCATGCAGCACTACCAGGCTCGTGGCATGCGTGATGTTACCGTGGAAGATATTTACATCGGCAATGGTGTATCGGAGCTTATCGTTCAGGCAATGCAGGCATTG‑CTGAACAGCGGGGACGAAATGTTGGTTCCTGCACCAGATTACCCACTCTGGACCGCGGCGGTTTCGCTTTCCAGCGGTAAAGCGGTGCATTATCTTTGCGATGAATCCTCTGACTGGTTCCCGGACCTCGATGATATTC  >  NC_000913/2407755‑2408044

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