Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,422,372 G→A 5.8% A63T (GCC→ACC)  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,3720GA 5.8% 228.4 / 2.6 69A63T (GCC→ACC) ulaD3‑keto‑L‑gulonate 6‑phosphate decarboxylase
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base A (2/2);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

TATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTG  >  NC_000913/4422227‑4422512
                                                                                                                                                 |                                                                                                                                            
tATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGcc                                                                                                                                            >  3:87968/1‑148 (MQ=255)
  tGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCa                                                                                                                                           >  8:100710/1‑147 (MQ=255)
      tAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAAtt                                                                                                                                      >  3:198432/1‑148 (MQ=255)
                           ccTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGAAGCCAAAATTGCCGATGCAGGCAAAATCCtt                                                                                                                 >  8:375202/1‑148 (MQ=255)
                                            tCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGa                                                                                                <  7:375202/148‑1 (MQ=255)
                                             agACATTATCGAAGTGGGCACGATTCTGTGCGTGGGCGAAGTCGTGCGTTCGGTTCGTGACCTGAAAGCGCTCTACCCGCCCAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCATTTCGCGTATGTGCTTCGat                                                                                               >  8:163348/2‑147 (MQ=255)
                                               aCATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCtt                                                                                                                 <  8:332568/128‑1 (MQ=255)
                                               aCATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCtt                                                                                                                 >  7:332568/1‑128 (MQ=255)
                                               aCATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTAtgt                                                                                                       <  4:332012/138‑1 (MQ=255)
                                               aCATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTAtgt                                                                                                       >  3:332012/1‑138 (MQ=255)
                                                   tATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         >  4:321908/1‑148 (MQ=255)
                                                   tATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGAGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATAGTACTGGCAGCCGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         >  7:116255/1‑148 (MQ=255)
                                                      cGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGCCCTGAAAGCGCTCTACCCGCACAAAATCGTACTCGCAGACGCCCAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCt                                                                                      >  5:278413/1‑148 (MQ=255)
                                                       gAAGTGTGCAGCATTCTGTGCGTGGGCGAAGTGGTGCGTGCGGTGCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTTGCAGACGCCAAAATTGCCGATGCAGGAAAAATCCTTTCTCGTATGTGCTTCGAAGCGAACGCTg                                                                                     <  3:321908/148‑1 (MQ=255)
                                                       gAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTg                                                                                     <  1:21811/148‑1 (MQ=255)
                                                       gAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTg                                                                                     <  3:342548/148‑1 (MQ=255)
                                                         aGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGAc                                                                                   <  3:136080/148‑1 (MQ=255)
                                                            gggCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTgg                                                                                >  7:59543/1‑148 (MQ=255)
                                                            gggCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTgg                                                                                >  5:236340/1‑148 (MQ=255)
                                                             ggCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTgg                                                                                <  7:100710/147‑1 (MQ=255)
                                                                  cATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGa                                                                                                <  7:128212/126‑1 (MQ=255)
                                                                  cATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGa                                                                                                >  8:128212/1‑126 (MQ=255)
                                                                  cATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGGCg                                                                          >  3:147457/1‑148 (MQ=255)
                                                                   aTTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACACCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTg                                                                                     <  7:128515/136‑1 (MQ=255)
                                                                   aTTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACACCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTg                                                                                     >  8:128515/1‑136 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCa                                                                                           >  3:173193/1‑127 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCa                                                                                           <  4:173193/127‑1 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         <  5:286223/129‑1 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         >  6:291327/1‑129 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         >  6:286223/1‑129 (MQ=255)
                                                                      cTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         <  5:291327/129‑1 (MQ=255)
                                                                        gtgCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGa                                                                                    <  8:123408/132‑1 (MQ=255)
                                                                        gtgCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGa                                                                                    >  7:123408/1‑132 (MQ=255)
                                                                         tgCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAAttt                                                                   >  4:236369/1‑148 (MQ=255)
                                                                          gCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTg                                                                  >  3:181044/1‑148 (MQ=255)
                                                                           cGTGGGCGAAGGCGTGCGTGCGGGTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCTAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGc                                                                 <  1:274477/148‑1 (MQ=255)
                                                                            gTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCt                                                                <  6:236340/148‑1 (MQ=255)
                                                                              gggCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGt                                                                                                           >  2:70455/1‑103 (MQ=255)
                                                                              gggCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGt                                                                                                           <  1:70455/103‑1 (MQ=255)
                                                                                   aaGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGa                                                         <  3:236369/148‑1 (MQ=255)
                                                                                     ggCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGata                                                       >  6:31791/1‑148 (MQ=255)
                                                                                        gtgcgtgcGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCa                                                    <  4:181044/148‑1 (MQ=255)
                                                                                          gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCg                                                                                                                                  >  2:337583/1‑68 (MQ=255)
                                                                                          gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCg                                                                                                                                  <  1:337583/68‑1 (MQ=255)
                                                                                          gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGaa                                                                                               >  7:344401/1‑103 (MQ=255)
                                                                                          gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGaa                                                                                               <  8:344401/103‑1 (MQ=255)
                                                                                          gcgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCaa                                                   >  3:348646/1‑147 (MQ=255)
                                                                                           cgtgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAac                                                  <  5:31791/147‑1 (MQ=255)
                                                                                             tgcgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGGCGGTAATTTGCTGTGCGGATATCAACAcc                                               <  4:147457/148‑1 (MQ=255)
                                                                                               cgGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCg                                              <  4:198432/147‑1 (MQ=255)
                                                                                                ggTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGcc                                            <  8:59543/148‑1 (MQ=255)
                                                                                                       gtccTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGgcg                                     >  5:49222/3‑148 (MQ=255)
                                                                                                        aCCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGgcgc                                    >  5:160211/1‑148 (MQ=255)
                                                                                                            gAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         <  7:243445/91‑1 (MQ=255)
                                                                                                            gAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAAc                                                                                         >  8:243445/1‑91 (MQ=255)
                                                                                                            gAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTg                                >  7:162778/1‑148 (MQ=255)
                                                                                                            gAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTg                                >  4:84722/1‑148 (MQ=255)
                                                                                                               cgcgcTCTACCGGCACAAAATCGTACTGGAAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGGTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGTTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGAc                             <  6:278413/147‑1 (MQ=255)
                                                                                                                     ctACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGc                        >  2:384587/1‑147 (MQ=255)
                                                                                                                               aaaaTCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTa             <  1:384587/148‑1 (MQ=255)
                                                                                                                                 aaTCGTACTGGCAGACACCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCaaa                     <  1:339957/138‑1 (MQ=255)
                                                                                                                                 aaTCGTACTGGCAGACACCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCaaa                     >  2:339957/1‑138 (MQ=255)
                                                                                                                                      tACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGa      <  6:49222/148‑1 (MQ=255)
                                                                                                                                          ggCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTa             >  3:326379/1‑137 (MQ=255)
                                                                                                                                          ggCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTa             <  4:326379/137‑1 (MQ=255)
                                                                                                                                          ggCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTg  >  2:122029/1‑148 (MQ=255)
                                                                                                                                          ggCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTg  >  5:193199/1‑148 (MQ=255)
                                                                                                                                                cGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAag                   <  3:107380/125‑1 (MQ=255)
                                                                                                                                                cGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAag                   >  4:107380/1‑125 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
TATGGATAGCGCCTACGAAACCACTCGCCTGATTGCCGAAGAAGTCGACATTATCGAAGTGGGCACCATTCTGTGCGTGGGCGAAGGCGTGCGTGCGGTTCGTGACCTGAAAGCGCTCTACCCGCACAAAATCGTACTGGCAGACGCCAAAATTGCCGATGCAGGCAAAATCCTTTCGCGTATGTGCTTCGAAGCCAACGCTGACTGGGTGACGGTAATTTGCTGTGCGGATATCAACACCGCCAAAGGCGCGCTGGACGTGGCAAAAGAGTTTAACGGCGACGTG  >  NC_000913/4422227‑4422512

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