Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	712,150	A→T	100%	L218Q (CTG→CAG)	ybfF ←	acyl‑CoA esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	712,150	0	A	T	97.9%	320.8 / ‑1.8	96	L218Q (CTG→CAG)	ybfF	acyl‑CoA esterase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (38/56);  total (39/57)
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.41e-01

GACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCT  >  NC_000913/712055‑712245                                                                                                |                                                                                                gACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTaa                                                                                             >  1:7008/1‑100 (MQ=255)    ggATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatc                                                                                          >  2:2425876/1‑100 (MQ=255)    ggATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatc                                                                                          <  1:581444/100‑1 (MQ=255)    ggATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatc                                                                                          <  1:576940/100‑1 (MQ=255)    ggATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatc                                                                                          >  2:559420/1‑100 (MQ=255)      aTAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAatcat                                                                                        >  1:282325/1‑100 (MQ=255)        aGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAAtcatca                                                                                      <  1:244179/100‑1 (MQ=255)         gcgcGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCGGTGGAAATTGAGCCTGTAAATCATCAc                                                                                     >  2:1032882/1‑100 (MQ=255)         gcgcGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCCGTGGAAATTGAGCCTGTAAatcat                                                                                        >  2:1273327/1‑97 (MQ=255)          cgcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACg                                                                                    <  1:240539/100‑1 (MQ=255)           gcgTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACgg                                                                                   <  1:942308/100‑1 (MQ=255)             gTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTa                                                                                 >  2:1608840/1‑100 (MQ=255)              tAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGATCCTGTAAATCATCACGGTAc                                                                                <  2:1560844/100‑1 (MQ=255)              tAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTAc                                                                                <  1:233496/100‑1 (MQ=255)                 tACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGc                                                                             >  2:1868317/1‑100 (MQ=255)                 tACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGc                                                                             <  1:690702/100‑1 (MQ=255)                  aCCGCATCCGGTTTTTCAGCATGTACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCt                                                                            <  1:1559813/100‑1 (MQ=255)                  aCCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCt                                                                            >  1:356124/1‑100 (MQ=255)                  aCCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCt                                                                            >  1:434943/1‑100 (MQ=255)                   ccGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTc                                                                           <  2:374164/100‑1 (MQ=255)                    cGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCg                                                                          <  2:998269/100‑1 (MQ=255)                     gCATCCGGTTTTTCGGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGc                                                                         <  2:312818/100‑1 (MQ=255)                     gCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGc                                                                         <  2:1193843/100‑1 (MQ=255)                     gCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGc                                                                         >  2:2014291/1‑100 (MQ=255)                       aTCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTa                                                                       >  1:1000334/1‑100 (MQ=255)                        tCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTaa                                                                      <  1:1888267/100‑1 (MQ=255)                          cGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACa                                                                    <  2:115384/100‑1 (MQ=255)                           ggTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACAt                                                                   >  1:237369/1‑100 (MQ=255)                            gTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACCTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATa                                                                  <  1:599148/100‑1 (MQ=255)                             tttttCAGCATGGACCCAGTTTCCCCCGCTTGCAATCACATGCGCTCTTTCCTGTGGAAATTGAGCCTGTAAATCATCACGCTACTACTCGCTACCATAc                                                                 <  2:410336/100‑1 (MQ=255)                                ttCATCATGGACCCAGTGTCTCGCTCCTGCAATTACTTGCTCCCGTTTCTGTGGAAATTGAGCCTGTAAATCATCACGTTACTGCTCGCTAACATACgg                                                               <  2:2033502/99‑1 (MQ=255)                                ttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGga                                                              <  1:2042993/100‑1 (MQ=255)                                ttCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGga                                                              >  2:1226169/1‑100 (MQ=255)                                 tCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGgag                                                             <  1:1873149/100‑1 (MQ=255)                                      aTGGACGCAGTGAGCCGCGGCTGTAATTACATGCGCCTGTGCCTGTGGAAATTTTGCCTGTAAATCATCACGTTACTGCTCGCTAACATACGGAGAATTg                                                        <  2:386466/100‑1 (MQ=255)                                         gACCCAGTGACCCGCGCCTGCAATCACATGCGCCCTTGCCTGTGAAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgccg                                                     <  1:748888/100‑1 (MQ=255)                                         gACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTATAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTgccg                                                     >  2:1635301/1‑100 (MQ=255)                                             cAGTGCCCCGCCCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCCCGGTACTGCTCGCTAACATACGGACCATTgccgcc                                                   >  2:770527/1‑98 (MQ=255)                                              aGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAgg                                                <  1:131835/100‑1 (MQ=255)                                              aGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAgg                                                <  1:1903350/100‑1 (MQ=255)                                              aGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAgg                                                <  2:216128/100‑1 (MQ=255)                                               gTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAggg                                               <  1:828497/100‑1 (MQ=255)                                               gTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAgg                                                <  2:216250/99‑1 (MQ=255)                                                 gACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGAt                                             >  1:234834/1‑100 (MQ=255)                                                 gACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGAt                                             <  2:1830138/100‑1 (MQ=255)                                                  aCCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATa                                            <  1:931315/100‑1 (MQ=255)                                                   cccGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATaa                                           >  1:2122263/1‑100 (MQ=255)                                                   cccGAGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAggg                                               >  2:1409671/1‑96 (MQ=255)                                                      gcgcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAca                                        <  2:1380137/100‑1 (MQ=255)                                                      gcgcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAca                                        <  2:2111134/100‑1 (MQ=255)                                                        gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                      <  2:934338/100‑1 (MQ=255)                                                        gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                      <  1:522393/100‑1 (MQ=255)                                                        gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                      >  1:1089180/1‑100 (MQ=255)                                                        gcCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagg                                      <  1:1340081/100‑1 (MQ=255)                                                            gCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcag                                  <  1:714153/100‑1 (MQ=255)                                                              aaTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAAcagggcaggg                                <  2:1692071/100‑1 (MQ=255)                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                               >  1:163409/1‑100 (MQ=255)                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                               >  1:1385609/1‑100 (MQ=255)                                                               aTCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGt                               >  1:2220574/1‑100 (MQ=255)                                                                tCACATGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTg                              <  2:1898704/100‑1 (MQ=255)                                                                  acaTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGAt                            >  1:1636219/1‑100 (MQ=255)                                                                   caTGCGCCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGAt                            >  2:1084784/1‑99 (MQ=255)                                                                      gcgcCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCa                        <  1:1273335/100‑1 (MQ=255)                                                                        gcCCGTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATg                      <  1:2410893/100‑1 (MQ=255)                                                                            gTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCgg                  <  1:1852093/100‑1 (MQ=255)                                                                            gTGCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCgg                  <  1:1942776/100‑1 (MQ=255)                                                                            gTGCCTGTGGAAATTGAGCCTGTAAATCAGCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCgg                  <  2:7008/100‑1 (MQ=255)                                                                             tGCCTGTGGAAATTGTGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCggg                 >  1:684865/1‑100 (MQ=255)                                                                              gCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGa                <  1:897477/100‑1 (MQ=255)                                                                              gCCTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGa                <  2:1266275/100‑1 (MQ=255)                                                                                cTGTGGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAtt              <  1:1186163/100‑1 (MQ=255)                                                                                 tgtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAttt             >  1:41041/1‑100 (MQ=255)                                                                                 tgtgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGAttt             <  1:1504944/100‑1 (MQ=255)                                                                                   tgGAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtc           >  1:480221/1‑100 (MQ=255)                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct          >  2:194488/1‑100 (MQ=255)                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct          >  1:96522/1‑100 (MQ=255)                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct          >  1:876601/1‑100 (MQ=255)                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtct          <  1:2183510/100‑1 (MQ=255)                                                                                    ggAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCAGCCAGGGATAAACAGGGCAGGGTGATGCCATGCCGGGATTTtct          >  2:437902/1‑100 (MQ=255)                                                                                     gAAATTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTtctc         <  2:1257563/100‑1 (MQ=255)                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGTTAACATACGGAGAATTGCCGCCAGGCATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc        <  2:275435/100‑1 (MQ=255)                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc        <  2:587815/100‑1 (MQ=255)                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc        >  2:871112/1‑100 (MQ=255)                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc        <  1:1863605/100‑1 (MQ=255)                                                                                      aaaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTcc        <  1:879743/100‑1 (MQ=255)                                                                                       aaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTccc       <  1:786556/100‑1 (MQ=255)                                                                                       aaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTccc       >  1:2379745/1‑100 (MQ=255)                                                                                       aaTTGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTccc       >  1:107336/1‑100 (MQ=255)                                                                                         ttGAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCaa     <  1:392666/100‑1 (MQ=255)                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc   <  2:1310756/100‑1 (MQ=255)                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc   >  2:1113267/1‑100 (MQ=255)                                                                                           gAGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcc   <  2:595314/100‑1 (MQ=255)                                                                                            aGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATGAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAAcgt  <  1:1730863/100‑3 (MQ=255)                                                                                            aGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCt  >  1:474737/1‑100 (MQ=255)                                                                                            aGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCt  >  1:1613230/1‑100 (MQ=255)                                                                                            aGCCTGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCt  >  1:1595999/1‑100 (MQ=255)                                                                                                |                                                                                                GACGGATAGCGCGTAATACCGCATCCGGTTTTTCAGCATGGACCCAGTGACCCGCGCCTGCAATCACATGCGCCCGTGCCTGTGGAAATTGAGCCAGTAAATCATCACGGTACTGCTCGCTAACATACGGAGAATTGCCGCCAGGGATAAACAGGGCAGGGTGATCCCATGCCGGGATTTTCTCCCAACCT  >  NC_000913/712055‑712245

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