Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 364,024 C→T R761H (CGC→CAC)  lacZ ← beta‑D‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913364,0240CT94.4% 280.5 / 3.4 89R761H (CGC→CAC) lacZbeta‑D‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (38/46);  total (40/49)
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.46e-01

GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCG  >  NC_000913/363927‑364117
                                                                                                 |                                                                                             
gTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACtggtgg                                                                                             <  2:164336/100‑1 (MQ=255)
gTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACtggtgg                                                                                             <  1:322228/100‑1 (MQ=255)
  cGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGtt                                                                                           <  2:162388/100‑1 (MQ=255)
    ttATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTaa                                                                                         <  2:337063/100‑1 (MQ=255)
    ttATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTaa                                                                                         <  1:702734/100‑1 (MQ=255)
     tATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTaaa                                                                                        >  2:521226/1‑100 (MQ=255)
     tATCCAGCGGTGACCGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTCGTTTTCTAACGCAAATCCACATCTGTGAAAGAAAGCCTGTCTAGTGGTTaaa                                                                                        >  2:1937112/1‑100 (MQ=255)
       tCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAAtt                                                                                      <  2:571196/100‑1 (MQ=255)
        ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGCAAGAAAGCCTGACTGGTGGTTAGATTg                                                                                     >  1:1573098/1‑100 (MQ=255)
        ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTg                                                                                     >  1:1294315/1‑100 (MQ=255)
        ccAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTg                                                                                     <  2:109262/100‑1 (MQ=255)
         cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTg                                                                                     <  1:934729/99‑1 (MQ=255)
         cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGc                                                                                    <  1:390606/100‑1 (MQ=255)
         cAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGc                                                                                    >  1:1841624/1‑100 (MQ=255)
           gCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCa                                                                                  <  2:1633841/100‑1 (MQ=255)
           gCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCa                                                                                  <  1:1217408/100‑1 (MQ=255)
              gTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACg                                                                               >  1:1374045/1‑100 (MQ=255)
              gTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACg                                                                               <  2:456104/100‑1 (MQ=255)
                gCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCt                                                                             >  1:1796900/1‑100 (MQ=255)
                 cACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCGCATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGc                                                                              >  1:317247/1‑98 (MQ=255)
                 cACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCtt                                                                            <  2:1204903/100‑1 (MQ=255)
                   cGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttat                                                                          >  2:1413716/1‑100 (MQ=255)
                    gggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatt                                                                         >  1:1931160/1‑100 (MQ=255)
                    gggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatt                                                                         >  2:1012739/1‑100 (MQ=255)
                     ggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatta                                                                        <  2:1772541/100‑1 (MQ=255)
                     ggTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatta                                                                        >  2:428138/1‑100 (MQ=255)
                      gTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCttatta                                                                        >  2:80425/1‑99 (MQ=255)
                      gTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAc                                                                       <  2:1742452/100‑1 (MQ=255)
                       tGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAcc                                                                      <  1:842276/100‑1 (MQ=255)
                       tGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAcc                                                                      >  2:1897972/1‑100 (MQ=255)
                       tGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAcc                                                                      >  2:826158/1‑100 (MQ=255)
                        gAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTAccc                                                                     <  1:1034245/100‑1 (MQ=255)
                         aaCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCa                                                                    >  2:211763/1‑100 (MQ=255)
                          aCTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAg                                                                   <  1:1543046/100‑1 (MQ=255)
                            tGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCt                                                                 >  1:1600544/1‑100 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTc                                                                <  1:25092/100‑1 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTc                                                                >  1:177061/1‑100 (MQ=255)
                             gATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTc                                                                <  2:991561/100‑1 (MQ=255)
                               tCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGa                                                              <  2:1178029/100‑1 (MQ=255)
                                cgcgCAGCGGCGTCAGCAGTTGTTTTTTATCGCTAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGAt                                                             <  2:46666/100‑1 (MQ=255)
                                 gcgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGAtn                                                            <  1:1880723/100‑2 (MQ=255)
                                 gcgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATg                                                            >  1:318441/1‑100 (MQ=255)
                                  cgcAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGc                                                           <  2:957563/100‑1 (MQ=255)
                                    cAGCGGCGTCAGCAGTTGTTTTTTATCGCGAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTTAATTGCCAACGCTTATTACCCAGCTCGATGCaa                                                         <  1:1173279/100‑1 (MQ=255)
                                      gcggcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCaaaa                                                       >  1:147800/1‑100 (MQ=255)
                                         gcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATc                                                    <  2:776211/100‑1 (MQ=255)
                                         gcgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATc                                                    <  2:442911/100‑1 (MQ=255)
                                          cgTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATcc                                                   <  2:965179/100‑1 (MQ=255)
                                            tCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAt                                                 <  2:1563401/100‑1 (MQ=255)
                                             cagcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAtt                                                <  2:799811/100‑1 (MQ=255)
                                             cagcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCAtt                                                >  1:1484222/1‑100 (MQ=255)
                                               gcagTTTTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTc                                              <  1:332930/100‑1 (MQ=255)
                                               gcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGTTTATTACCCAGCTCGATGCAAAAATCCATTTc                                              <  2:1841636/100‑1 (MQ=255)
                                               gcagTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTc                                              >  2:1082875/1‑100 (MQ=255)
                                                  gttgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCt                                           <  2:1307251/100‑1 (MQ=255)
                                                    tgttTTTTATTGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTACTACCCAGCTCGATGCAAAAATCCATTTCGCtgg                                         <  1:1274323/100‑1 (MQ=255)
                                                    tgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtgg                                         >  2:1274328/1‑100 (MQ=255)
                                                    tgttTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtgg                                         <  2:689586/100‑1 (MQ=255)
                                                      ttttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtg                                       >  2:283605/1‑100 (MQ=255)
                                                       tttttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCtggtgg                                      <  2:177062/100‑1 (MQ=255)
                                                         tttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTc                                    <  1:810227/100‑1 (MQ=255)
                                                         tttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTc                                    <  1:1665234/100‑1 (MQ=255)
                                                          ttATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCa                                   <  1:957242/100‑1 (MQ=255)
                                                              cGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATg                               <  2:1324623/100‑1 (MQ=255)
                                                               gttatgccACATCTCTGAAAGAAATCCTTACTTTCGGTTAAATGGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGc                              <  2:1330047/94‑1 (MQ=255)
                                                               gacAATCCACATCCGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGc                              >  2:1308836/3‑100 (MQ=255)
                                                               gCCAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGc                              >  2:1009318/1‑100 (MQ=255)
                                                                ccAATCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCg                             <  1:1897959/100‑1 (MQ=255)
                                                                  aaTCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTCCCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCggg                           <  2:362504/100‑1 (MQ=255)
                                                                   aTCCACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGa                          <  1:1128807/100‑1 (MQ=255)
                                                                     ccACATCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATg                        >  1:675359/1‑100 (MQ=255)
                                                                          tCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTg                   <  2:456068/100‑1 (MQ=255)
                                                                          tCTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGATTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCg                     >  1:799007/1‑98 (MQ=255)
                                                                           cTGTGAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTgg                  >  2:532483/1‑100 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                 <  1:283603/100‑1 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                 >  2:117472/1‑100 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                 >  2:1528714/1‑100 (MQ=255)
                                                                            tgtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTggg                 >  1:190266/1‑100 (MQ=255)
                                                                             gtgAAAGAAAGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGa                >  2:1250173/1‑100 (MQ=255)
                                                                               gaaagaaaGCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcg              >  1:188680/1‑100 (MQ=255)
                                                                                aaagaaagCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGAcgc             >  2:700619/1‑100 (MQ=255)
                                                                                  agaaagGCGGACGGGTGGGTAAAGTGCCAACGGAAATTACCAAGCTGGATGCAAAAATCCATATCGCTGGTGGTCAGATGCGGGAAGGCGTGGTACgcgg           >  1:1911122/1‑100 (MQ=255)
                                                                                   gaaagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAGAAATCCATTGCGCTGATGGTCAGAGGCGGGATGGCGTGGGACgcggc          >  1:1312789/1‑100 (MQ=255)
                                                                                     aagCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGcggcgg        <  2:1374053/100‑1 (MQ=255)
                                                                                      agCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCggg       <  2:1796911/100‑1 (MQ=255)
                                                                                       gCCTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCgggg      <  1:1713735/100‑1 (MQ=255)
                                                                                         cTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg    <  2:205155/100‑1 (MQ=255)
                                                                                         cTGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAg    >  1:118389/1‑100 (MQ=255)
                                                                                          tGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGc   >  2:1106326/1‑100 (MQ=255)
                                                                                          tGACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGc   >  1:798202/1‑100 (MQ=255)
                                                                                           gACTGGTGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATGTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCg  <  1:1082874/100‑1 (MQ=255)
                                                                                                 |                                                                                             
GTCGTTATCCAGCGGTGCACGGGTGAACTGATCGCGCAGCGGCGTCAGCAGTTGTTTTTTATCGCCAATCCACATCTGTGAAAGAAAGCCTGACTGGCGGTTAAATTGCCAACGCTTATTACCCAGCTCGATGCAAAAATCCATTTCGCTGGTGGTCAGATGCGGGATGGCGTGGGACGCGGCGGGGAGCG  >  NC_000913/363927‑364117

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