Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,953,902 T→C Q377R (CAG→CGG)  b1872 ← biotin sulfoxide reductase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,953,9020TC99.5% 610.7 / ‑4.7 195Q377R (CAG→CGG) b1872biotin sulfoxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (110/84);  total (111/84)
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.00e-01

Only 100 of 205 total aligned reads displayed.

TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953813‑1953991
                                                                                         |                                                                                         
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                        >  1:657933/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                        >  2:100134/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                        >  1:3118854/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                        >  1:1787584/1‑93 (MQ=255)
 ttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGccc                                                                                       <  1:1603950/93‑1 (MQ=255)
   ccGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTa                                                                                     <  1:3372108/93‑1 (MQ=255)
     ggcggcACCCCGCAAACCGGCGTCGGGTTGCCACGTTTGAAGTATTGATAAGAGAAACC‑AAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  <  2:6771800/93‑1 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAAc                                                                                   >  1:6513598/1‑93 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTCAc                                                                                   >  1:2864466/1‑93 (MQ=255)
      gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  <  1:5056429/93‑1 (MQ=255)
        ctcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGATAAACCAAATCCGCCGCCCGGTGTACCAATCCTCCCTAACATc                                                                                <  2:4853586/91‑1 (MQ=255)
          caCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgc                                                                              <  2:2572795/93‑1 (MQ=255)
           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgcc                                                                             >  2:4151317/1‑93 (MQ=255)
              ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  1:2553420/93‑1 (MQ=255)
              ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  1:1908760/93‑1 (MQ=255)
              ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  1:6043590/93‑1 (MQ=255)
                 ccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAACCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGc                                                                       <  1:3691191/93‑1 (MQ=255)
                 ccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGc                                                                       <  1:4908906/93‑1 (MQ=255)
                 ccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGc                                                                       <  2:1194421/93‑1 (MQ=255)
                   aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGt                                                                     >  1:352598/1‑93 (MQ=255)
                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGTCCGGTGTTCCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                  >  2:614877/1‑93 (MQ=255)
                       cgcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAcc                                                                 <  2:817558/93‑1 (MQ=255)
                          gTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGc                                                              >  2:203441/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                             >  1:5162702/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAAACAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                             >  2:7799811/1‑93 (MQ=255)
                           tCGGGTTGCCAACGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                             >  1:1396837/1‑93 (MQ=255)
                            cGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCAt                                                            >  1:4280063/1‑93 (MQ=255)
                             gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           >  2:2685648/1‑93 (MQ=255)
                             gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           >  2:4574173/1‑93 (MQ=255)
                             gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           >  1:487475/1‑93 (MQ=255)
                              ggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                          <  2:7249088/93‑1 (MQ=255)
                               gTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCa                                                         >  1:3572874/1‑93 (MQ=255)
                                 tGCCACCGTTGGAGTAATGATAAGAGAAACCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       <  2:4855204/93‑1 (MQ=255)
                                 tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       <  2:4354710/93‑1 (MQ=255)
                                    cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                    >  2:1258562/1‑93 (MQ=255)
                                      ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                  <  2:1698895/93‑1 (MQ=255)
                                      ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                  >  2:2237068/1‑93 (MQ=255)
                                      ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                  <  1:1331742/93‑1 (MQ=255)
                                       cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                 >  2:6553441/1‑93 (MQ=255)
                                       cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                 <  1:4474155/93‑1 (MQ=255)
                                        gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                >  1:782093/1‑93 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTGCCAATCCGCCCTAACATCCCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                               <  2:6688342/93‑1 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                               >  1:212102/1‑93 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                               >  1:7461461/1‑93 (MQ=255)
                                           ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTc                                             >  2:5898865/1‑93 (MQ=255)
                                               tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  2:5362987/1‑93 (MQ=255)
                                               tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  2:1967535/1‑93 (MQ=255)
                                               tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  2:4351244/1‑93 (MQ=255)
                                                aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                        <  1:2313016/93‑1 (MQ=255)
                                                 aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                       <  1:3864867/93‑1 (MQ=255)
                                                 aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGAGTTTCTGTTCACCata                                       >  1:7037893/1‑93 (MQ=255)
                                                  tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      >  2:3281180/1‑93 (MQ=255)
                                                  tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      >  2:1642570/1‑93 (MQ=255)
                                                     tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                   >  1:232765/1‑93 (MQ=255)
                                                     tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                   <  1:6840894/93‑1 (MQ=255)
                                                     tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                   >  1:6023759/1‑93 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACAAGCATCCAGTGTTTCTGTTCGCCATATTGCt                                  <  2:5487293/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  1:5332258/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  2:2997779/93‑1 (MQ=255)
                                                        gataaaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                <  1:5858851/90‑1 (MQ=255)
                                                         agaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgc                               <  2:407505/93‑1 (MQ=255)
                                                          gaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGCTTCTGTTCACCATATTGCTGgcg                              >  1:2034694/1‑93 (MQ=255)
                                                               cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                         <  1:4407215/93‑1 (MQ=255)
                                                                aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                        <  1:937945/93‑1 (MQ=255)
                                                                aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                        <  1:276489/93‑1 (MQ=255)
                                                                 aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                       <  1:364360/93‑1 (MQ=255)
                                                                  aaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                      <  1:2616502/93‑1 (MQ=255)
                                                                  aaCCGCCCCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCCTCCCCAGCCTCCGGTGTTTCTTTTCACCATATTGCTGGCGCTGAATTc                      <  2:5560081/93‑1 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     >  1:1070440/1‑93 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     >  2:6641142/1‑93 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     >  1:7588699/1‑93 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     <  1:4811462/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGTCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                    <  2:5660070/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                    >  1:5469356/1‑93 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                    <  2:4726110/93‑1 (MQ=255)
                                                                     cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  2:4117836/93‑1 (MQ=255)
                                                                     cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  2:5280396/93‑1 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  1:654145/93‑1 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  <  2:1369081/93‑1 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  >  2:6304325/1‑93 (MQ=255)
                                                                       ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                 <  1:7900732/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc               <  1:5675128/93‑1 (MQ=255)
                                                                          cccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt              <  2:7358794/93‑1 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             >  1:4411392/1‑93 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             <  2:7533896/93‑1 (MQ=255)
                                                                            cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            <  2:2534989/93‑1 (MQ=255)
                                                                            cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  1:1166817/1‑93 (MQ=255)
                                                                             ggTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGcc           >  2:6702375/1‑93 (MQ=255)
                                                                              gtgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa          <  1:7766969/93‑1 (MQ=255)
                                                                              gtgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa          >  2:3336137/1‑93 (MQ=255)
                                                                               tgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt         >  2:4332088/1‑93 (MQ=255)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAGTTCCCCAGCCTGCCATCa       >  1:4622713/1‑93 (MQ=255)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  2:5748705/93‑1 (MQ=255)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  2:5828577/93‑1 (MQ=255)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  1:1003645/93‑1 (MQ=255)
                                                                                  aCCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAg      <  1:2422516/93‑1 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc     >  1:3684906/1‑93 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc     >  2:3640808/1‑93 (MQ=255)
                                                                                     aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   <  2:1792093/93‑1 (MQ=255)
                                                                                     aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   >  1:1699394/1‑93 (MQ=255)
                                                                                         |                                                                                         
TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953813‑1953991

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