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,9020TC96.5% 347.0 / 1.5 118Q377R (CAG→CGG) b1872biotin sulfoxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (3/1);  new base C (54/58);  total (57/61)
Fisher's exact test for biased strand distribution p-value = 3.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01

Only 100 of 125 total aligned reads displayed.

TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAT  >  NC_000913/1953813‑1953990
                                                                                         |                                                                                        
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  2:2245044/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  1:2138454/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       <  2:2245043/93‑1 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  2:2685911/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  2:3164603/1‑93 (MQ=255)
tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                       >  1:3226506/1‑93 (MQ=255)
    cGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTaa                                                                                   >  2:1688574/1‑93 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                  >  1:4533467/1‑93 (MQ=255)
     ggCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAc                                                                                  >  2:1572830/1‑93 (MQ=255)
       cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                <  2:4085077/93‑1 (MQ=255)
        gacaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAATCCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATc                                                                               <  2:3122183/91‑1 (MQ=255)
        agcaCGCCGCCAACACGCGTCGTTTTGCCACCGTTGGAGTAATGATAAGAGAACCCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               <  2:4717585/93‑1 (MQ=255)
        agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               >  2:3522581/1‑93 (MQ=255)
        agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               <  2:355880/93‑1 (MQ=255)
        agcaCGCCGCCAACACGAGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                               >  1:4296456/1‑93 (MQ=255)
         gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                              >  1:738633/1‑93 (MQ=255)
          caCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgc                                                                             <  1:1331805/93‑1 (MQ=255)
           aCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgcc                                                                            <  1:1783775/93‑1 (MQ=255)
            cgccgcCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgccg                                                                           >  1:4568294/1‑93 (MQ=255)
             gccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcgccgc                                                                          <  1:2841678/93‑1 (MQ=255)
               cgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCa                                                                        <  1:3574010/93‑1 (MQ=255)
                   aaCACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGt                                                                    <  2:69178/93‑1 (MQ=255)
                    acacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                   >  1:904858/1‑93 (MQ=255)
                     cacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCa                                                                  <  2:3803073/93‑1 (MQ=255)
                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                 >  1:1675169/1‑93 (MQ=255)
                      acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                 <  1:340722/93‑1 (MQ=255)
                       cgcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAcc                                                                <  1:1008799/93‑1 (MQ=255)
                        gcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCa                                                               >  2:2440505/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                            >  1:1814200/1‑93 (MQ=255)
                           tCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCa                                                            <  2:2626385/93‑1 (MQ=255)
                              ggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATcc                                                         <  1:1553002/93‑1 (MQ=255)
                                 tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                      <  1:4225181/93‑1 (MQ=255)
                                  gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                     >  2:35996/1‑93 (MQ=255)
                                  gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                     >  2:2900200/1‑93 (MQ=255)
                                   ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                    <  2:1070704/93‑1 (MQ=255)
                                   ccACCGGTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                    <  2:4534237/93‑1 (MQ=255)
                                    cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                   >  2:3708609/1‑93 (MQ=255)
                                      ccGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTc                                                 <  2:4599681/93‑1 (MQ=255)
                                       cGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCt                                                <  2:3904043/93‑1 (MQ=255)
                                        gTTGGAGTAATGATAAGAGAAACCGAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                               >  2:4554996/1‑93 (MQ=255)
                                        gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                               <  1:3716307/93‑1 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                              <  1:460033/93‑1 (MQ=255)
                                         ttGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGt                                              >  2:2885484/1‑93 (MQ=255)
                                          tGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                             <  2:847575/93‑1 (MQ=255)
                                           ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTc                                            <  1:2657613/93‑1 (MQ=255)
                                           ggAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTc                                            >  2:3177407/1‑93 (MQ=255)
                                            gAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGCGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCa                                           >  1:1891094/1‑93 (MQ=255)
                                              gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                         >  1:2781113/1‑93 (MQ=255)
                                              gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                         <  1:2233343/93‑1 (MQ=255)
                                               tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                        >  1:1156337/1‑93 (MQ=255)
                                                aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                       >  2:2644905/1‑93 (MQ=255)
                                                aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                       >  1:1756066/1‑93 (MQ=255)
                                                  tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                     <  1:4326466/93‑1 (MQ=255)
                                                   gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                    >  1:1317585/1‑93 (MQ=255)
                                                    aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                   <  2:10282/93‑1 (MQ=255)
                                                     aaaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                  <  1:3591295/92‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 <  2:127711/93‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 >  2:648307/1‑93 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 <  1:2511149/93‑1 (MQ=255)
                                                      aaGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAAAATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCt                                 >  1:2430788/1‑93 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                <  2:1658021/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                <  1:2398592/93‑1 (MQ=255)
                                                       agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                >  2:2089836/1‑93 (MQ=255)
                                                         agaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTCTCTGTTCACCATATTGCTGgc                              >  1:1752553/1‑93 (MQ=255)
                                                          gaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcg                             >  1:2826875/1‑93 (MQ=255)
                                                            aaCCAAACCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                           <  2:3351765/93‑1 (MQ=255)
                                                            aaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCt                           <  2:3595843/93‑1 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCta                          >  1:1018494/1‑92 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                          >  2:1566675/1‑93 (MQ=255)
                                                             aCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTg                          >  1:1018493/1‑93 (MQ=255)
                                                               cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                        <  1:783479/93‑1 (MQ=255)
                                                                aaaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAt                       <  2:4380332/93‑1 (MQ=255)
                                                                 aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                      <  1:914275/93‑1 (MQ=255)
                                                                  aaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTc                     >  1:4309892/1‑93 (MQ=255)
                                                                   aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                    >  2:3826033/1‑93 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCGGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                   <  1:1337102/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACTATATTGCTGGCGCTGAATTccc                   <  2:2990198/93‑1 (MQ=255)
                                                                    ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGAATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                   >  2:2567723/1‑93 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 >  2:645055/1‑93 (MQ=255)
                                                                      gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 <  1:2216718/93‑1 (MQ=255)
                                                                      cccgccCGGTGTACCAATCCGCCCTAGCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                 <  1:1070634/92‑1 (MQ=255)
                                                                       ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                >  1:1625303/1‑93 (MQ=255)
                                                                       ccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAg                <  1:3658545/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  1:1765622/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  2:475953/93‑1 (MQ=255)
                                                                         gccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGcc              <  2:475951/93‑1 (MQ=255)
                                                                          cccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCt             >  1:4071992/1‑93 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg            >  1:2111999/1‑93 (MQ=255)
                                                                           ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg            <  1:1497715/93‑1 (MQ=255)
                                                                            cGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc           >  2:4006645/1‑93 (MQ=255)
                                                                             ggTGTCCCAATCCTCCCTAACATCGCCGCCAGCGTCACCAGCATCAAGTGTTTCTGTTCCCCATATTGCTGGCGCTGAATTCCCCAGCCTGcc          <  2:3709513/93‑1 (MQ=25)
                                                                                 tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa      <  1:3122384/93‑1 (MQ=255)
                                                                                   gcAATCCGCCCTAACATCCGCGCCAGCGTACCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    <  1:1564988/92‑1 (MQ=25)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    >  1:3118841/1‑93 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    >  2:1423065/1‑93 (MQ=255)
                                                                                   ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc    <  1:1005745/93‑1 (MQ=255)
                                                                                    cAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa   >  2:4051882/1‑93 (MQ=255)
                                                                                    cAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCa   <  1:2531912/93‑1 (MQ=255)
                                                                                     aaTCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt  >  2:3998792/1‑93 (MQ=255)
                                                                                     aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt  <  2:2770844/93‑1 (MQ=255)
                                                                                         |                                                                                        
TTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAT  >  NC_000913/1953813‑1953990

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