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% 665.9 / ‑5.7 205Q377R (CAG→CGG) b1872biotin sulfoxide reductase 2
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (112/92);  total (112/93)
Fisher's exact test for biased strand distribution p-value = 4.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953812‑1953991
                                                                                          |                                                                                         
aTTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGc                                                                                         >  1:1754818/1‑93 (MQ=255)
aTTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGc                                                                                         >  1:2611047/1‑93 (MQ=255)
 tttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGcc                                                                                        >  1:6316671/1‑93 (MQ=255)
  ttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGccc                                                                                       <  2:2351171/93‑1 (MQ=255)
  ttCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGccc                                                                                       <  1:4454726/93‑1 (MQ=255)
   tCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCt                                                                                      <  1:1675677/93‑1 (MQ=255)
    ccGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTa                                                                                     <  2:3742429/93‑1 (MQ=255)
     cGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTaa                                                                                    <  2:2969291/93‑1 (MQ=255)
       gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  >  1:2961024/1‑93 (MQ=255)
       gcagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACa                                                                                  >  1:2758439/1‑93 (MQ=255)
        cagcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACAt                                                                                 <  2:3635665/93‑1 (MQ=255)
         agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                                >  1:704016/1‑93 (MQ=255)
         agcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATc                                                                                >  2:5643667/1‑93 (MQ=255)
          gcaCGCCGCCAACACGCGTCGGGTTGCGACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCCAACATcc                                                                               >  1:4617864/1‑92 (MQ=255)
          gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                               >  2:6694891/1‑93 (MQ=255)
          gcaCGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATcg                                                                               >  2:4102368/1‑93 (MQ=255)
               ccgccAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          >  2:3715230/1‑93 (MQ=255)
               ccgccAACACGCGTCGGGTTACCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCgccgcc                                                                          <  1:456084/93‑1 (MQ=255)
                   cAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCg                                                                      >  2:6402101/1‑93 (MQ=255)
                     acacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                    >  2:2192968/1‑93 (MQ=255)
                     acacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTc                                                                    <  2:2694745/93‑1 (MQ=255)
                      cacGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCa                                                                   >  1:2806094/1‑93 (MQ=255)
                       acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                  >  1:505114/1‑93 (MQ=255)
                       acGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCAc                                                                  <  1:5675677/93‑1 (MQ=255)
                         gcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCa                                                                <  2:2256377/93‑1 (MQ=255)
                         gcgTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCa                                                                >  2:1921967/1‑93 (MQ=255)
                             cGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCAt                                                            >  2:1237762/1‑93 (MQ=255)
                              gggTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATc                                                           >  2:1780934/1‑93 (MQ=255)
                                gTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCa                                                         >  2:3676112/1‑93 (MQ=255)
                                gTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCa                                                         <  1:5143409/93‑1 (MQ=255)
                                  tGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgt                                                       <  2:281676/93‑1 (MQ=255)
                                   gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                      >  2:5950451/1‑93 (MQ=255)
                                   gCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtg                                                      >  2:5559027/1‑93 (MQ=255)
                                    ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                     <  2:4464819/93‑1 (MQ=255)
                                    ccACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAgtgt                                                     >  2:6453339/1‑93 (MQ=255)
                                     cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                    >  1:3922138/1‑93 (MQ=255)
                                     cACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGtt                                                    <  1:2225434/93‑1 (MQ=255)
                                         gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                >  2:520090/1‑93 (MQ=255)
                                         gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                <  2:2520034/93‑1 (MQ=255)
                                         gTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTg                                                <  2:4539520/93‑1 (MQ=255)
                                           tGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                              >  2:6277524/1‑93 (MQ=255)
                                           tGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGtt                                              <  1:2107999/93‑1 (MQ=255)
                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          <  2:612846/93‑1 (MQ=255)
                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          >  2:504398/1‑93 (MQ=255)
                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          <  1:1410942/93‑1 (MQ=255)
                                               gTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCAcc                                          <  1:3011140/93‑1 (MQ=255)
                                                tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         <  2:3428956/93‑1 (MQ=255)
                                                tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  1:5336585/1‑93 (MQ=255)
                                                tAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCa                                         >  2:5086892/1‑93 (MQ=255)
                                                 aaTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCat                                        >  1:5114359/1‑93 (MQ=255)
                                                  aTGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCata                                       >  1:5224231/1‑93 (MQ=255)
                                                   tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      >  2:5301466/1‑93 (MQ=255)
                                                   tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      <  2:1499920/93‑1 (MQ=255)
                                                   tGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCatat                                      >  1:5427092/1‑93 (MQ=255)
                                                    gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                     <  2:76377/93‑1 (MQ=255)
                                                    gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                     <  1:3706194/93‑1 (MQ=255)
                                                    gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                     <  1:978414/93‑1 (MQ=255)
                                                    gATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATAtt                                     <  2:1441160/93‑1 (MQ=255)
                                                     aTAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTg                                    >  2:3926445/1‑93 (MQ=255)
                                                      tAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGc                                   >  2:3536509/1‑93 (MQ=255)
                                                        agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTACCATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  1:3261712/93‑1 (MQ=255)
                                                        agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  2:6219884/93‑1 (MQ=255)
                                                        agagAAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTg                                 <  1:2388088/93‑1 (MQ=255)
                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                >  1:1630658/1‑93 (MQ=255)
                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                >  1:3377587/1‑93 (MQ=255)
                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                <  1:858157/93‑1 (MQ=255)
                                                         gagaAACCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTgg                                >  1:5153881/1‑93 (MQ=255)
                                                            aaaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcgc                             >  1:1971863/1‑93 (MQ=255)
                                                            aaaCCAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGgcgc                             <  2:5683294/93‑1 (MQ=255)
                                                                cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                         >  2:724572/1‑93 (MQ=255)
                                                                cAAAACCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGaa                         <  2:1624635/93‑1 (MQ=255)
                                                                  aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGACtt                       >  1:1001082/1‑93 (MQ=255)
                                                                  aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                       >  1:2196526/1‑93 (MQ=255)
                                                                  aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                       >  1:2359625/1‑93 (MQ=255)
                                                                  aaaCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtt                       >  1:6153444/1‑93 (MQ=255)
                                                                    aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAAtttc                     >  1:1790549/1‑91 (MQ=255)
                                                                    aCCGCCGCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcc                     <  1:2713610/93‑1 (MQ=255)
                                                                     ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                    >  1:4351543/1‑93 (MQ=255)
                                                                     ccgccgCCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTccc                    <  2:2898162/93‑1 (MQ=255)
                                                                      cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  2:4627517/93‑1 (MQ=255)
                                                                      cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   <  2:6338356/93‑1 (MQ=255)
                                                                      cgccgcCCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTcccc                   >  1:891410/1‑93 (MQ=255)
                                                                       gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCa                  >  1:5931248/1‑93 (MQ=255)
                                                                       gccgccCGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGGGTTTCTGTTCACCATATTGTTGGCGCTGAATTCCCCa                  >  2:3262415/1‑93 (MQ=255)
                                                                            ccGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTg             >  2:3301599/1‑93 (MQ=255)
                                                                             cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTGACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  2:4897425/1‑93 (MQ=255)
                                                                             cGGTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGc            >  1:382371/1‑93 (MQ=255)
                                                                              ggTGTACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGcc           <  1:3358193/93‑1 (MQ=255)
                                                                               gtgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCa          >  1:654926/1‑93 (MQ=255)
                                                                                tgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt         >  1:1420499/1‑93 (MQ=255)
                                                                                tgtACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCAt         >  2:959634/1‑93 (MQ=255)
                                                                                 gtACCAATCCGCCCTAACATCGCCGCCAGCGTCCCCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATc        >  1:4702053/1‑93 (MQ=255)
                                                                                  tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  2:1272036/93‑1 (MQ=255)
                                                                                  tACCAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCa       <  1:3733623/93‑1 (MQ=255)
                                                                                    ccAATCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGc     <  1:4670070/93‑1 (MQ=255)
                                                                                      aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   <  1:191449/93‑1 (MQ=255)
                                                                                      aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAt   <  2:4565851/93‑1 (MQ=255)
                                                                                      aaTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCAGCCTGCCATCAGCAtt  <  2:1097213/93‑1 (MQ=255)
                                                                                       aTCCGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCAtt  <  2:4704852/93‑1 (MQ=255)
                                                                                       aTCCGCCCTAAAGTTGCCGGCAGCGTCTTCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCCCTGAATTCCCCAGCCTGCCATCAGCAtt  <  2:5315329/93‑1 (MQ=21)
                                                                                          |                                                                                         
ATTTCCGGCAGCACGCCGCCAACACGCGTCGGGTTGCCACCGTTGGAGTAATGATAAGAGAAACCAAAACCGCCGCCCGGTGTACCAATCTGCCCTAACATCGCCGCCAGCGTCACCAGCATCCAGTGTTTCTGTTCACCATATTGCTGGCGCTGAATTCCCCAGCCTGCCATCAGCATT  >  NC_000913/1953812‑1953991

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