Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,222 Δ1 bp 14.7% intergenic (+30/‑111) udp → / → rmuC uridine phosphorylase/DNA recombination protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,017,2220T.14.7% 217.6 / 31.9 68intergenic (+30/‑111)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base . (5/5);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.81e-01

TAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTT  >  NC_000913/4017077‑4017370
                                                                                                                                                  |                                                                                                                                                    
tAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCg                                                                                                                                                   <  1:176581/149‑1 (MQ=255)
     ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                              <  3:541866/149‑1 (MQ=255)
     ggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtt                                                                                                                                              <  1:463627/149‑1 (MQ=255)
      gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGACAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCttt                                                                                                                                             >  4:359680/1‑149 (MQ=255)
       tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                                                                            <  3:657002/149‑3 (MQ=255)
       tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCtttt                                                                                                                                            <  1:89384/149‑1 (MQ=255)
          tATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTa                                                                                                                                         <  1:113424/149‑1 (MQ=255)
           aTCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAt                                                                                                                                        >  2:519363/1‑149 (MQ=255)
           aTCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAt                                                                                                                                        >  4:582486/1‑149 (MQ=255)
            tCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTAtt                                                                                                                                       >  4:299366/1‑149 (MQ=255)
             cGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCGATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                                                                             <  3:607364/142‑2 (MQ=255)
             cGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTATCTGAAGGCCGAACGCG‑TCGGCCttc                                                                                                                                             >  4:607364/1‑141 (MQ=255)
               ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                                                                            >  3:99758/1‑139 (MQ=255)
               ttAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCG‑TCGGCCttca                                                                                                                                            <  4:99758/141‑3 (MQ=255)
                         ccAGTAAGAGATCCCGAATGCTGAAACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                          >  1:587790/1‑149 (MQ=255)
                           aGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGTTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                        <  2:569253/149‑1 (MQ=255)
                           aGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                        >  2:684352/1‑149 (MQ=255)
                              aaGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCa                                                                                                                     <  4:361719/149‑1 (MQ=255)
                                gagaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                                                                                   >  1:139912/1‑149 (MQ=255)
                                   aTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                                                                                <  3:299137/149‑1 (MQ=255)
                                      ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCttt                                                                                                                                             >  2:607539/1‑117 (MQ=255)
                                      ccGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCttt                                                                                                                                             <  1:607539/117‑1 (MQ=255)
                                         aaTGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCtt                                                                                                          <  4:562142/149‑1 (MQ=255)
                                            gCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTcc                                                                                                       <  1:177906/149‑1 (MQ=255)
                                                      tGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGt                                                                                             >  3:19898/1‑149 (MQ=255)
                                                      tGAAACAAACCGAAAGCCATGCGGTGAAAAACGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGt                                                                                             >  1:509832/1‑149 (MQ=255)
                                                        aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  2:179708/149‑1 (MQ=255)
                                                        aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  4:24188/149‑1 (MQ=255)
                                                           caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTa                                                                                        <  2:293293/149‑1 (MQ=255)
                                                            aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                                                                 <  1:135399/123‑1 (MQ=255)
                                                            aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaa                                                                                                                 >  2:135399/1‑123 (MQ=255)
                                                            aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                       >  3:348442/1‑149 (MQ=255)
                                                                    aGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAc                                                                                               <  2:411420/133‑1 (MQ=255)
                                                                    aGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAc                                                                                               >  1:411420/1‑133 (MQ=255)
                                                                            ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                                                                       <  4:19898/149‑1 (MQ=255)
                                                                                   aTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg                                                                <  2:139912/149‑1 (MQ=255)
                                                                                   aTCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCggg                                                                >  3:276215/1‑149 (MQ=255)
                                                                                            gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgtt                                                       <  3:296380/149‑1 (MQ=255)
                                                                                                 ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTc                                                  >  2:66628/1‑149 (MQ=255)
                                                                                                  gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                                 >  4:237870/1‑149 (MQ=255)
                                                                                                  gcgcGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCa                                                 <  3:582486/149‑1 (MQ=255)
                                                                                                     cgtcgtCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                              <  1:457601/149‑1 (MQ=255)
                                                                                                      gtcgtcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCaggag                                             <  3:440246/149‑1 (MQ=255)
                                                                                                          tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGt                                         <  2:453016/149‑1 (MQ=255)
                                                                                                            tgctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtg                                       >  2:198644/1‑149 (MQ=255)
                                                                                                             gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgt                                      <  1:198644/149‑1 (MQ=255)
                                                                                                             gctgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgt                                      <  1:684352/149‑1 (MQ=255)
                                                                                                               tgTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                    >  2:330527/1‑149 (MQ=255)
                                                                                                                 tAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGtt                                        >  1:528896/1‑143 (MQ=255)
                                                                                                                 tAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGtt                                        <  2:528896/143‑1 (MQ=255)
                                                                                                                   aTTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAttt                                >  1:545573/1‑149 (MQ=255)
                                                                                                                     tctcttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtct                              >  4:92241/1‑149 (MQ=255)
                                                                                                                         ttctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAAt                          <  1:66628/149‑1 (MQ=255)
                                                                                                                           ctcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCa                        >  4:427721/1‑149 (MQ=255)
                                                                                                                            tcCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                       <  2:587790/149‑1 (MQ=255)
                                                                                                                             ccTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATg                      >  1:318236/1‑149 (MQ=255)
                                                                                                                               tGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGt                    >  4:437597/1‑149 (MQ=255)
                                                                                                                                   tGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTAc                >  2:134703/1‑149 (MQ=255)
                                                                                                                                    gAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACGGCACAATTCTACTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg               >  2:494211/1‑149 (MQ=255)
                                                                                                                                    gAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg               >  3:500098/1‑149 (MQ=255)
                                                                                                                                    aaaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACg               >  3:341565/2‑149 (MQ=255)
                                                                                                                                     aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa                                                             <  4:218226/102‑1 (MQ=255)
                                                                                                                                     aaGGCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa                                                             >  3:218226/1‑102 (MQ=255)
                                                                                                                                       ggCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                              >  1:46873/1‑115 (MQ=255)
                                                                                                                                       ggCCGAACGCG‑TCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                              <  2:46873/115‑1 (MQ=255)
                                                                                                                                             aCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt       <  3:437597/149‑1 (MQ=255)
                                                                                                                                              cgcgTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTg      >  3:279/1‑149 (MQ=255)
                                                                                                                                               gcgTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGc     >  2:693352/1‑149 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt       >  1:487230/1‑144 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt       <  2:487230/144‑1 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGtt  <  4:500098/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCG‑ACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTT  >  NC_000913/4017077‑4017370

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