Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,471,357 C→T 15.5% intergenic (‑28/+43) tufA ← / ← fusA translation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,471,3570CT15.5% 364.2 / 47.1 148intergenic (‑28/+43)tufA/fusAtranslation elongation factor EF‑Tu 1/protein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base C (58/66);  new base T (10/13);  total (68/80)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 160 total aligned reads displayed.

ACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATA  >  NC_000913/3471258‑3471456
                                                                                                   |                                                                                                   
aCCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTcc                                                                                                     <  2:2598829/100‑1 (MQ=255)
 ccGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTccc                                                                                                    >  1:1103829/1‑100 (MQ=255)
   gTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTAAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTccc                                                                                                    >  1:2854440/1‑98 (MQ=255)
     ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAACTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                >  1:622016/1‑100 (MQ=255)
     ggTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTc                                                                                                <  1:2425154/100‑1 (MQ=255)
       tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                              <  1:1080817/100‑1 (MQ=255)
       tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                              <  2:16588/100‑1 (MQ=255)
       tCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAg                                                                                              >  2:2533376/1‑100 (MQ=255)
        cAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAgg                                                                                             >  2:1603950/1‑100 (MQ=255)
         aaCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTCCAgg                                                                                             >  1:2280841/1‑99 (MQ=255)
          aCGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgag                                                                                           >  1:277909/1‑100 (MQ=255)
            gTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagag                                                                                         <  1:2647400/100‑1 (MQ=255)
             tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                        <  1:1121775/100‑1 (MQ=255)
             tGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGgagaga                                                                                        <  1:2526151/100‑1 (MQ=255)
               gCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGc                                                                                      <  2:2359851/100‑1 (MQ=255)
                ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCa                                                                                     <  2:69945/100‑1 (MQ=21)
                ccGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCa                                                                                     >  1:3380127/1‑100 (MQ=255)
                ccGATAGTACCAACGTCAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCa                                                                                     <  1:2231337/100‑1 (MQ=21)
                 cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                    <  1:744818/100‑1 (MQ=25)
                 cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                    <  1:280967/100‑1 (MQ=25)
                 cGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAAAGAGCa                                                                                     >  1:3363836/1‑99 (MQ=255)
                  gATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAct                                                                                   <  2:1554790/100‑2 (MQ=25)
                   aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTGAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCAc                                                                                    <  2:454170/98‑1 (MQ=25)
                   aTAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACg                                                                                   <  2:638244/99‑1 (MQ=255)
                    tAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCTCCTGAAGGGGAGAGCActat                                                                                 <  2:262395/100‑4 (MQ=25)
                     aGTACCAACGTTAACGTGCGGCTTTGTACGTTCAAATTTTTCTTTAGACACGCCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACggg                                                                                 <  1:1635460/99‑1 (MQ=255)
                      gTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGAc                                                                               <  1:2072836/100‑1 (MQ=255)
                       tACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACt                                                                              >  2:1108133/1‑100 (MQ=255)
                         ccAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACttt                                                                            >  1:2698896/1‑100 (MQ=255)
                          cAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTCCTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTg                                                                           >  1:409630/1‑100 (MQ=255)
                           aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                          <  1:2284882/100‑1 (MQ=255)
                           aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                          <  2:2212109/100‑1 (MQ=255)
                           aaCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                          <  2:376910/100‑1 (MQ=255)
                            aCGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTGGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                          <  1:2984400/99‑1 (MQ=255)
                             cGTTAACGTGCTGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                        <  2:102218/100‑1 (MQ=255)
                             cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTgg                                                                          <  1:1596802/98‑1 (MQ=255)
                             cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTa                                                                        <  2:474545/100‑1 (MQ=255)
                             cGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCCCCCCTTCAGGAGAGAGCTCGGGACTTTggca                                                                        <  2:2568441/100‑3 (MQ=255)
                               ttAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTAtt                                                                      >  2:3193833/1‑100 (MQ=255)
                                tAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTa                                                                     <  1:358739/100‑1 (MQ=255)
                                  aCGTGCTGTTTTGTACGTTCAAATTTTTCTTTCGACACGGCTATATTCAGTACTATAGTGCTCTCCACTTCAGGAGAGTGCACGGGACTTTGGTATTAAc                                                                   >  1:2595987/1‑100 (MQ=38)
                                  aCGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAAc                                                                   <  2:1332040/100‑1 (MQ=255)
                                     tGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCt                                                                >  1:1979282/1‑100 (MQ=255)
                                         gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                            >  1:2009931/1‑100 (MQ=255)
                                         gTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAgg                                                            <  1:1499857/100‑1 (MQ=255)
                                           tttGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCt                                                          >  1:953956/1‑100 (MQ=255)
                                             tGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTa                                                        <  1:2158085/100‑1 (MQ=255)
                                             tGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTACCCTTCAGAAGAGGGCACGGGACTTTGGTATTAACCCTTCGGCTTa                                                        >  1:187862/1‑100 (MQ=255)
                                                aCGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTAttt                                                     >  1:1757445/1‑100 (MQ=255)
                                                   ttCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTTTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAcc                                                  >  1:510577/1‑100 (MQ=255)
                                                   ttCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTAcc                                                  <  2:1567341/100‑1 (MQ=255)
                                                      aaaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACg                                               >  2:2969310/1‑100 (MQ=255)
                                                       aaTTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACgg                                              <  2:3405575/100‑1 (MQ=255)
                                                           tttCTTGAGACACGGCTATATTCCTTACTATAGTGCTCTCCACTTCAGGAGAGAGCACGGGACTTTGATATTAACCCTTAGGCTTATTTACCACGGGCt                                           <  1:1454861/99‑1 (MQ=255)
                                                             tCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCa                                        >  2:2489706/1‑100 (MQ=255)
                                                               tttAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAAt                                      >  1:1245555/1‑100 (MQ=255)
                                                                    acacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCCTAGGCTTATTTACCACGGGCTTCAATTACgg                                 <  2:2403296/100‑1 (MQ=255)
                                                                     cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                >  2:33108/1‑100 (MQ=255)
                                                                     cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                >  1:3433287/1‑100 (MQ=255)
                                                                     cacGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGc                                >  2:3549973/1‑100 (MQ=255)
                                                                      acGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGcc                               <  2:2803305/100‑1 (MQ=255)
                                                                        ggCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTg                             <  1:3299440/100‑1 (MQ=255)
                                                                        ggCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTg                             >  1:1269626/1‑100 (MQ=255)
                                                                              atTCCTTACTTTAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAAc                       <  1:1820319/100‑1 (MQ=255)
                                                                               ttCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACg                      <  2:1466648/100‑1 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACTGGCTTCAATTACGGCCTGAGCAACGt                     <  2:1305944/100‑1 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                     <  1:476226/100‑1 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                     <  1:3343469/100‑1 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                     <  2:3311950/100‑1 (MQ=255)
                                                                                tCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGt                     <  1:2166698/100‑1 (MQ=255)
                                                                                 ccTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGtt                    >  2:927713/1‑100 (MQ=255)
                                                                                 ccTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGtt                    <  2:3427971/100‑1 (MQ=255)
                                                                                  cTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTa                   >  2:3048566/1‑100 (MQ=255)
                                                                                    tACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACt                 <  1:176780/100‑1 (MQ=255)
                                                                                     aCTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTc                >  2:2327045/1‑100 (MQ=255)
                                                                                     aCTATAGTGCCCATACCTTCAAGAGAGAGTACGGGACTTTTCCATTATCCCTTAGGCTTATTTACCACCGTCTTCAATTACGGCCTGAGCAACGTTAc                  <  1:1347220/98‑1 (MQ=255)
                                                                                        cccgTGCTCTCTCCTGAAGGGGAGAGCACAGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  2:3759679/4‑100 (MQ=37)
                                                                                        ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCgg              >  2:22810/1‑99 (MQ=255)
                                                                                        ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  2:2393423/1‑100 (MQ=255)
                                                                                        ataGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  2:2230453/1‑100 (MQ=255)
                                                                                        agaGTGCTCTTCCCTTCAGGAGGGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGATTCAATTACGGCCTGAGCAACGTTACTc                <  1:112355/95‑1 (MQ=255)
                                                                                         ccgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcg            >  1:961037/3‑100 (MQ=38)
                                                                                          cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgcgc           >  2:1210210/2‑100 (MQ=38)
                                                                                          cgTGCTCTCTCCTGAAGGGGAGAGCACGGGACCTTGGTATTAGCCCTTAGGCTTATTTACCACGTGCTTCAATTACGGCGTGAGCAACGTTACTCGgc             >  1:3422479/2‑98 (MQ=25)
                                                                                           gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt          <  1:3534040/100‑1 (MQ=255)
                                                                                           gTGCTCTCCCCTTCAGGAGAGAGCACGGGACTGTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCt          <  1:1580963/100‑1 (MQ=255)
                                                                                            tGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGgc             >  1:1188074/1‑96 (MQ=39)
                                                                                            tGCTCTCTCCTGAAGGGGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTAGGCGCtt         >  2:873824/1‑100 (MQ=255)
                                                                                             gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc        >  2:2643584/1‑100 (MQ=255)
                                                                                             gCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtc        >  1:152964/1‑100 (MQ=255)
                                                                                              ctctCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca       <  1:3645724/100‑1 (MQ=255)
                                                                                              ctctCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca       >  1:3766551/1‑100 (MQ=255)
                                                                                              ctatcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtca       <  2:920369/97‑1 (MQ=255)
                                                                                               tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat      >  2:2731508/1‑100 (MQ=255)
                                                                                               tctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcat      <  1:1558646/100‑1 (MQ=255)
                                                                                                ctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc     <  1:3720410/100‑1 (MQ=255)
                                                                                                ctcCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTtcatc     >  2:283592/1‑100 (MQ=255)
                                                                                                  ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat   >  2:573741/1‑100 (MQ=255)
                                                                                                  ccccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTcatcat   <  2:2698909/100‑1 (MQ=255)
                                                                                                   cccTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATa  >  2:1127983/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ACCGTGGTCAACGTGGCCGATAGTACCAACGTTAACGTGCGGTTTTGTACGTTCAAATTTTTCTTTAGACACGGCTATATTCCTTACTATAGTGCTCTCCCCTTCAGGAGAGAGCACGGGACTTTGGTATTAACCCTTAGGCTTATTTACCACGGGCTTCAATTACGGCCTGAGCAACGTTACTCGGCGCTTCATCATA  >  NC_000913/3471258‑3471456

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