Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,514,021 C→T 5.3% L123L (TTG→TTA gph ← phosphoglycolate phosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,514,0210CT 5.3% 309.7 / 4.6 94L123L (TTG→TTAgphphosphoglycolate phosphatase
Reads supporting (aligned to +/- strand):  ref base C (47/42);  new base T (3/2);  total (50/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.31e-01

CGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGCGGTTTACCCATT  >  NC_000913/3513886‑3514165
                                                                                                                                       |                                                                                                                                                
cGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAAcgc                                                                                                                                             <  6:57889/141‑1 (MQ=255)
cGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAAcgc                                                                                                                                             >  5:57889/1‑141 (MQ=255)
cGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAAc                                                                                                                                      <  7:338521/148‑1 (MQ=255)
cGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGATGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAAc                                                                                                                                      <  8:238782/148‑1 (MQ=255)
 ggTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACg                                                                                                                                     >  7:80921/1‑148 (MQ=255)
  gTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTa                                                                                                                                   >  8:246368/1‑149 (MQ=255)
      tttgttttgCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCg                                                                                                                                <  8:291756/148‑1 (MQ=255)
             tgCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCGTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGt                                                                                                                         <  5:22605/148‑1 (MQ=255)
                 catcatCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCggg                                                                                                                    <  1:383388/149‑1 (MQ=255)
                  atcatcACCCCCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCggg                                                                                                                    >  4:32141/1‑148 (MQ=255)
                    catcaCCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCg                                                                                                                                >  8:96561/1‑134 (MQ=255)
                    catcaCCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCg                                                                                                                                <  7:96561/134‑1 (MQ=255)
                      tcaCCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACa                                                                                                                >  3:351603/1‑148 (MQ=255)
                        accaccAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACaaaa                                                                                                             <  1:177442/149‑1 (MQ=255)
                        accaccAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACaaa                                                                                                              >  6:254472/1‑148 (MQ=255)
                         ccaccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACaaaaa                                                                                                            >  4:131978/1‑149 (MQ=255)
                          caccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGTAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAAc                                                                                                           <  2:142028/149‑1 (MQ=255)
                          caccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGTAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAAc                                                                                                           >  1:142028/1‑149 (MQ=255)
                          caccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAAc                                                                                                           <  5:217974/149‑1 (MQ=255)
                          caccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAAc                                                                                                           <  4:351603/149‑1 (MQ=255)
                           accaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACg                                                                                                          >  3:270766/1‑149 (MQ=255)
                            ccaATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGt                                                                                                         >  2:42753/1‑149 (MQ=255)
                             caATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTc                                                                                                        <  7:1996/149‑1 (MQ=255)
                             caATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTc                                                                                                        >  6:125060/1‑149 (MQ=255)
                             caATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTc                                                                                                        >  5:34025/1‑149 (MQ=255)
                              aaTCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTcc                                                                                                       >  7:9355/1‑149 (MQ=255)
                              aaTCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTcc                                                                                                       <  8:310183/149‑1 (MQ=255)
                               aTCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTccc                                                                                                      <  7:391797/149‑1 (MQ=255)
                                   ccacGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCctct                                                                                                  >  8:134810/1‑149 (MQ=255)
                                    cacGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCtt                                                                                                 >  2:233209/1‑149 (MQ=255)
                                      cGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCATCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGGAACGTGCGGGAACAAAAACGTCCCCTCTTCg                                                                                               <  3:32141/149‑1 (MQ=255)
                                       gCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCgg                                                                                              <  8:80921/149‑1 (MQ=255)
                                        cTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAAc                                                                                                                 >  8:250217/1‑129 (MQ=255)
                                        cTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAAc                                                                                                                 <  7:250217/129‑1 (MQ=255)
                                         tGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGc                                                                                             >  3:175067/1‑148 (MQ=255)
                                          gAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCaa                                                                                           >  2:170865/1‑149 (MQ=255)
                                          gAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCaa                                                                                           >  6:147005/1‑149 (MQ=255)
                                              tATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTc                                                                                       <  1:170865/149‑1 (MQ=255)
                                                tttGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCg                                                                                      <  3:131978/148‑1 (MQ=255)
                                                tttGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGc                                                                                     <  6:34025/149‑1 (MQ=255)
                                                  tGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCa                                                                                   >  1:192007/1‑149 (MQ=255)
                                                   ggCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTc                                                                                                        <  4:197028/127‑1 (MQ=255)
                                                      gATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCAtag                                                                                >  3:251453/1‑148 (MQ=255)
                                                       atatCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCAtagta                                                                              <  4:175067/149‑1 (MQ=255)
                                                       atatCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGCGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCAtag                                                                                >  6:7030/1‑147 (MQ=255)
                                                          tCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGc                                                                            >  1:355499/1‑148 (MQ=255)
                                                          tCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGc                                                                            >  2:379908/1‑148 (MQ=255)
                                                          tCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGc                                                                            <  7:134810/148‑1 (MQ=255)
                                                           cTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTc                                                                                                        >  8:105375/1‑119 (MQ=255)
                                                            tAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGAt                                                                         >  1:204901/1‑149 (MQ=255)
                                                            tAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGAt                                                                         >  7:328378/1‑149 (MQ=255)
                                                               ggCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGTACCGCGCCCAACGTATCGGCAACGTGCGGGAACaaaaa                                                                                                            >  1:255223/1‑111 (MQ=255)
                                                               ggCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGTAACGCGCCCAACGTATCGGCAACGTGCGGGAACaaaaa                                                                                                            <  2:255223/111‑1 (MQ=255)
                                                                  ttCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGaa                                                                     >  3:168838/1‑147 (MQ=255)
                                                                     gAGCAGCGGCTCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGt                                                                 <  3:379846/148‑1 (MQ=255)
                                                                     gAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGCGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGt                                                                 <  5:7030/148‑1 (MQ=255)
                                                                      atcagcGGCGCGACTAACGGAGTTGGTTTGTTGCTGACCAGGCCTAGGGGCAGTCCTTTAGCCTGCAACGCGCCCAATGTATCGGCTAAGTGCGGTAACATAAACCTCGCCTCTTCGGCAACCTAGCCATAGTAGCTATCGAACAGtt                                                                <  8:348466/146‑1 (MQ=255)
                                                                      agcagaGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGACCCCTCTTCGGCAACCACGCCATAGTAGCGATCGAAccgat                                                                >  6:372758/1‑144 (MQ=255)
                                                                       gcagcGGCGCGACGAACGGCGTCGGGTTGTTGGTGCCCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTa                                                              <  2:102763/149‑1 (MQ=255)
                                                                           cggcgCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCa                                                          <  4:54342/149‑1 (MQ=255)
                                                                                 gacgaACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATAc                                                    >  2:349718/1‑149 (MQ=255)
                                                                                     aaCGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCg                                                                                               >  7:142660/1‑102 (MQ=255)
                                                                                     aaCGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCg                                                                                               <  8:142660/102‑1 (MQ=255)
                                                                                      aCGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGTAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTAcc                                               >  8:150212/1‑149 (MQ=255)
                                                                                          cGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGtt                                           <  7:107889/149‑1 (MQ=255)
                                                                                                tttGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTg                                     >  6:64048/1‑149 (MQ=255)
                                                                                                  tgttgGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGcc                                   >  4:372876/1‑149 (MQ=255)
                                                                                                      ggTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGaa                               >  2:37548/1‑149 (MQ=255)
                                                                                                      ggTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGa                                <  2:355499/148‑1 (MQ=255)
                                                                                                       gTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAAt                              >  6:102938/1‑149 (MQ=255)
                                                                                                       gTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCATCGGCAACCTAGCCATACTAGCGATCGAACAGTATACGCAGAATACATACCTTTTCATCTGCCAGAAt                              >  2:413647/1‑149 (MQ=255)
                                                                                                          aCCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgtc                           <  1:379908/149‑1 (MQ=255)
                                                                                                          aCCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgtc                           <  5:125060/149‑1 (MQ=255)
                                                                                                          aCCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgt                            <  1:42753/148‑1 (MQ=255)
                                                                                                            cAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgtcgt                         >  2:337167/1‑149 (MQ=255)
                                                                                                            cAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgtcgt                         >  8:208475/1‑149 (MQ=255)
                                                                                                            cAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATgtcgt                         <  2:204901/149‑1 (MQ=255)
                                                                                                               gCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGtcat                      >  2:60565/1‑149 (MQ=255)
                                                                                                               gCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGtcat                      >  7:267733/1‑149 (MQ=255)
                                                                                                                  tAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATg                             >  8:292413/1‑139 (MQ=255)
                                                                                                                  tAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATg                             <  7:292413/139‑1 (MQ=255)
                                                                                                                    gCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACg                 <  4:270766/149‑1 (MQ=255)
                                                                                                                      ggCAGGCCTTTAGCCTGCAACGCGCCCACCGTATCGCCAACGAGCGGCAACAAACACGTCCCCACTTCGGCAACCACGCCATAGTAGCGAACGAACAGTTTACGCAGAAAACGAACCTGTTCTTCTGCCGGAATGTCGtcatcacccgg               >  1:200729/1‑144 (MQ=255)
                                                                                                                      ggCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACggg               >  4:263374/1‑149 (MQ=255)
                                                                                                                       gCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGc              >  3:275313/1‑149 (MQ=255)
                                                                                                                         aGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGCgg            <  1:349718/149‑1 (MQ=255)
                                                                                                                           gCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGCGGtt          >  7:345557/1‑149 (MQ=255)
                                                                                                                             cTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGGCATCAACGGGCGGTTTa        <  5:254472/149‑1 (MQ=255)
                                                                                                                                tAGCCTGCAACGCGCCCTACGTATCTGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCATCCTCGCCATAGTAGCGATCGAACAGTTTTCGCTGAATACGTACCTGTTCTTCTGCCGGATTGTCGTCATCAACGGGCGGTTTAc       <  5:372758/147‑1 (MQ=255)
                                                                                                                                  gCCTGCTACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAg                                                                  <  8:111396/86‑1 (MQ=255)
                                                                                                                                  gCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAg                                                                  >  7:111396/1‑86 (MQ=255)
                                                                                                                                   ccTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGCGGTTTACCCAtt  <  5:102938/149‑1 (MQ=255)
                                                                                                                                      gCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATc                                                                        <  8:122746/76‑1 (MQ=255)
                                                                                                                                      gCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATc                                                                        >  7:122746/1‑76 (MQ=255)
                                                                                                                                       |                                                                                                                                                
CGGTTTTTTGTTTTGCACATCATCACCACCAATCACCACGCTGAAGTATTTGGCGATATCTAAGGCTTCGAGCAGCGGCGCGACGAACGGCGTCGGTTTGTTGGTGACCAGGCCTAGCGGCAGGCCTTTAGCCTGCAACGCGCCCAACGTATCGGCAACGTGCGGGAACAAAAACGTCCCCTCTTCGGCAACCTCGCCATAGTAGCGATCGAACAGTTTACGCAGAATACGTACCTGTTCTTCTGCCGGAATGTCGTCATCAACGGGCGGTTTACCCATT  >  NC_000913/3513886‑3514165

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: