Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,554,516 A→G 6.4% intergenic (+155/‑60) uxuB → / → uxuR D‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,554,5160AG 6.4% 185.2 / 3.3 64intergenic (+155/‑60)uxuB/uxuRD‑mannonate oxidoreductase, NAD‑dependent/fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base G (2/2);  total (32/32)
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.98e-01

TTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCT  >  NC_000913/4554370‑4554662
                                                                                                                                                  |                                                                                                                                                  
ttACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAg                                                                                                                                                  <  1:766006/149‑1 (MQ=255)
 taaccTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGCCCTGGTTTCGGTTTTTTGCCCACCTGTACGTGCCAACTACCAGTGCTAATGGTAAAGAATGAGATTAACGGGGGCCGTAAAATTGCACGTTGTAGGCCGGATAAgg                                                                                                                                                 >  1:339129/4‑149 (MQ=255)
  aCCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTTTAGGCCGGATAATGc                                                                                                                                                >  1:874067/1‑149 (MQ=255)
   cccTACAGACTTACTGGTCAATCAAACTGATATTTGGTTCACCAGTTTTCGTTATTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAcgcg                                                                                                                                               >  1:685865/1‑149 (MQ=255)
      tACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCCCCTCTACGTGCCAACGTCCAGTGCTAATGGTATAGTTTGAGATTACCGGTGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTc                                                                                                                                            >  2:736309/1‑149 (MQ=255)
      tACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGtga                                                                                                                                            >  2:521324/1‑147 (MQ=255)
      tACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTc                                                                                                                                            <  1:349366/149‑1 (MQ=255)
       aCAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                           >  2:302323/1‑149 (MQ=255)
       aCAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCa                                                                                                                                           >  2:392609/1‑149 (MQ=255)
        cAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAc                                                                                                                                          >  2:325483/1‑149 (MQ=255)
          caCTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgc                                                                                                                                        <  2:430032/148‑1 (MQ=255)
             ttACTGGTCAATCAAACTGAGATTTGGTTGACCAGTTTTCGTTTTTTTTCCCACCTGACCGTGCCAACTTCCAGTGCAAATGGTATAGTTTGAGATTAACGGGGGCCTTAAAATTGCCCGTTGTAGTCCGTATAAGGCGTTCAcgccgc                                                                                                                                     >  1:983429/1‑149 (MQ=255)
             ttACTGGTAACTCAACCTGATATTTGGTTGACCAGTTTTCGATTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTCTAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCAcgccgc                                                                                                                                     <  1:581949/149‑1 (MQ=255)
               aCTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAAcgccttat                                                                                                                                   <  1:521324/149‑5 (MQ=255)
               aCTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACATCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCAt                                                                                                                                   <  2:367135/149‑1 (MQ=255)
                cTGGGCACTCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATAAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATc                                                                                                                                  <  1:879694/149‑1 (MQ=255)
                   gTCAATCAAACTGATATTTTGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGTGGCCGTAAAATTGACCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                               >  2:815330/1‑149 (MQ=255)
                   gTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                               >  2:376310/1‑149 (MQ=255)
                    ccactcaaACTGAAATTTGCTTGACCAGTTTTCGTTTTATAGCCCACCTGTACGTGCCAACTTCCAGGGCTAATGGTATAGTTTGAGATTAACGGCGGCCGTAAAATTGCCCGTTGTAGGCCGGATGCGGCGTGAACGCCTTATCCGGc                                                                                                                              <  1:671479/145‑1 (MQ=255)
                      aatcaaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCaa                                                                                                                            >  2:704903/1‑149 (MQ=255)
                       atcaaACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGATAATGTTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCTGATAAGGCGTTCACGCCGCATCCGGCaaa                                                                                                                           >  2:638755/1‑149 (MQ=255)
                          aaaCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAt                                                                                                                        <  2:49392/149‑1 (MQ=255)
                            aCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATTGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                      <  2:884500/149‑1 (MQ=255)
                            aCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                      <  2:355798/149‑1 (MQ=255)
                            aCTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAAttt                                                                                                                      <  1:449103/149‑1 (MQ=255)
                                 tatTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTa                                                                                                                 >  1:575391/1‑149 (MQ=255)
                                     tGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCg                                                                                                             <  1:686924/149‑1 (MQ=255)
                                                     gTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaa                                                                                             <  1:916574/149‑1 (MQ=255)
                                                      tttttttGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACacaac                                                                                            >  2:365595/1‑149 (MQ=255)
                                                            tGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATg                                                                                      <  1:956116/149‑1 (MQ=255)
                                                                cACCTGTACGTGCCAACTTCCAGTGCTAATGGTTTAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCAGCCGGCAAAAAGTTGATTAGCCGCACCTTAGGGTCGCAGCACCATGAAAt                                                                                  >  2:1017759/1‑149 (MQ=38)
                                                                cACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGatat                                                                                  <  2:983429/149‑1 (MQ=255)
                                                                        cGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCAcc                                                                          <  1:983398/149‑1 (MQ=255)
                                                                         gTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCccccg                                                                         >  1:1024212/1‑148 (MQ=255)
                                                                              aaCTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCACACCTAACGGACACAACACCATGAAATCTGCCACCTCTgcg                                                                    <  1:913625/149‑1 (MQ=255)
                                                                               aCTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAAGTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAATTTGAGTAACCGCCCCTACCGGACACAACACCATGAAATCTGCCACCTCTgcgc                                                                   >  2:195983/1‑149 (MQ=39)
                                                                                   ccAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAg                                                               <  1:63045/149‑1 (MQ=255)
                                                                                   caaGTGCTAATGGTTTAGTTTTAGATTAACGGGGGCCGTAAAATTGACCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATACGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAg                                                               <  1:914783/147‑1 (MQ=255)
                                                                                     aGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGAc                                                             <  1:1012129/149‑1 (MQ=255)
                                                                                       tGCTAATGGTATATTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCt                                                           <  1:704903/149‑1 (MQ=255)
                                                                                           aaTGGTATAGTTGGAGATTAACTGGGGCCGTAAAATTGCCCGTTGTAGGCCGGAGATGGCGTTCACGCCGCATCAGGCAAAATTTTGATTATCCGCACCTAACGGACCCAACACAATGAATTCTGCCACCGCTGCTCATAGACCTTAcc                                                       >  1:107636/1‑149 (MQ=255)
                                                                                              ggTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAgg                                                    <  1:263111/149‑1 (MQ=255)
                                                                                                tataGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGaa                                                  >  1:6764/1‑149 (MQ=255)
                                                                                                tataGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGaa                                                  <  2:575391/149‑1 (MQ=255)
                                                                                                tataGTGTGACATTACCGGGGGCCGTAAAATTGCCCGTTGTTGGCCGGCTTCGGCGTTCACGCCGCATCCGGCAAAAATTTGATTACCCGCACCTAACGGACCCAACACCATGAAATCTGCCACCTCTGCGCACAGACCTTACCAGGaa                                                  <  2:838989/149‑1 (MQ=255)
                                                                                                    gTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCATGATAAGGCGATCACGCCGCATCCGGCAAAAATTTGTTTACACGCACCTAAGGGACACAACACCATGAAATCTTCCTCCTCTGCGCAAAGACCTTACCAGGAAGTCg                                              <  2:361554/149‑1 (MQ=255)
                                                                                                     tttGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCAGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGCCACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCgg                                             >  2:566342/1‑149 (MQ=255)
                                                                                                       tGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACCACACCAAGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCgggg                                           >  1:908834/1‑149 (MQ=255)
                                                                                                        gagaTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTACCCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTGACCAGGAAGTCGGGGc                                          >  1:721533/1‑149 (MQ=255)
                                                                                                           aTTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgat                                       <  2:272070/149‑1 (MQ=255)
                                                                                                             tAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTCACGGACACAACACAATGAAATCTGCCACCTCTGCGCAAAGCCCTTACCAGGAAGTCGGGGCgatga                                     >  1:449033/1‑149 (MQ=255)
                                                                                                              aaCGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCgatgat                                    >  1:314888/1‑149 (MQ=255)
                                                                                                               aCGGGGGCCGTAAAATTGCCCTGTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGAACCTAACCGAGACAACACCATGAAATCTGCCACCTCTGCGCAAAGAACTTACCAGGAAGTCGGGGCgatgatg                                   <  2:685865/149‑2 (MQ=38)
                                                                                                                       tgTAAAATTGCCCGTTGTAGTCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCTCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCt                           <  1:687526/148‑1 (MQ=255)
                                                                                                                             aTTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCAt                     <  2:314888/149‑1 (MQ=255)
                                                                                                                                gCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACCACACCATGAAATCTGCCCCCTCTGCGCAAAGCCCTTACCAGGATGTCGGGGCGATGATCCGCGAACGGAGCATaaa                  >  2:43855/1‑149 (MQ=38)
                                                                                                                                  ccGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGATGTCGGGGCGATGATCCGCGATCTGATCATAAAGa                >  1:417958/1‑149 (MQ=255)
                                                                                                                                      tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACGTCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGTCGATGATCCGCTATCTGATCATAAAGACGcc            <  1:815330/149‑1 (MQ=255)
                                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTATCGGACAAAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTa         >  1:145489/1‑149 (MQ=255)
                                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTa         >  1:381564/1‑149 (MQ=255)
                                                                                                                                          ggCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTAc        >  2:551108/1‑149 (MQ=255)
                                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACa       >  1:349213/1‑149 (MQ=255)
                                                                                                                                            ccGGATGCGGCGTGAACGCGTTATCCGGCAAAAATTTGATTAACCGCACCTCAAGGTCACAACCCCATGAAATCTGCCACCTCTGCGCAAAGACCGTACCAGGAAGTAGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTAcca      >  2:602544/1‑147 (MQ=255)
                                                                                                                                                aTAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCt  <  1:736309/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TTACCCTACAGACTTACTGGTCAATCAAACTGATATTTGGTTGACCAGTTTTCGTTTTTTTGCCCACCTGTACGTGCCAACTTCCAGTGCTAATGGTATAGTTTGAGATTAACGGGGGCCGTAAAATTGCCCGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAAAAATTTGATTAACCGCACCTAACGGACACAACACCATGAAATCTGCCACCTCTGCGCAAAGACCTTACCAGGAAGTCGGGGCGATGATCCGCGATCTGATCATAAAGACGCCGTACAATCCT  >  NC_000913/4554370‑4554662

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑