Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 909,422 G→A 41.8% A543V (GCC→GTC)  poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913909,4220GA41.8% 43.7 / 75.3 67A543V (GCC→GTC) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base G (20/19);  new base A (14/14);  total (34/33)
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

CTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGC  >  NC_000913/909278‑909569
                                                                                                                                                |                                                                                                                                                   
cTGCGTAAATCAATCATGGCATGTCCTTATTATGACCGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTAATCACCGAGTCCGCTGATGATTGCGCGCAGCATATACAGTCTGAAAACTTTGACCt                                                                                                                                                  >  4:88151/1‑148 (MQ=255)
         tcaatcaTTGCATGTCCTTATTATGACGTGAAATGCCTCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGtt                                                                                                                                         <  6:257164/148‑1 (MQ=255)
         tcaatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTATTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGtt                                                                                                                                         >  4:90398/1‑148 (MQ=255)
          caatcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGttt                                                                                                                                        >  7:331629/1‑148 (MQ=255)
            atcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGa                                                                                                                                      >  7:192816/1‑148 (MQ=255)
            atcaTGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGa                                                                                                                                      <  7:105649/148‑1 (MQ=255)
                  gCATGTCCTTATTATGACGGGAACTGCCACCCTTTTTACATTAGCCAGTTTGTTTTAACCAGTTCGATCACTTCATCACCGCTTCAGCTGATGATTGCGCGCAGCATATACAGGATGAAACCTTTGGCCTGTTCGAGTTTGATCTGct                                                                                                                                >  6:303236/1‑147 (MQ=255)
                  gCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGc                                                                                                                                 <  7:318229/147‑1 (MQ=255)
                   cATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCgg                                                                                                                               <  1:233104/148‑1 (MQ=255)
                    aTGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGt                                                                                                                              <  8:331629/148‑1 (MQ=255)
                      gTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTg                                                                                                                             >  1:281008/1‑147 (MQ=255)
                       tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           <  7:149434/148‑1 (MQ=255)
                       tCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           <  4:105612/148‑1 (MQ=255)
                         cTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGaa                                                                                                                          <  8:270909/147‑1 (MQ=255)
                         cTTAATATGACTCGAAATGAAACACTTTTTACCTTAGCCAGTTTGTTTTCGCCATTTCTATCACTTCAGCACCGCGTCATCTGATGATTGAGCGCAGCATATACAGGATGAAACCATTGACCTGTTCGAGTTTGATATGAGGTGTaa                                                                                                                          <  5:61509/147‑1 (MQ=255)
                          ttattaTGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           <  3:91098/145‑1 (MQ=255)
                          ttattaTGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTACATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           >  4:91098/1‑145 (MQ=255)
                           tattatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTg                                                                                                                                       <  2:158218/132‑1 (MQ=255)
                           tattatGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTg                                                                                                                                       >  1:158218/1‑132 (MQ=255)
                               atGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTa                                                                                                                    <  2:281008/147‑1 (MQ=255)
                                  aCGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAActc                                                                                                                <  4:56182/148‑1 (MQ=255)
                                            cgaCCCTTTTTACCTTAGCCAGTTTGTTATCACGAGGTCGAGCCCTTCATCCCCGCGTCCGCTGATGATTGCGCGCAGCATATACTGGCTGATACCTTTGGCCTGTTCGTGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTgg                                                                                                        <  8:282818/144‑1 (MQ=255)
                                              aCCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGac                                                                                                    >  1:250077/1‑148 (MQ=255)
                                                ccTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           >  2:198160/1‑123 (MQ=255)
                                                ccTTTTTACCTTAGCCAGTTTGTGTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGa                                                                                                                           <  1:198160/123‑1 (MQ=255)
                                                             gCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATccacca                                                                                      <  3:90398/147‑1 (MQ=255)
                                                                aGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAacac                                                                                  >  7:147729/1‑148 (MQ=255)
                                                                  tttgtttTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCg                                                                                >  6:109687/1‑148 (MQ=255)
                                                                   ttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacc                                                                                                   <  2:192159/128‑1 (MQ=255)
                                                                   ttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGacc                                                                                                   >  1:192159/1‑128 (MQ=255)
                                                                   ttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCg                                                                                                                                <  2:186955/99‑1 (MQ=255)
                                                                   ttgtttTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCg                                                                                                                                >  1:186955/1‑99 (MQ=255)
                                                                   ttgtttTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTATGGCGACCACCACATCCACCAACACCg                                                                                >  1:162998/1‑147 (MQ=255)
                                                                    tgtttTCACCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTATTCTTTGTCGACCACCACATCCACCAACACCgg                                                                               <  8:84853/147‑1 (MQ=255)
                                                                        ttCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGt                                                                          >  5:210373/1‑148 (MQ=255)
                                                                          caccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCttt                                                                                                          >  8:300068/3‑114 (MQ=255)
                                                                          caccAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCttt                                                                                                          <  7:300068/112‑1 (MQ=255)
                                                                           gCCAGTTCGATCACTTCATCACCGCGGCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGa                                                                       <  5:97404/148‑1 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                                                              >  8:216020/1‑148 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                                                              >  8:206194/1‑148 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                                                              <  8:147729/148‑1 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgag                                                                  >  2:199804/1‑144 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGgag                                                                  <  1:199804/144‑1 (MQ=255)
                                                                                    aTCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAgg                                                              >  8:174520/1‑148 (MQ=255)
                                                                                     tCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGgc                                                             >  8:328645/1‑148 (MQ=255)
                                                                                       aCTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGTTTTTTTTTTATTTTTCGTCGACCACATAGATCAACACCGGACCGTCGATGGAGAAGgcgc                                                           <  8:192816/148‑1 (MQ=255)
                                                                                         ttCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATgg                                                                    <  5:290557/137‑1 (MQ=255)
                                                                                         ttCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATgg                                                                    >  6:290557/1‑137 (MQ=255)
                                                                                             tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACcaca                                                                                             >  1:179117/1‑108 (MQ=255)
                                                                                             tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGttt                                                                                                                                        >  8:58173/1‑65 (MQ=255)
                                                                                             tcaCCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGttt                                                                                                                                        <  7:58173/65‑1 (MQ=255)
                                                                                                            atgatTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGTGGTTCATCAACtt                                      >  6:252334/1‑148 (MQ=255)
                                                                                                             tgatTGCGCGCAGCATATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACtt                                      >  8:239551/1‑147 (MQ=255)
                                                                                                                    gagcagcaTATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGAGTTGCAGGGCTTCATCAACTTCAGACGCt                              <  3:88151/146‑1 (MQ=255)
                                                                                                                          caTATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACtt                                      <  6:171369/134‑1 (MQ=255)
                                                                                                                          caTATACAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACtt                                      >  5:171369/1‑134 (MQ=255)
                                                                                                                             ataCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                                                                             >  2:54415/1‑92 (MQ=255)
                                                                                                                             ataCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGAc                                                                             <  1:54415/92‑1 (MQ=255)
                                                                                                                             ataCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                     >  2:66521/1‑148 (MQ=255)
                                                                                                                              taCAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                    >  3:196397/1‑148 (MQ=255)
                                                                                                                              taCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTacac                    >  2:73915/1‑148 (MQ=255)
                                                                                                                              taCAGGCTGAAACCTTTGACCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTaca                     <  7:174520/147‑1 (MQ=255)
                                                                                                                                cAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACgg                  <  7:216020/148‑1 (MQ=255)
                                                                                                                                       aaaCCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg                                                       >  5:258791/1‑104 (MQ=255)
                                                                                                                                       aaaCCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTg                                                       <  6:258791/104‑1 (MQ=255)
                                                                                                                                         aCCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTa         <  2:162998/148‑1 (MQ=255)
                                                                                                                                         aCCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATAAACTTCAGACGCTTTTTCTACACGGATACCCGt          >  2:306019/1‑147 (MQ=255)
                                                                                                                                                gCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGc  <  7:328645/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
CTGCGTAAATCAATCATGGCATGTCCTTATTATGACGGGAAATGCCACCCTTTTTACCTTAGCCAGTTTGTTTTCGCCAGTTCGATCACTTCATCACCGCGTCCGCTGATGATTGCGCGCAGCATATACAGGCTGAAACCTTTGGCCTGTTCGAGTTTGATCTGCGGTGGAATGGCTAACTCTTCTTTGGCGACCACCACATCCACCAACACCGGACCGTCGATGGAGAAGGCGCGTTGCAGGGCTTCATCAACTTCAGACGCTTTTTCTACACGGATACCCGTAATGCCGC  >  NC_000913/909278‑909569

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: