Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,325,037:1 +T 7.3% coding (282/1323 nt) atoE → short chain fatty acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,325,0371.T 7.3% 379.5 / 27.0 96coding (282/1323 nt)atoEshort chain fatty acid transporter
Reads supporting (aligned to +/- strand):  ref base . (46/43);  new base T (3/4);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.59e-01

CATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCT  >  NC_000913/2324890‑2325185
                                                                                                                                                    |                                                                                                                                                    
cATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATTCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑g                                                                                                                                                     >  2:912019/1‑149 (MQ=255)
    gTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACa                                                                                                                                                 <  1:167797/149‑1 (MQ=255)
      gAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAgg                                                                                                                                               <  3:177892/149‑1 (MQ=255)
            gtgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTct                                                                                                                                         <  4:239998/149‑2 (MQ=255)
            gtgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCa                                                                                                                                         >  1:91846/1‑149 (MQ=255)
             tgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTctt                                                                                                                                        <  1:271449/149‑3 (MQ=255)
             tgGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCAt                                                                                                                                        >  3:568004/1‑149 (MQ=255)
              ggggTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTAC‑GGGCGTcttt                                                                                                                                       <  4:116228/149‑4 (MQ=255)
                 gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTg                                                                                                                                    >  2:838851/1‑149 (MQ=255)
                 gTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTg                                                                                                                                    >  2:956848/1‑149 (MQ=255)
                  tGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTgg                                                                                                                                   >  3:570665/1‑149 (MQ=255)
                      ggTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCAc                                                                                                                               <  2:5022/149‑1 (MQ=255)
                       gTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACt                                                                                                                              >  2:649109/1‑149 (MQ=255)
                        tttCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACtt                                                                                                                             >  2:453511/1‑149 (MQ=255)
                         ttCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACttt                                                                                                                            <  2:363757/149‑1 (MQ=255)
                             ggAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                        <  4:639091/149‑1 (MQ=255)
                             ggAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                        <  3:516372/149‑1 (MQ=255)
                              gAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCg                                                                                                                       <  1:1021081/149‑1 (MQ=255)
                              gAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCg                                                                                                                       <  3:934271/149‑1 (MQ=255)
                                    gCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAg                                                                                                                 <  3:250747/149‑1 (MQ=255)
                                      tGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTc                                                                                                               >  1:628297/1‑149 (MQ=255)
                                          gTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCtt                                                                                                           <  3:104543/149‑1 (MQ=255)
                                            ttGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTtgt                                                                                                         <  3:1040683/149‑1 (MQ=255)
                                              ggTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTCGTCACTTTCTTCGGTTCAGTCGCTtgtct                                                                                                       <  3:898084/149‑3 (MQ=255)
                                                tATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCACAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtca                                                                                                     <  4:758942/149‑1 (MQ=255)
                                                  tGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                        <  4:191817/128‑1 (MQ=255)
                                                  tGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTttcttc                                                                                                                        >  3:191817/1‑128 (MQ=255)
                                                  tGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGtcatc                                                                                                   <  2:716663/149‑1 (MQ=255)
                                                    cAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCaa                                                                                                 >  2:334477/1‑149 (MQ=255)
                                                      gATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACt                                                                                               <  4:803192/149‑1 (MQ=255)
                                                            gcTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAt                                                                                         >  1:336931/1‑149 (MQ=255)
                                                               tATCATCGTTACCGGCCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttg                                                                                      <  4:387365/149‑1 (MQ=255)
                                                                 tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggt                                                                                    >  3:533067/1‑149 (MQ=255)
                                                                 tcatcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggt                                                                                    >  3:934882/1‑149 (MQ=255)
                                                                   atcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                  >  4:59278/1‑149 (MQ=255)
                                                                   atcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggttt                                                                                  <  4:568004/149‑1 (MQ=255)
                                                                    tcGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtttg                                                                                 >  3:359059/1‑149 (MQ=255)
                                                                     cGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGAtttggtttgg                                                                                <  1:128058/149‑1 (MQ=255)
                                                                       ttACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCtt                                                                                                           >  2:583503/1‑120 (MQ=255)
                                                                       ttACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCtt                                                                                                           <  1:583503/120‑1 (MQ=255)
                                                                       ttACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTtggtttggtt                                                                              <  1:1010156/149‑1 (MQ=255)
                                                                         aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTA‑AGGGCGTcttt                                                                                                                                       >  2:66973/1‑87 (MQ=255)
                                                                         aCCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCCTGTA‑AGGGCGTcttt                                                                                                                                       <  1:66973/90‑4 (MQ=255)
                                                                             gTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgc                                                                        <  4:245606/149‑1 (MQ=255)
                                                                             gTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgc                                                                        <  3:13735/149‑1 (MQ=255)
                                                                             gTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGgc                                                                        >  2:529694/1‑149 (MQ=255)
                                                                                aTGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCa                                                                     >  3:838510/1‑149 (MQ=255)
                                                                                  gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTg                                                                                                                                    <  4:852346/84‑1 (MQ=255)
                                                                                  gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTg                                                                                                                                    >  3:852346/1‑84 (MQ=255)
                                                                                  gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACtt                                                                                                                             >  4:752429/1‑91 (MQ=255)
                                                                                  gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACtt                                                                                                                             <  3:752429/91‑1 (MQ=255)
                                                                                  gCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAAt                                                                   >  3:325456/1‑149 (MQ=255)
                                                                                   cccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCACTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATg                                                                  <  4:853625/149‑1 (MQ=255)
                                                                                    ccTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGt                                                                 >  4:262592/1‑149 (MQ=255)
                                                                                     cTTGCCAGCTCTGCCCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGtt                                                                >  1:798502/1‑149 (MQ=255)
                                                                                        gCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTTCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGGTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTTATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGc                                                             <  2:605101/149‑1 (MQ=255)
                                                                                        gCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCGCTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGc                                                             <  4:74502/149‑1 (MQ=255)
                                                                                            gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGt                                                         >  4:245365/1‑149 (MQ=255)
                                                                                            gCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGt                                                         >  2:593876/1‑149 (MQ=255)
                                                                                              tctGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGa                                                       >  1:317471/1‑149 (MQ=255)
                                                                                              tctGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGa                                                       >  1:732221/1‑149 (MQ=255)
                                                                                              tctGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGa                                                       <  4:838510/149‑1 (MQ=255)
                                                                                                tGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAg                                                     <  3:105667/149‑1 (MQ=255)
                                                                                                  cTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTc                                                   <  3:662460/149‑1 (MQ=255)
                                                                                                  cTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTc                                                   >  1:158030/1‑149 (MQ=255)
                                                                                                      ggTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGccc                                               <  4:36539/149‑1 (MQ=255)
                                                                                                       gTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCg                                              >  4:429922/1‑149 (MQ=255)
                                                                                                        tGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCgg                                             >  2:972826/1‑149 (MQ=255)
                                                                                                         gAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTgg                                                                                             >  4:406489/1‑100 (MQ=255)
                                                                                                         gAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTgg                                                                                             <  3:406489/100‑1 (MQ=255)
                                                                                                           aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGa                                          <  2:474979/149‑1 (MQ=255)
                                                                                                           aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGa                                          >  4:562892/1‑149 (MQ=255)
                                                                                                           aaaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGa                                          <  2:842024/149‑1 (MQ=255)
                                                                                                            aaGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAg                                         >  4:642835/1‑149 (MQ=255)
                                                                                                              gTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCAATTTCTTCGGTTCAGTTGCTTGTGTCAGCAACCGGGGGTTTGGTTTGGTTGTCGGCGCAAAGTTTGCCCGTGAAGTCGCCCGGCGAGTc                                       >  2:836232/1‑149 (MQ=255)
                                                                                                               tttGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCAGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTcc                                      >  4:74013/1‑149 (MQ=255)
                                                                                                                     gcGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGtt                                >  4:504312/1‑149 (MQ=255)
                                                                                                                      cGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTc                               <  4:612416/149‑1 (MQ=255)
                                                                                                                       gTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCt                              <  1:522152/149‑1 (MQ=255)
                                                                                                                           tGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGt                                                                            >  4:820352/1‑99 (MQ=255)
                                                                                                                           tGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGt                                                                            <  3:820352/99‑1 (MQ=255)
                                                                                                                             ccgccTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTAt                        >  4:165154/1‑149 (MQ=255)
                                                                                                                              cgccTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATc                       <  3:562892/149‑1 (MQ=255)
                                                                                                                                  tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATc                       >  2:1022417/1‑145 (MQ=255)
                                                                                                                                  tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATc                       <  1:1022417/145‑1 (MQ=255)
                                                                                                                                  tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt                   >  3:197192/1‑149 (MQ=255)
                                                                                                                                  tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt                   <  4:533067/149‑1 (MQ=255)
                                                                                                                                  tCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGtt                   <  2:338759/149‑1 (MQ=255)
                                                                                                                                       cagcatgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCa              >  1:527818/7‑149 (MQ=255)
                                                                                                                                        gcAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCAt             >  3:611730/1‑149 (MQ=255)
                                                                                                                                        gcAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCAt             >  3:491470/1‑149 (MQ=255)
                                                                                                                                         cAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCACCGGTTCTGATTATCCGTTGCTCAtt            <  1:705567/149‑1 (MQ=255)
                                                                                                                                          catgACGCCCTGTAC‑GGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTg           >  4:957675/4‑149 (MQ=255)
                                                                                                                                          aaaGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTg           >  3:626723/1‑149 (MQ=255)
                                                                                                                                                gCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg     >  1:528547/1‑149 (MQ=255)
                                                                                                                                                gCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg     >  1:610206/1‑149 (MQ=255)
                                                                                                                                                gCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcg     >  3:32998/1‑149 (MQ=255)
                                                                                                                                                 ccc‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTgcgc    <  3:642835/149‑1 (MQ=255)
                                                                                                                                                   c‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCt  <  1:278241/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CATGGTGAAAATGTGGGGTGACGGTTTCTGGAACTTGCTGGCGTTTGGTATGCAGATGGCGCTTATCATCGTTACCGGTCATGCCCTTGCCAGCTCTGCTCCGGTGAAAAGTTTGCTGCGTACTGCCGCCTCCGCCGCAAAGACGCCC‑GTACAGGGCGTCATGCTGGTCACTTTCTTCGGTTCAGTCGCTTGTGTCATCAACTGGGGATTTGGTTTGGTTGTCGGCGCAATGTTTGCCCGTGAAGTCGCCCGGCGAGTCCCCGGTTCTGATTATCCGTTGCTCATTGCCTGCGCCT  >  NC_000913/2324890‑2325185

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