Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,660,346 A→C M760L (ATG→CTG)  ynfE → putative selenate reductase, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,660,3460AC94.1% 436.2 / 6.2 138M760L (ATG→CTG) ynfEputative selenate reductase, periplasmic
Reads supporting (aligned to +/- strand):  ref base A (5/3);  new base C (60/68);  total (65/73)
Fisher's exact test for biased strand distribution p-value = 4.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 138 total aligned reads displayed.

GATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCA  >  NC_000913/1660249‑1660441
                                                                                                 |                                                                                               
gATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTg                                                                                               >  2:2609458/1‑100 (MQ=255)
  ttGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTggg                                                                                             >  2:909622/1‑100 (MQ=255)
  ttGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTggg                                                                                             <  2:673953/100‑1 (MQ=255)
  ctGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTggg                                                                                             <  1:451416/99‑1 (MQ=255)
    gCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGcc                                                                                           >  2:507074/1‑100 (MQ=255)
    gCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGcc                                                                                           <  1:2458509/100‑1 (MQ=255)
     ccAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCa                                                                                          >  2:3032807/1‑100 (MQ=255)
     ccAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCa                                                                                          >  2:121341/1‑100 (MQ=255)
      cAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAg                                                                                         <  1:2669743/100‑1 (MQ=255)
       aaCGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAgg                                                                                        <  1:83095/100‑1 (MQ=255)
           gCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGc                                                                                    <  1:636690/100‑1 (MQ=255)
           gCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGc                                                                                    >  1:66318/1‑100 (MQ=255)
            cGATATGGTGCGGGTGTTTGAACATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGcc                                                                                   >  1:3485776/1‑100 (MQ=255)
            cGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGcc                                                                                   >  2:2669616/1‑100 (MQ=255)
             gATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCt                                                                                  >  1:320020/1‑100 (MQ=255)
               tatGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTgg                                                                                <  1:2850944/100‑1 (MQ=255)
               cttGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTgg                                                                                <  1:1136167/98‑1 (MQ=255)
                atGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGc                                                                               <  1:2710838/100‑1 (MQ=255)
                 tGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTTCCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCa                                                                              >  2:3171010/1‑100 (MQ=255)
                 tGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCa                                                                              <  1:70537/100‑1 (MQ=255)
                 tGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCa                                                                              <  1:2604386/100‑1 (MQ=255)
                 tGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCa                                                                              <  1:2604108/100‑1 (MQ=255)
                  ggTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCAc                                                                             >  1:1545301/1‑100 (MQ=255)
                  ggTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCAc                                                                             <  2:1435065/100‑1 (MQ=255)
                  ggTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCAc                                                                             >  2:1970669/1‑100 (MQ=255)
                   gTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACg                                                                            >  1:2672139/1‑100 (MQ=255)
                     gCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGTCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAg                                                                          <  1:1417976/100‑1 (MQ=255)
                     gCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAg                                                                          <  1:2384505/100‑1 (MQ=255)
                     gCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAg                                                                          >  1:3983683/1‑100 (MQ=255)
                     gCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAg                                                                          <  1:1544889/100‑1 (MQ=255)
                     gCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAg                                                                          <  1:1080276/100‑1 (MQ=255)
                      cGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAgg                                                                         <  2:2982585/100‑1 (MQ=255)
                       gggTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGc                                                                        <  2:3416554/100‑1 (MQ=255)
                        ggTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGcc                                                                       <  1:3271867/100‑1 (MQ=255)
                        ggTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGcc                                                                       >  2:2625655/1‑100 (MQ=255)
                        ggTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCAAAAGTAACACCACGTATACTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGcc                                                                       >  2:3745874/1‑100 (MQ=255)
                         gtgtTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCa                                                                      <  1:1756995/100‑1 (MQ=255)
                          tgtTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCaa                                                                     >  1:1855157/1‑100 (MQ=255)
                           gtTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAat                                                                    <  1:2975658/100‑1 (MQ=255)
                            tttAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAata                                                                   <  1:3396967/100‑1 (MQ=255)
                             ttAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAatat                                                                  >  2:1377425/1‑100 (MQ=255)
                              tAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATg                                                                 <  2:1125743/100‑1 (MQ=255)
                              tAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATg                                                                 >  1:1197060/1‑100 (MQ=255)
                               aaCCAACGTGGCGAAGTACGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCCCGAGGCCAacatat                                                                >  2:3103250/1‑98 (MQ=255)
                                aCCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTc                                                               <  1:1336165/100‑1 (MQ=255)
                                 ccATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTTTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCt                                                              <  2:3896240/100‑1 (MQ=255)
                                 ccATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCt                                                              >  1:2076667/1‑100 (MQ=255)
                                  cATCGTGGCGAAGTTCGGGTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTg                                                             <  1:3826510/100‑1 (MQ=255)
                                  cATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTg                                                             <  2:531688/100‑1 (MQ=255)
                                   aTCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTgg                                                            >  2:2230339/1‑100 (MQ=255)
                                    tCGTGGCGAAGTTCGGCTACCAGTCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGc                                                           <  1:2322933/100‑1 (MQ=255)
                                    tCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGc                                                           <  2:1237323/100‑1 (MQ=255)
                                      gTGGCGAAGTTCGGCGACCAGCGAAATTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGAGCCAGGGAGCCTGGCGCGTGGCCACTATGTCTGGCGa                                                         >  2:826349/1‑100 (MQ=255)
                                        ggCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATa                                                       <  2:1555688/100‑1 (MQ=255)
                                        ggCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATa                                                       >  1:1089988/1‑100 (MQ=255)
                                        ggCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATa                                                       >  2:295786/1‑100 (MQ=255)
                                          cGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATaaa                                                     >  2:2316920/1‑100 (MQ=255)
                                            aaGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAAt                                                   <  1:3360511/100‑1 (MQ=255)
                                             aGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATc                                                  >  1:233202/1‑100 (MQ=255)
                                             aGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATc                                                  <  2:1052979/100‑1 (MQ=255)
                                             aGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATc                                                  >  1:477435/1‑100 (MQ=255)
                                              gTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCg                                                 >  2:1885858/1‑100 (MQ=255)
                                              gTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCg                                                 >  1:149526/1‑100 (MQ=255)
                                                tCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGAc                                               >  2:415899/1‑100 (MQ=255)
                                                tCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGAc                                               <  2:3810222/100‑1 (MQ=255)
                                                tCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGAc                                               <  2:2179839/100‑1 (MQ=255)
                                                 cGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGAcc                                              <  2:3969542/100‑1 (MQ=255)
                                                 cGGCCACCAGCCACAGTAACACCACGTACTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGAcc                                              <  1:1254506/100‑1 (MQ=255)
                                                  ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCa                                             >  1:3508113/1‑100 (MQ=255)
                                                  ggCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCa                                             <  2:696888/100‑1 (MQ=255)
                                                   gCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCAt                                            <  1:1762268/100‑1 (MQ=255)
                                                     tACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATgg                                          <  1:1188885/100‑1 (MQ=255)
                                                       ccAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATggcg                                        >  1:2274586/1‑100 (MQ=255)
                                                        ccgCCAAAGTAACACCACGTCTTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATggcgg                                       <  1:2232290/98‑1 (MQ=255)
                                                        cAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATggcgg                                       <  2:2421778/100‑1 (MQ=255)
                                                        cAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATggcgg                                       >  2:447734/1‑100 (MQ=255)
                                                          gCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCt                                     >  2:723056/1‑100 (MQ=255)
                                                            cAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCtgt                                   <  2:902737/100‑1 (MQ=255)
                                                            cAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCCGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCtgt                                   >  1:2285522/1‑100 (MQ=255)
                                                               aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCtgtgtg                                >  2:111169/1‑100 (MQ=255)
                                                               aGTAACACCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCtgtgtg                                <  1:3707044/100‑1 (MQ=255)
                                                                   acacCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATa                            <  2:816090/100‑1 (MQ=255)
                                                                   acacCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATa                            <  1:208351/100‑1 (MQ=255)
                                                                   acacCACGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATa                            <  1:2160774/100‑1 (MQ=255)
                                                                     accacGTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACg                          >  1:3160513/1‑100 (MQ=255)
                                                                     acaacGTATTCTCACTGAAGTTAGCGCTATGGGCCAGGGAGCGTAGCACGAAGGCAATAAGTCTGGCGATAAAATCGACCATGGCGGATATGTGCATACg                          >  2:3017028/4‑100 (MQ=255)
                                                                          gTATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAAc                     <  1:1885804/100‑1 (MQ=255)
                                                                           tATTCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAAcc                    <  1:2173268/100‑1 (MQ=255)
                                                                             ttCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCAc                  >  2:380615/1‑100 (MQ=255)
                                                                             ttCTCCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCAc                  <  1:2218574/100‑1 (MQ=255)
                                                                               ctcCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCActc                >  1:846016/1‑100 (MQ=255)
                                                                               ctcCCTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCActc                >  1:93979/1‑100 (MQ=255)
                                                                                 cccTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTg              >  1:2101004/1‑100 (MQ=255)
                                                                                 cccTGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTg              >  1:2013/1‑100 (MQ=255)
                                                                                    tGGAGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTgcgc           <  1:2390943/100‑1 (MQ=255)
                                                                                       aGTTAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCt        >  2:3951237/1‑100 (MQ=255)
                                                                                         ttAGCGCTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATCAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTc      <  1:3878502/100‑1 (MQ=255)
                                                                                           aGCGCTCTGGGCCAGGGGGCCGGGCACGAGGCCAATATTTCTGGCGATAAAATCGACCATGGCGGCTGTGGGAATACGCTAACCACTCTGCGCCCTTcac    <  2:3118120/100‑1 (MQ=255)
                                                                                            gcgcTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcacc   <  2:2974162/100‑1 (MQ=255)
                                                                                             cgcTCTGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTcacca  <  2:2175889/100‑1 (MQ=255)
                                                                                                 |                                                                                               
GATTGCCAACGGCGATATGGTGCGGGTGTTTAACCATCGTGGCGAAGTTCGGCTACCAGCCAAAGTAACACCACGTATTCTCCCTGGAGTTAGCGCTATGGGCCAGGGAGCCTGGCACGAGGCCAATATGTCTGGCGATAAAATCGACCATGGCGGCTGTGTGAATACGCTAACCACTCTGCGCCCTTCACCA  >  NC_000913/1660249‑1660441

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