Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,422,793 C→T A203V (GCC→GTC)  ulaD → 3‑keto‑L‑gulonate 6‑phosphate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,422,7930CT98.8% 276.4 / NA 86A203V (GCC→GTC) ulaD3‑keto‑L‑gulonate 6‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (45/40);  minor base G (1/0);  total (46/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGG  >  NC_000913/4422673‑4422927
                                                                                                                        |                                                                                                                                      
gCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGCAGTCGCACGTCAGTTCAAACg                                                                                                                      >  2:423198/1‑139 (MQ=255)
    cAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTcc                                                                                                                  >  2:650617/1‑139 (MQ=255)
    cAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTcc                                                                                                                  >  1:1166594/1‑139 (MQ=255)
           accgtcaccgGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTg                                                                                                           >  1:1167844/1‑139 (MQ=255)
                caccgGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                                      <  1:885839/139‑1 (MQ=255)
                 accgGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCaaa                                                                                                                        <  2:460097/120‑1 (MQ=255)
                 accgGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCaaa                                                                                                                        >  1:460097/1‑120 (MQ=255)
                       ggCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCt                                                                                               <  1:650617/139‑1 (MQ=255)
                            ggCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGa                                                                                          >  1:550572/1‑139 (MQ=255)
                             gcgcTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAg                                                                                         <  2:550572/139‑1 (MQ=255)
                                          tGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCt                                                                                               >  1:1563620/1‑120 (MQ=255)
                                          tGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCt                                                                                               <  2:1563620/120‑1 (MQ=255)
                                                  ttCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCcaaacaaa                                                                    >  1:380549/1‑139 (MQ=255)
                                                        ggTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGc                                                              <  2:380549/139‑1 (MQ=255)
                                                          tATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCg                                                                                       <  2:1205121/112‑1 (MQ=255)
                                                          tATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCg                                                                                       >  1:1205121/1‑112 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGt                                                                                                                                 >  2:501522/1‑68 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGt                                                                                                                                 <  1:501522/68‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCaaa                                                                                                                        >  1:734001/1‑77 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCaaa                                                                                                                        <  2:734001/77‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCTAACGTTCCATCGCTGAACtgtg                                                                                                    <  1:572395/97‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCTTCGCTGCACTATGGGGCTAAGGAGCGGATATGTTGTACAAACAAATCTCGCTTgg                                                          >  2:174058/1‑139 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtgtg                                                                                                    >  2:572395/1‑97 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtgt                                                                                <  1:845874/117‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtgt                                                                                >  2:845874/1‑117 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtg                                                                                 <  2:604394/116‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATAtg                                                                                 >  1:604394/1‑116 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCca                                                                          <  1:195310/123‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCca                                                                          >  2:195310/1‑123 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTgg                                                          <  2:1167844/139‑1 (MQ=255)
                                                            ttCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTgg                                                          >  1:81772/1‑139 (MQ=255)
                                                             tCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                         <  1:730533/139‑1 (MQ=255)
                                                                gATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATc                                                      <  1:414271/139‑1 (MQ=255)
                                                                gATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATc                                                      <  2:1593911/139‑1 (MQ=255)
                                                                 aTTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         >  2:1622787/1‑87 (MQ=255)
                                                                 aTTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         <  1:1622787/87‑1 (MQ=255)
                                                                  ttCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTa                                                    >  2:810868/1‑139 (MQ=255)
                                                                   tCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                         <  2:1097727/133‑1 (MQ=255)
                                                                   tCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGATTGGc                                                         >  1:1097727/1‑133 (MQ=255)
                                                                    cACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCt                                                                                                            <  1:748290/81‑1 (MQ=255)
                                                                    cACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCt                                                                                                            >  2:748290/1‑81 (MQ=255)
                                                                      cGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGaa                                                >  2:1288564/1‑139 (MQ=255)
                                                                       gTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGaaa                                               >  2:57840/1‑139 (MQ=255)
                                                                       gTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGaaa                                               >  2:846091/1‑139 (MQ=255)
                                                                             aTCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAgcgc                                         <  1:57840/139‑1 (MQ=255)
                                                                             aTCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAgcgc                                         <  1:846091/139‑1 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                                   <  2:1554638/59‑1 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                                   <  1:1631542/59‑1 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                                   >  2:1631542/1‑59 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         >  2:330257/1‑69 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         <  1:330257/69‑1 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccg                                   >  2:1128421/1‑139 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTATCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCccg                                   >  1:860265/1‑139 (MQ=255)
                                                                                   ggCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGGCGCACGTCAGTTCAAACGTTc                                                                                                                   >  1:1554638/1‑59 (MQ=255)
                                                                                        tagtaTCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCggg                              >  1:1366874/1‑139 (MQ=255)
                                                                                             tCCGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTg                         <  1:1288564/139‑1 (MQ=255)
                                                                                              ccGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAg                                                                                         <  2:1665416/74‑1 (MQ=255)
                                                                                              ccGTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAg                                                                                         >  1:1665416/1‑74 (MQ=255)
                                                                                                gTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                         >  2:103551/1‑104 (MQ=255)
                                                                                                gTGATGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGc                                                         <  1:103551/104‑1 (MQ=255)
                                                                                                   aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                                      <  1:525943/56‑1 (MQ=255)
                                                                                                   aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                                      >  2:525943/1‑56 (MQ=255)
                                                                                                   aTGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgctggct                   >  2:364056/1‑139 (MQ=255)
                                                                                                    tGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTg                                                                                                           <  1:263314/50‑1 (MQ=255)
                                                                                                    tGCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTg                                                                                                           >  2:263314/1‑50 (MQ=255)
                                                                                                     gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCt                                                                                               <  2:896836/61‑1 (MQ=255)
                                                                                                     gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCt                                                                                               >  1:896836/1‑61 (MQ=255)
                                                                                                     gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc                        <  1:837078/132‑1 (MQ=255)
                                                                                                     gCCGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTgc                        >  2:837078/1‑132 (MQ=255)
                                                                                                      ccGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGa                <  1:174058/139‑1 (MQ=255)
                                                                                                       cGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                                   <  2:11499/39‑1 (MQ=255)
                                                                                                       cGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTc                                                                                                                   >  1:11499/1‑39 (MQ=255)
                                                                                                       cGCTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGaa               <  2:1166594/139‑1 (MQ=255)
                                                                                                         cTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAg                                                                                         <  2:1321688/63‑1 (MQ=255)
                                                                                                         cTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAg                                                                                         >  1:1321688/1‑63 (MQ=255)
                                                                                                         cTTCTCCGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACg             >  2:222821/1‑139 (MQ=255)
                                                                                                              ccGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGc        >  2:1498504/1‑139 (MQ=255)
                                                                                                               cGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         <  2:1219635/41‑1 (MQ=255)
                                                                                                               cGGTGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGaa                                                                                                         >  1:1219635/1‑41 (MQ=255)
                                                                                                                ggTGGAAGTCGCACGTCAGTTAAAACGTTCCATCGCTGAACTGTGGGGCTa                                                                                              <  1:1130682/51‑1 (MQ=255)
                                                                                                                ggTGGAAGTCGCACGACAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTa                                                                                              >  2:1130682/1‑51 (MQ=255)
                                                                                                                  tGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                                      >  1:159935/1‑41 (MQ=255)
                                                                                                                  tGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACtg                                                                                                      <  2:159935/41‑1 (MQ=255)
                                                                                                                  tGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGaa               <  2:1607044/128‑1 (MQ=255)
                                                                                                                  tGGAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGaa               >  1:1607044/1‑128 (MQ=255)
                                                                                                                    gAAGTCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTgg  >  1:260710/1‑139 (MQ=255)
                                                                                                                        |                                                                                                                                      
GCTTCAAAGTCACCGTCACCGGAGGCCTGGCGCTGGAAGATCTGCCGCTGTTCAAGGGTATTCCGATTCACGTCTTTATCGCGGGCCGTAGTATCCGTGATGCCGCTTCTCCGGTGGAAGCCGCACGTCAGTTCAAACGTTCCATCGCTGAACTGTGGGGCTAAGGAGCGGATATGTTGTCCAAACAAATCCCGCTTGGCATCTATGAAAAAGCGCTCCCCGCCGGGGAGTGCTGGCTGGAACGCCTGCAACTGG  >  NC_000913/4422673‑4422927

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