Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,888,301 G→T 5.3% L4L (CTC→CTA cysH ← phosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,888,3010GT 5.3% 614.8 / 18.5 169L4L (CTC→CTAcysHphosphoadenosine phosphosulfate reductase; PAPS reductase, thioredoxin dependent
Reads supporting (aligned to +/- strand):  ref base G (83/77);  new base T (5/4);  total (88/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.22e-01

Only 100 of 174 total aligned reads displayed.

GGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAA  >  NC_000913/2888157‑2888449
                                                                                                                                                |                                                                                                                                                    
    gaaAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGa                                                                                                                                              <  2:367150/149‑1 (MQ=255)
     aaaGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGAc                                                                                                                                             >  4:170302/1‑149 (MQ=255)
          aCATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGa                                                                                                                                        >  4:609766/1‑149 (MQ=255)
           cATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGAt                                                                                                                                       >  3:1450443/1‑149 (MQ=255)
           cATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGAt                                                                                                                                       <  3:1090757/149‑1 (MQ=255)
            atatTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGa                                                                                                                                              >  1:1113711/1‑141 (MQ=255)
              atTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTcc                                                                                                                                    >  3:843679/1‑149 (MQ=255)
                  aCCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCAc                                                                                                                                >  1:509843/1‑149 (MQ=255)
                   ccGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACt                                                                                                                               >  1:66139/1‑149 (MQ=255)
                       gCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGttt                                                                                                                           <  1:423458/149‑1 (MQ=255)
                         aGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgc                                                                                                                         <  3:170302/149‑1 (MQ=255)
                           aTTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGtt                                                                                                                                                  >  3:1836543/1‑122 (MQ=255)
                           aTTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcct                                                                                                                       <  2:66139/149‑1 (MQ=255)
                           aTTATCCAACGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcct                                                                                                                       >  4:969169/1‑149 (MQ=255)
                              aTCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcct                                                                                                                       <  3:1130729/146‑1 (MQ=255)
                              aTCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcct                                                                                                                       >  4:1130729/1‑146 (MQ=255)
                               tCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGttt                                                                                                                           <  1:1133026/141‑1 (MQ=255)
                               tCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGttt                                                                                                                           >  2:1133026/1‑141 (MQ=255)
                                 cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGa                                                                                                                                              >  4:1318079/1‑120 (MQ=255)
                                 cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGa                                                                                                                                              <  3:1318079/120‑1 (MQ=255)
                                 cAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTg                                                                                                                 >  4:1246357/1‑149 (MQ=255)
                                  aGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCTAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGtt                                                                                                                            <  4:1179502/137‑1 (MQ=255)
                                   gcgcCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGACCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               <  2:1109725/149‑1 (MQ=255)
                                     gcCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGc                                                                                                             <  3:969169/149‑1 (MQ=255)
                                        cAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               >  3:1288046/1‑144 (MQ=255)
                                             tACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgc                                                                                                                         <  4:1011524/129‑1 (MQ=255)
                                                     cTTCAGCGTCAAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGAAATAGATTCCTCACTGTTTGCCTTGCCTGATGCGa                                                                                                           <  2:651231/135‑1 (MQ=255)
                                                      ttCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgccttgcc                                                                                                                   <  1:510847/126‑1 (MQ=255)
                                                      ttCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgccttgc                                                                                                                    <  1:651771/125‑1 (MQ=255)
                                                      ttCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               >  2:1156142/1‑130 (MQ=255)
                                                      ttCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               <  1:1156142/130‑1 (MQ=255)
                                                      ttCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCtt                                                                                            >  3:60323/1‑149 (MQ=255)
                                                          gCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCAcgc                                                                                                  <  1:1261891/139‑1 (MQ=255)
                                                          gCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCAcgc                                                                                                  >  2:1261891/1‑139 (MQ=255)
                                                             tCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                >  1:1344152/1‑90 (MQ=255)
                                                             tCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGc                                                                                     >  2:519455/1‑149 (MQ=255)
                                                               cAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGa                                                                                                                <  4:1883562/120‑1 (MQ=255)
                                                                 gTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTAc                                                                                 <  3:1246357/149‑1 (MQ=255)
                                                                      tCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATa                                                                                                                                          <  3:1427787/87‑1 (MQ=255)
                                                                       ccAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGc                                                                                                        >  4:283449/1‑120 (MQ=255)
                                                                       ccAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGc                                                                                                        <  3:283449/120‑1 (MQ=255)
                                                                          gTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCggg                                                                        <  4:1438576/149‑1 (MQ=255)
                                                                           ttCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACaa                                                                               >  1:1408213/1‑141 (MQ=255)
                                                                           ttCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACaa                                                                               <  2:1408213/141‑1 (MQ=255)
                                                                           ttCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCt                                                                                             <  4:723099/127‑1 (MQ=255)
                                                                               gCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                <  1:1823980/72‑1 (MQ=255)
                                                                               gCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                <  3:552796/72‑1 (MQ=255)
                                                                               gCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                >  4:552796/1‑72 (MQ=255)
                                                                               gCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTg                                                                                                                                                >  2:1823980/1‑72 (MQ=255)
                                                                                cGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  <  1:907103/149‑1 (MQ=255)
                                                                                 gTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGa                                                                    <  2:942121/146‑1 (MQ=255)
                                                                                  ttAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATgg                                                                 >  3:1809060/1‑148 (MQ=255)
                                                                                  ttAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATgg                                                                 <  4:1809060/148‑1 (MQ=255)
                                                                                  ttAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGt                                                                <  3:691230/149‑1 (MQ=255)
                                                                                  ttAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGt                                                                <  1:675297/149‑1 (MQ=255)
                                                                                      tttCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAAt                                                            >  1:337207/1‑149 (MQ=255)
                                                                                      tttCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAAt                                                            >  2:840447/1‑149 (MQ=255)
                                                                                         ccGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               <  1:63282/95‑1 (MQ=255)
                                                                                         ccGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGAt                                                                                                               >  2:63282/1‑95 (MQ=255)
                                                                                          cgccagcgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCca                                                        <  2:509843/149‑1 (MQ=255)
                                                                                          cgccagcgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCca                                                        <  2:699054/149‑1 (MQ=255)
                                                                                           gccagcgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  <  2:1332627/138‑1 (MQ=255)
                                                                                              agcgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACaaa                                                    <  1:840447/149‑1 (MQ=255)
                                                                                                cgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  <  4:1606834/133‑1 (MQ=255)
                                                                                                cgccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  >  3:1606834/1‑133 (MQ=255)
                                                                                                  ccagAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCATGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCAc                                                <  3:1717591/149‑1 (MQ=255)
                                                                                                     gAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCAcgcg                                             >  1:1734133/1‑149 (MQ=255)
                                                                                                       aTGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATCGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  <  3:111161/126‑1 (MQ=255)
                                                                                                       aTGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATg                                                                  >  4:111161/1‑126 (MQ=255)
                                                                                                       aTGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATc                                                           <  3:1170515/133‑1 (MQ=255)
                                                                                                            aTCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATc                                      >  4:1104698/1‑149 (MQ=255)
                                                                                                               tACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGt                                                                >  2:1693921/1‑120 (MQ=255)
                                                                                                               tACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGt                                                                >  1:521502/1‑120 (MQ=255)
                                                                                                               tACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGt                                                                <  2:257642/120‑1 (MQ=255)
                                                                                                                aCCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGc                                  >  2:90635/1‑149 (MQ=255)
                                                                                                                 ccTTCGGCAGTTCGTTCAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCa                                 >  2:1443724/1‑149 (MQ=255)
                                                                                                                  cTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcct                                                                                                                       <  1:194505/62‑1 (MQ=255)
                                                                                                                  cTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCAc                                >  2:1592394/1‑149 (MQ=255)
                                                                                                                   ttCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACGTAGATTCCTCACTGTttgc                                                                                                                         >  1:1043587/1‑59 (MQ=255)
                                                                                                                   ttCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACGTAGATTCCTCACTGTttgc                                                                                                                         <  2:1043587/59‑1 (MQ=255)
                                                                                                                   ttCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCa                                                 >  4:1714287/1‑131 (MQ=255)
                                                                                                                         aGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcctt                                                                                                                      <  2:1102767/56‑1 (MQ=255)
                                                                                                                         aGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTttgcctt                                                                                                                      >  1:1102767/1‑56 (MQ=255)
                                                                                                                         aGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCa                                                                                                     >  3:185816/1‑73 (MQ=255)
                                                                                                                          gttcgttcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTa                                                                                           >  4:1199901/1‑82 (MQ=255)
                                                                                                                          gttcgttcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTa                                                                                           <  3:1199901/82‑1 (MQ=255)
                                                                                                                           ttcgttcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAat                       >  1:1691289/1‑149 (MQ=255)
                                                                                                                            tcgttcAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatg                      >  3:793986/1‑149 (MQ=255)
                                                                                                                              gttcAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatgat                    <  4:1416393/149‑1 (MQ=255)
                                                                                                                              gttcAGGGCGTTTAGATCTAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatgat                    <  4:793986/149‑1 (MQ=255)
                                                                                                                              gttcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAatgat                    <  4:1850306/149‑1 (MQ=255)
                                                                                                                                tcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCt                                                                                             >  3:904197/1‑74 (MQ=255)
                                                                                                                                tcAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCt                                                                                             <  4:904197/74‑1 (MQ=255)
                                                                                                                                     gCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCcgca             <  1:90635/149‑1 (MQ=255)
                                                                                                                                           aGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCAcc       <  2:614435/149‑1 (MQ=255)
                                                                                                                                             aTCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTaa                                                             >  4:1378870/1‑93 (MQ=255)
                                                                                                                                             aTCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTaa                                                             <  3:1378870/93‑1 (MQ=255)
                                                                                                                                              tCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTa    >  1:1031702/1‑149 (MQ=255)
                                                                                                                                                gAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTaaa  <  2:337207/149‑1 (MQ=255)
                                                                                                                                                gAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTaaa  <  3:1104698/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGAAGAAAGTACATATTCACCGGGCAGATTATCCAGCGCCCAGGCTACGCGGCCTTCAGCGTCCAGTTTTTCCAGTTCGGCGTTAGTTTCCGCCAGCGCCAGAATGCGATCTACCTTCGGCAGTTCGTTCAGGGCGTTTAGATCGAGTTTGGACATAGATTCCTCACTGTTTGCCTTGCCTGATGCGACGCTCACGCGTTCTTATCAGGCCTACAAGACCGGGCTGATGGTTAATCCCACAAATCACGCGCCGGATCGAGCACCGGGCGAATGATGCCCGCACGCACCGTAAA  >  NC_000913/2888157‑2888449

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: