Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,049,975	C→A	A231A (GCG→GCT)	gcvT ←	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,049,975	0	C	A	98.5%	233.6 / ‑2.8	67	A231A (GCG→GCT)	gcvT	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (33/33);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

CCATGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGG  >  NC_000913/3049838‑3050111                                                                                                                                          |                                                                                                                                         ccATGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGc                                                                                                                                        >  1:45884/1‑140 (MQ=255)    tGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTc                                                                                                                                     >  1:1000692/1‑140 (MQ=255)    tGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTc                                                                                                                                     >  1:70809/1‑140 (MQ=255)    tGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTc                                                                                                                                     >  1:638636/1‑140 (MQ=255)     gCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTcc                                                                                                                                    >  1:555157/1‑140 (MQ=255)        cACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGa                                                                                                                                 <  1:1128004/140‑1 (MQ=255)        cACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGa                                                                                                                                 <  2:855274/140‑1 (MQ=255)             tGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAg                                                                                                                            <  2:463878/140‑1 (MQ=255)               cACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCg                                                                                                                          >  1:1223329/1‑140 (MQ=255)                    ccAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCa                                                                                                                     >  2:1170207/1‑140 (MQ=255)                         gCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCAcgcgc                                                                                                                <  1:352752/140‑1 (MQ=255)                            tCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGTTGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGcc                                                                                                             >  2:265791/1‑140 (MQ=255)                            tCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGcc                                                                                                             >  2:757586/1‑140 (MQ=255)                            tCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGcc                                                                                                             >  1:343450/1‑140 (MQ=255)                             cACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGccc                                                                                                            <  1:179311/140‑1 (MQ=255)                              aCGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACgcgcgcgca                                                                                                           <  2:70809/140‑4 (MQ=255)                                  ccGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGcc                                                                                                       >  1:13653/1‑140 (MQ=255)                                   cGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCca                                                                                                      >  2:276190/1‑140 (MQ=255)                                     aTAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcaca                                                                                                    <  2:5513/140‑1 (MQ=255)                                     aTAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcaca                                                                                                    >  1:861423/1‑140 (MQ=255)                                     aTAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCcaca                                                                                                    <  2:570885/140‑1 (MQ=255)                                      tAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACAt                                                                                                   >  1:1119948/1‑140 (MQ=255)                                        aaGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATgg                                                                                                 <  2:555157/140‑1 (MQ=255)                                             aCGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTaa                                                                                            >  2:770994/1‑140 (MQ=255)                                                aTCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAacac                                                                                         <  2:643228/140‑1 (MQ=255)                                                aTCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAacac                                                                                         <  2:1035688/140‑1 (MQ=255)                                                aTCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAacac                                                                                         <  1:746768/140‑1 (MQ=255)                                                aTCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAacac                                                                                         <  2:557207/140‑1 (MQ=255)                                                             cAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAg                                                                            >  2:1258633/1‑140 (MQ=255)                                                               ggCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAgcg                                                                          <  2:799346/140‑1 (MQ=255)                                                                   aTGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACg                                                                      >  1:512004/1‑140 (MQ=255)                                                                     ggTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGcc                                                                    >  1:733715/1‑140 (MQ=255)                                                                        ccAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGa                                                                 >  1:1200013/1‑140 (MQ=255)                                                                               aTGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGc                                                          <  1:770994/140‑1 (MQ=255)                                                                                tGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGcc                                                         >  2:165313/1‑140 (MQ=255)                                                                                 gTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCg                                                        >  1:1000885/1‑140 (MQ=255)                                                                                    ggcggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtt                                                     <  1:731009/140‑1 (MQ=255)                                                                                    ggcggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtt                                                     <  2:512004/140‑1 (MQ=255)                                                                                    ggcggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtt                                                     <  1:168878/140‑1 (MQ=255)                                                                                    ggcggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtt                                                     <  2:13653/140‑1 (MQ=255)                                                                                      cggcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAATACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCtttt                                                   >  1:1130239/1‑140 (MQ=255)                                                                                        gcTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTc                                                 >  2:996259/1‑140 (MQ=255)                                                                                              ggAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTggg                                           <  1:1253147/140‑1 (MQ=255)                                                                                                     tGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCa                                    <  2:1200013/140‑1 (MQ=255)                                                                                                     tGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCa                                    <  2:736910/140‑1 (MQ=255)                                                                                                         ttCGTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAAttt                                >  2:716618/1‑140 (MQ=255)                                                                                                            gTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCAt                             >  2:526991/1‑140 (MQ=255)                                                                                                            gTCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCAt                             >  2:575805/1‑140 (MQ=255)                                                                                                             tCCATCTCCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGACACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATa                            <  2:854580/140‑1 (MQ=255)                                                                                                                 tctcCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGccc                        >  2:170858/1‑140 (MQ=255)                                                                                                                  ctcCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAgcg                                      >  2:943266/1‑125 (MQ=255)                                                                                                                  ctcCTGACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAgcg                                      <  1:943266/125‑1 (MQ=255)                                                                                                                     ttGACCATAAAGATTCATGCCAGCTTCCTGACGCAGCGTGCCACGCGCGCACAAGCCACATGCCTTAACACCCGCTTCCCCCAGCGCACGCCAGAAATCGGCCGTTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCtt                    <  2:861423/139‑1 (MQ=255)                                                                                                                       gACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCa                  <  2:917769/140‑1 (MQ=255)                                                                                                                       gACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCa                  <  2:1158450/140‑1 (MQ=255)                                                                                                                       gACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCa                  <  2:1119948/140‑1 (MQ=255)                                                                                                                       gACCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCa                  <  1:1170207/140‑1 (MQ=255)                                                                                                                        aCCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCAc                 >  1:1065426/1‑140 (MQ=255)                                                                                                                        aCCATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCAc                 <  2:343450/140‑1 (MQ=255)                                                                                                                          cATAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCg               >  1:649609/1‑140 (MQ=255)                                                                                                                            tAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGt             <  1:265791/140‑1 (MQ=255)                                                                                                                            tAAAGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGt             <  1:663180/140‑1 (MQ=255)                                                                                                                             aaaGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGta            >  2:758423/1‑140 (MQ=255)                                                                                                                             aaaGATTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGta            >  1:975687/1‑140 (MQ=255)                                                                                                                                 aTTCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAATACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAAc        <  2:1130239/140‑1 (MQ=255)                                                                                                                                  ttCATGCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAAcc       >  2:27702/1‑140 (MQ=255)                                                                                                                                    cATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAg     >  1:654625/1‑140 (MQ=255)                                                                                                                                       gCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtgg  <  2:1000692/140‑1 (MQ=255)                                                                                                                                       gCCAGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAgtgg  <  2:223134/140‑1 (MQ=255)                                                                                                                                          |                                                                                                                                         CCATGCTCACGCTGCACTTCCAGGGCTTCACGACCGATAAAGTCACGATCTGCCGGTTCCCAGGCGATGGTCCAGCCCATGTTGGCGGCTAAAGGAGAGATGGTTTCGTCCATCTCCTGACCATAAAGATTCATGCCCGCTTCCAGACGCAGCGTGTCACGCGCGCCCAAGCCACATGGCTTAACACCCGCTTCCACCAGCGCACGCCAGAAATCGGCCGCTTTTTCATTGGGCAGCGCAATTTCATAGCCCGCTTCACCGGTATAACCAGTGG  >  NC_000913/3049838‑3050111

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