Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,062,255	T→C	V469A (GTA→GCA)	scpA →	methylmalonyl‑CoA mutase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,062,255	0	T	C	98.6%	238.6 / ‑3.4	69	coding (1406/2145 nt)	scpA	methylmalonyl‑CoA mutase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (32/36);  minor base . (0/1);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

AAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAAC  >  NC_000913/3062119‑3062382                                                                                                                                         |                                                                                                                                aaaGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACt                                                                                                                              >  2:300663/1‑140 (MQ=255) aaaGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACt                                                                                                                              >  2:860456/1‑140 (MQ=255)    gCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTga                                                                                                                           <  2:399713/140‑1 (MQ=255)      gATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTgaga                                                                                                                         <  2:966122/140‑1 (MQ=255)      gATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTgaga                                                                                                                         <  1:1618642/140‑1 (MQ=255)       aTCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGAt                                                                                                                        <  2:373763/140‑1 (MQ=255)        tCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATc                                                                                                                       <  2:563341/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  2:225144/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  1:860456/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  2:1260232/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  1:300663/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  2:244921/140‑1 (MQ=255)          gAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  2:700220/140‑1 (MQ=255)          gAAGCAGGTCCGCCCAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGa                                                                                                                     <  2:1407254/140‑1 (MQ=255)            aGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACa                                                                                                                   <  2:383386/140‑1 (MQ=255)               aGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACg                                                                                                                <  2:711419/140‑1 (MQ=255)                 gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTg                                                                                                              <  2:1626043/140‑1 (MQ=255)                 gTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTg                                                                                                              <  1:1042118/140‑1 (MQ=255)                  tCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGa                                                                                                             <  1:1312348/140‑1 (MQ=255)                  tCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGa                                                                                                             <  1:698083/140‑1 (MQ=255)                  tCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGAACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGa                                                                                                             <  1:949248/140‑1 (MQ=255)                          aaCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGt                                                                                                     <  2:34670/140‑1 (MQ=255)                           aCGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTa                                                                                                    >  1:415608/1‑140 (MQ=255)                                        gagGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGc                                                                                       >  1:1088121/1‑140 (MQ=255)                                         agGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATG‑CCTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCtt                                                                                     <  1:1225037/140‑1 (MQ=255)                                           gCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTc                                                                                    >  2:274409/1‑140 (MQ=255)                                            ccTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCg                                                                                   <  1:287731/140‑1 (MQ=255)                                            ccTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCg                                                                                   >  1:396604/1‑140 (MQ=255)                                             cTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGc                                                                                  >  1:1297523/1‑140 (MQ=255)                                              tCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCt                                                                                 >  2:901227/1‑140 (MQ=255)                                                 gcgcgcGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGa                                                                              >  1:1075628/1‑140 (MQ=255)                                                     gcgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGc                                                                          >  2:1346048/1‑140 (MQ=255)                                                      cgAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCa                                                                         <  1:758417/140‑1 (MQ=255)                                                        aaCAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCAtt                                                                       <  1:341603/140‑1 (MQ=255)                                                            gTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTcgcg                                                                   <  1:901227/140‑1 (MQ=255)                                                             tCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCgcgc                                                                  <  2:396604/140‑1 (MQ=255)                                                             tCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCgcgc                                                                  <  1:1346048/140‑1 (MQ=255)                                                              cGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGcc                                                                 <  2:1088121/140‑1 (MQ=255)                                                                   aTCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCg                                                            >  2:393437/1‑140 (MQ=255)                                                                       aCCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCgcgc                                                                  >  2:1355419/1‑130 (MQ=255)                                                                       aCCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCgcgc                                                                  <  1:1355419/130‑1 (MQ=255)                                                                       aCCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGtgat                                                        >  2:1370201/1‑140 (MQ=255)                                                                          aGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTgatgat                                                     <  1:203118/140‑1 (MQ=255)                                                                           gggCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGatgatg                                                    >  2:1377499/1‑140 (MQ=255)                                                                             gCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATgcc                                                  >  1:887407/1‑140 (MQ=255)                                                                             gCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATgcc                                                  >  2:435266/1‑140 (MQ=255)                                                                                aGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATgccgcc                                               >  2:1329531/1‑140 (MQ=255)                                                                                  cGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGt                                             >  2:782276/1‑140 (MQ=255)                                                                                  cGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGt                                             >  2:435478/1‑140 (MQ=255)                                                                                      tCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAcc                                         >  2:97422/1‑140 (MQ=255)                                                                                      tCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAcc                                         <  1:1377499/140‑1 (MQ=255)                                                                                       cATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAAccg                                        <  1:782276/140‑1 (MQ=255)                                                                                           gTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACcgccgc                                    >  1:1299175/1‑140 (MQ=255)                                                                                                tgtCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGGGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGa                               >  1:1163284/1‑140 (MQ=255)                                                                                                 gtCAACAAGTACAAACTGGATCACGAAGACGAACCCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGaa                              <  2:1297523/140‑1 (MQ=255)                                                                                                   caacaaGTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCGGTGATGATGCCGCCGTAACCGCCGCGTTGAACg                            <  1:1129066/140‑1 (MQ=255)                                                                                                     acaaTTAAAAACTGGATCACGAAGACGAACCCGATGCTTTTGAGATCAACAACGTATTGGTGCGTAACGAGCAAATTGTTCCGCTGGACCGCATCCGCGCCGCCCGTGATAAGGCCGCCGTAACCGCCGCGTTGAACGcc                          <  2:813786/140‑1 (MQ=255)                                                                                                         gTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGa                      >  1:1525056/1‑140 (MQ=255)                                                                                                         gTACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGa                      >  2:630408/1‑140 (MQ=255)                                                                                                          tACAAACTGGGTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATCGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGAc                     >  2:20575/1‑140 (MQ=255)                                                                                                          tACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGAc                     >  2:867383/1‑140 (MQ=255)                                                                                                          tACAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGAc                     <  1:1455083/140‑1 (MQ=255)                                                                                                            cAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTc                   >  2:489435/1‑140 (MQ=255)                                                                                                            cAAACTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTc                   >  2:1440446/1‑140 (MQ=255)                                                                                                               ccTGGATCAGGAAGACAAAACCGATCCATTTGAGATCGAAAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                <  2:1299175/139‑1 (MQ=255)                                                                                                               aCTGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                <  1:274409/140‑1 (MQ=255)                                                                                                               aCTGAAACACAAAGACAAAACGGATGCATTTGAGATCGACAAGGTGATGGTGCGTAACGAGAAATTTGCTTCGCGGGACCGCTTTCGCGCCCCCCGTGATGAGGCCGCCGTAACCGCCGCGTTGAACGCCCCGACTCAcg                <  2:315839/140‑1 (MQ=255)                                                                                                                 tGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                >  2:434338/1‑138 (MQ=255)                                                                                                                 tGGATCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                <  1:434338/138‑1 (MQ=255)                                                                                                                  ggATCACGAAGACGAACCCGATGCACTGGAGATCGACACCGTGATGGTGCGTAACGAGAAAATTGCTTCGCTGGACCGCATCCGCCCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcgccg             <  2:1525056/140‑1 (MQ=255)                                                                                                                    aTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGca           <  2:1075628/140‑1 (MQ=255)                                                                                                                    aTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                >  1:857725/1‑135 (MQ=255)                                                                                                                    aTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCAcg                <  2:857725/135‑1 (MQ=255)                                                                                                                    aTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGca           <  1:867383/140‑1 (MQ=255)                                                                                                                    aTCACGAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGca           >  2:1265860/1‑140 (MQ=255)                                                                                                                     ccaggAAGACAAACCCGATGCATTTGAGATCGACAACGTGATGGTGCGTAACGAGAAAATTGCTTCGCTGGACCGCATTCGCGCCCCCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGcac          <  2:1163284/136‑1 (MQ=255)                                                                                                                         gAAGACGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCa      <  1:97422/140‑1 (MQ=255)                                                                                                                             aCGAAACCGATGCACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAAc  >  2:659886/1‑140 (MQ=255)                                                                                                                                         |                                                                                                                                AAAGCGATCGAAGCAGGTCTGCCAAAACGAATGATCGAAGAGGCCTCAGCGCGCGAACAGTCGCTGATCGACCAGGGCAAGCGTGTCATCGTTGGTGTCAACAAGTACAAACTGGATCACGAAGACGAAACCGATGTACTTGAGATCGACAACGTGATGGTGCGTAACGAGCAAATTGCTTCGCTGGAACGCATTCGCGCCACCCGTGATGATGCCGCCGTAACCGCCGCGTTGAACGCCCTGACTCACGCCGCACAGCATAAC  >  NC_000913/3062119‑3062382

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