Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,073,307 Δ1 bp 100% coding (705/1149 nt) rutA ← pyrimidine oxygenase, FMN‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,073,3050A.96.2% 449.9 / 8.0 104coding (707/1149 nt)rutApyrimidine oxygenase, FMN‑dependent
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (49/51);  total (52/52)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.92e-01

Only 100 of 109 total aligned reads displayed.

GTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAG  >  NC_000913/1073208‑1073400
                                                                                                 |                                                                                               
gtcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaag                                                                                               >  2:1053169/1‑99 (MQ=255)
 tcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaagc                                                                                              <  1:2875578/100‑3 (MQ=255)
 tcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaagc                                                                                              <  1:208115/100‑3 (MQ=255)
 tcgGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGGTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGaagc                                                                                              <  1:752742/100‑3 (MQ=255)
    gTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGt                                                                                          >  1:2794965/1‑100 (MQ=255)
    gTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGt                                                                                          >  1:3048679/1‑100 (MQ=255)
    gTTTCATCGGCAATCACCATAAACAACACACAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGt                                                                                          <  1:81628/100‑1 (MQ=255)
     tttCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTg                                                                                         <  2:125401/100‑1 (MQ=255)
       tCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTggg                                                                                       <  1:2990146/100‑1 (MQ=255)
       tCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTggg                                                                                       <  1:1579829/100‑1 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgt                                                                                      <  2:2162450/100‑1 (MQ=255)
        cATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgt                                                                                      <  2:122860/100‑1 (MQ=255)
         aTCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtg                                                                                     >  1:2166040/1‑100 (MQ=255)
         aTCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtg                                                                                     >  1:1972991/1‑100 (MQ=255)
         aTCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtg                                                                                     >  2:1504509/1‑100 (MQ=255)
         aTCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCGGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtg                                                                                     <  1:2486417/100‑1 (MQ=255)
          tCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtgt                                                                                    >  2:2264080/1‑100 (MQ=255)
           cGGCAATCACCATATACAACACATAAGATCCATCGTCGCGTCCGGTTTGCGCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtgtn                                                                                   <  1:2407735/100‑2 (MQ=255)
           cGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGgtgt                                                                                    <  2:1330394/99‑1 (MQ=255)
           cGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTa                                                                                   <  2:2813540/100‑1 (MQ=255)
            ggCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAt                                                                                  >  2:2622165/1‑100 (MQ=255)
             gCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAtt                                                                                 >  2:2963413/1‑100 (MQ=255)
             gCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAtt                                                                                 <  1:259001/100‑1 (MQ=255)
              cAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAttt                                                                                <  1:1745850/100‑1 (MQ=255)
              cAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAttt                                                                                >  1:1158441/1‑100 (MQ=255)
               aaTCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTa                                                                               <  2:625814/100‑1 (MQ=255)
                 tCACCATAAACAACACATATGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACg                                                                             <  1:2406795/100‑1 (MQ=255)
                 tCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACg                                                                             >  2:2085285/1‑100 (MQ=255)
                  cACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGc                                                                            <  2:675867/100‑1 (MQ=255)
                   actttAAACAACACATAAGATCCAACGTCGCTTCCGGTTTGCTCTGTTGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGcc                                                                           <  1:119929/96‑1 (MQ=255)
                   aCCATAAACAACACATAAGTTCCAACGTCGCGTTCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGcc                                                                           <  1:585288/100‑1 (MQ=255)
                   aCCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGcc                                                                           >  2:2280675/1‑100 (MQ=255)
                   aCCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGcc                                                                           >  1:1864267/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCtt                                                                         >  2:2295571/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCtt                                                                         >  1:2824257/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCtt                                                                         >  2:287249/1‑100 (MQ=255)
                     cATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCtt                                                                         >  1:367759/1‑100 (MQ=255)
                       tAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCttt                                                                        >  1:2227322/1‑100 (MQ=255)
                         aacaacACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGcc                                                                     <  2:2940988/100‑1 (MQ=255)
                          acaacaCATAAGAGTTAACGTCGCGTTCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCCGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCg                                                                    >  2:2628274/1‑100 (MQ=255)
                          acaacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCg                                                                    >  1:3018626/1‑100 (MQ=255)
                            aacaCATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGaa                                                                  <  1:1923032/100‑1 (MQ=255)
                             acacATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGaaa                                                                 >  1:2591769/1‑100 (MQ=255)
                              cacaTAAGTGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAAc                                                                <  1:1774695/100‑1 (MQ=255)
                              cacaTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAAc                                                                <  2:1742729/100‑1 (MQ=255)
                                caTAAGAGCCAACGAAGCGACCGGAATGCTAAGCAGCCTGTTTCATACGCGCAGCGGTCGGGGC‑CAAGACGTGGGCGTATTTACGCCTTTGCCGAAACAg                                                              >  2:286067/1‑100 (MQ=255)
                                 aTAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGa                                                             >  1:1551703/1‑100 (MQ=255)
                                  tAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGaaa                                                                 >  2:2082441/1‑95 (MQ=255)
                                  tAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGaa                                                            <  1:576359/100‑1 (MQ=255)
                                   aaGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAg                                                           <  1:2499030/100‑1 (MQ=255)
                                     gagCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGtt                                                         <  1:1292015/100‑1 (MQ=255)
                                      agCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTg                                                        >  2:1412443/1‑100 (MQ=255)
                                      agCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTg                                                        <  1:1354770/100‑1 (MQ=255)
                                        ccAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGaa                                                      >  1:2808577/1‑100 (MQ=255)
                                          aaCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAAt                                                    >  1:313906/1‑100 (MQ=255)
                                          aaCGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAAt                                                    >  1:318068/1‑100 (MQ=255)
                                             gTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCgg                                                 <  2:2972039/100‑1 (MQ=255)
                                             gTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCgg                                                 <  2:2825532/100‑1 (MQ=255)
                                             gTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCgg                                                 <  1:1412457/100‑1 (MQ=255)
                                               cgcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGTGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGc                                                <  1:1629753/100‑1 (MQ=255)
                                               cgcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCa                                               <  2:2875141/100‑1 (MQ=255)
                                                gcgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTAGTTACGCCTTTGCCGAAACAGAAGTTGAAATCCGCAt                                              >  2:2580578/1‑100 (MQ=255)
                                                 cgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATa                                             >  2:325506/1‑100 (MQ=255)
                                                 cgTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCATTTGCCGAAACAGAAGTTGAAATCGGCATa                                             >  1:2458254/1‑100 (MQ=255)
                                                 cgTCCGGTTTGCTCTGCGGCCTGTTCCATCCTCGCAGCGTTCGGGACG‑AAGCCGTGGGTGTCTTTACGCCTTTGCCTAAACAGAAGTTGAAATCGGc                                                <  2:1321953/97‑1 (MQ=255)
                                                  gTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATAc                                            <  1:287252/100‑1 (MQ=255)
                                                    ccGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCg                                          >  1:838056/1‑100 (MQ=255)
                                                     cGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCgg                                         >  2:695484/1‑100 (MQ=255)
                                                       gTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGc                                       >  1:811677/1‑100 (MQ=255)
                                                         ttGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCgg                                     <  2:579152/100‑1 (MQ=255)
                                                         ttGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCgg                                     <  1:2684324/100‑1 (MQ=255)
                                                          tGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCgg                                     >  2:1336259/1‑100 (MQ=255)
                                                          tGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGga                                    >  2:2412141/1‑100 (MQ=255)
                                                           gCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgag                                   >  2:2517975/1‑100 (MQ=255)
                                                            ctctGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGgaga                                  >  1:1519626/1‑100 (MQ=255)
                                                             tctGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGaa                                 >  1:2471370/1‑100 (MQ=255)
                                                              ctGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGAGGAGAAc                                <  1:1589988/100‑1 (MQ=255)
                                                                    ccTGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATg                          <  1:695488/100‑1 (MQ=255)
                                                                      tGTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGcc                        <  1:2306319/100‑1 (MQ=255)
                                                                       gTTTCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGcc                        >  1:506485/1‑99 (MQ=255)
                                                                          tCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                    <  1:1328525/100‑1 (MQ=255)
                                                                          tCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                    <  1:466959/100‑1 (MQ=255)
                                                                          tCATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCg                    <  2:2794933/100‑1 (MQ=255)
                                                                           cATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGt                   >  2:2262275/1‑100 (MQ=255)
                                                                           cATCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGt                   >  2:1187763/1‑100 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTc                  >  2:1319749/1‑100 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTc                  <  1:2264105/100‑1 (MQ=255)
                                                                            aTCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGACATACCGGGCGGAGAACGCCATGCCAGCGTc                  >  1:1289052/1‑100 (MQ=255)
                                                                             tCCGCGCAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCg                 <  2:782314/100‑1 (MQ=255)
                                                                                  gcAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTg            <  1:2232184/100‑1 (MQ=255)
                                                                                   cAGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGc           <  2:2157272/100‑1 (MQ=255)
                                                                                    aGCGGTCGGGGCG‑CAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGcc          >  2:3018044/1‑100 (MQ=255)
                                                                                    aGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCTATGCCAGCGTCGCTTTGcg          <  2:1341845/100‑2 (MQ=255)
                                                                                    aGCGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCAGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGcc          >  1:2701939/1‑100 (MQ=255)
                                                                                      cGGTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCg        <  2:1474067/100‑1 (MQ=255)
                                                                                        gTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCgg       <  2:2591737/99‑1 (MQ=255)
                                                                                        gTCGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACACCATGCCAGCGTCGCTTTGCCCGgc      >  2:509793/1‑100 (MQ=255)
                                                                                          cGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgcgc    >  1:654699/1‑100 (MQ=255)
                                                                                          cGGGGCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGgcgc    <  2:2399055/100‑1 (MQ=255)
                                                                                            gggCG‑AAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAg  >  1:2417131/1‑100 (MQ=255)
                                                                                                 |                                                                                               
GTCGGTTTCATCGGCAATCACCATAAACAACACATAAGAGCCAACGTCGCGTCCGGTTTGCTCTGCGGCCTGTTTCATCCGCGCAGCGGTCGGGGCGAAAGCCGTGGGTGTATTTACGCCTTTGCCGAAACAGAAGTTGAAATCGGCATACCGGGCGGAGAACGCCATGCCAGCGTCGCTTTGCCCGGCGCAG  >  NC_000913/1073208‑1073400

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑