Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,017,222 Δ1 bp 11.4% 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.11.4% 258.4 / 26.5 79intergenic (+30/‑111)udp/rmuCuridine phosphorylase/DNA recombination protein
Reads supporting (aligned to +/- strand):  ref base T (31/39);  new base . (4/5);  total (35/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 = 9.97e-01

TAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTT  >  NC_000913/4017077‑4017370
                                                                                                                                                  |                                                                                                                                                    
tAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCg                                                                                                                                                   <  1:76572/149‑1 (MQ=255)
   cGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGcc                                                                                                                                                <  2:365617/149‑1 (MQ=255)
    gggTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCt                                                                                                                                               <  4:294487/149‑1 (MQ=255)
      gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttc                                                                                                                                             <  1:491782/149‑2 (MQ=255)
      gtgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCttt                                                                                                                                             <  2:521290/149‑1 (MQ=255)
       tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                            <  2:260434/149‑3 (MQ=255)
       tgtTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttca                                                                                                                                            <  3:60421/149‑3 (MQ=255)
          tATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTAAAGGCCGAACGCGTCGGCCttcagac                                                                                                                                         <  3:326234/149‑6 (MQ=255)
           aTCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAt                                                                                                                                        <  2:199651/149‑1 (MQ=255)
              gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttt                                                                                                                                     <  1:595449/149‑1 (MQ=255)
              gTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTAtttt                                                                                                                                     >  2:89367/1‑149 (MQ=255)
                 aaCCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGc                                                                                                                                  <  2:254959/149‑1 (MQ=255)
                   ccGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGt                                                                                                                                >  2:402187/1‑149 (MQ=255)
                       aCCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAgcg                                                                                                                            >  4:546810/1‑149 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                                                                           <  3:176296/132‑4 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGAACGCGTCGGCCttcag                                                                                                                                           >  4:176296/1‑129 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                          >  1:451979/1‑149 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                          >  4:532753/1‑149 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                          >  4:524834/1‑149 (MQ=255)
                         ccAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGcc                                                                                                                          >  4:530882/1‑149 (MQ=255)
                           aGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTc                                                                                                                        <  4:249090/149‑1 (MQ=255)
                               agagATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAg                                                                                                                    >  2:162679/1‑149 (MQ=255)
                                gagaTCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAgg                                                                                                                   >  1:248897/1‑149 (MQ=255)
                                              tGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCaa                                                                                                     >  2:69179/1‑149 (MQ=255)
                                                 gaCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTg                                                                                                  >  2:575685/1‑149 (MQ=255)
                                                  aCGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTgg                                                                                                 >  1:524627/1‑149 (MQ=255)
                                                    gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                                                                                <  2:436855/132‑1 (MQ=255)
                                                    gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGaaa                                                                                                                >  1:436855/1‑132 (MQ=255)
                                                    gATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGAc                                                                                               >  2:535344/1‑149 (MQ=255)
                                                        aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  1:144291/149‑1 (MQ=255)
                                                        aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  4:202483/149‑1 (MQ=255)
                                                        aaacaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  3:535133/149‑1 (MQ=255)
                                                          acaaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGt                                                                                         >  4:18931/1‑149 (MQ=255)
                                                           caaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTa                                                                                        <  3:60708/149‑1 (MQ=255)
                                                            aaacCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTAc                                                                                       >  2:467068/1‑149 (MQ=255)
                                                               ccGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGc                                                                                    >  3:298752/1‑149 (MQ=255)
                                                                cGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGca                                                                                   <  4:298752/149‑1 (MQ=255)
                                                                        aTGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGACTTTCCAACTGGACGTTTGTACAGCACAATTCTa                                                                           >  2:165629/1‑149 (MQ=255)
                                                                         tGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAt                                                                          <  2:123106/149‑1 (MQ=255)
                                                                         tGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAt                                                                          <  4:336441/149‑1 (MQ=255)
                                                                          gCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtt                                                                         <  2:270623/149‑1 (MQ=255)
                                                                            ggTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAtttt                                                                       >  4:295754/1‑149 (MQ=255)
                                                                             gTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtg                                                                      <  3:18751/149‑1 (MQ=255)
                                                                              tGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTtgt                                                                     <  2:263930/149‑1 (MQ=255)
                                                                                      gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa                                                             <  4:338745/149‑1 (MQ=255)
                                                                                      gtggtgGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTaa                                                             <  2:525832/149‑1 (MQ=255)
                                                                                         gtggAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgtt                                                          <  4:531713/149‑1 (MQ=255)
                                                                                            gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgtt                                                       <  2:524627/149‑1 (MQ=255)
                                                                                            gAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAgttgtt                                                       <  3:387561/149‑1 (MQ=255)
                                                                                              aGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGc                                                     <  1:89367/149‑1 (MQ=255)
                                                                                                 ggcgCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTc                                                  <  3:530882/149‑1 (MQ=255)
                                                                                                   cgcgTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                                <  1:298624/149‑1 (MQ=255)
                                                                                                     cgtcgtCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCagga                                              <  3:18931/149‑1 (MQ=255)
                                                                                                          tcTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGt                                         >  3:489867/1‑149 (MQ=255)
                                                                                                           ctgctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGtt                                        <  1:105195/149‑1 (MQ=255)
                                                                                                              ctgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTtgtg                                     <  4:354890/149‑1 (MQ=255)
                                                                                                               tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGGCGTTTGTACAGCACATTTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTCGTgg                                    >  1:199233/1‑149 (MQ=255)
                                                                                                               tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                                                                 <  3:353416/120‑1 (MQ=255)
                                                                                                               tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCgg                                                                 >  4:353416/1‑120 (MQ=255)
                                                                                                               tgTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTgg                                    >  2:44369/1‑149 (MQ=255)
                                                                                                                gTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGa                                   <  1:208940/149‑1 (MQ=255)
                                                                                                                   aTTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGAttt                                <  1:165629/149‑1 (MQ=255)
                                                                                                                     tctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                                >  4:135900/1‑131 (MQ=255)
                                                                                                                     tctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCag                                                <  3:135900/131‑1 (MQ=255)
                                                                                                                     tctcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTtct                              <  3:546810/149‑1 (MQ=255)
                                                                                                                       tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                            >  1:37479/1‑149 (MQ=255)
                                                                                                                       tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                            >  4:270181/1‑149 (MQ=255)
                                                                                                                       tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                            >  3:59044/1‑149 (MQ=255)
                                                                                                                       tcttctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTATATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCa                            >  1:357524/1‑149 (MQ=255)
                                                                                                                          tctcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATc                         >  3:307768/1‑149 (MQ=255)
                                                                                                                            tcCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCAt                       >  2:9384/1‑149 (MQ=255)
                                                                                                                              cTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATgg                     <  4:307768/149‑1 (MQ=255)
                                                                                                                               tGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTAGTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGt                    <  3:493378/149‑1 (MQ=255)
                                                                                                                                   tGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTAc                <  1:575685/149‑1 (MQ=255)
                                                                                                                                     aaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc              >  2:46908/1‑149 (MQ=255)
                                                                                                                                     aaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc              >  4:261226/1‑149 (MQ=255)
                                                                                                                                     aaGGCCGAACGCGTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGc              >  2:123054/1‑149 (MQ=255)
                                                                                                                                        gccgaaCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAg            >  2:105389/6‑149 (MQ=255)
                                                                                                                                             aCGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTAtt       <  1:61371/149‑1 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           >  1:153429/1‑60 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGttt                                                                                           <  2:153429/60‑1 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGtt  <  1:402187/149‑1 (MQ=255)
                                                                                                                                                  ttCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGtt  <  1:162679/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAGCGGGTGTTATCGTTAACCGCACCCAGCAAGAGATCCCGAATGCTGAGACGATGAAACAAACCGAAAGCCATGCGGTGAAAATCGTGGTGGAAGCGGCGCGTCGTCTGCTGTAATTCTCTTCTCCTGTCTGAAGGCCGACGCGTTCGGCCTTTTGTATTTTTGCGTAGCGCCTCGCAGGAAATGCCTTTCCAACTGGACGTTTGTACAGCACAATTCTATTTTGTGCGGGTAAGTTGTTGCGTCAGGAGGCGTTGTGGATTTCTCAATCATGGTTTACGCAGTTATTGCGTT  >  NC_000913/4017077‑4017370

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