Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 302,227 A→T 6.0% L116H (CTT→CAT)  paoA ← PaoABC aldehyde oxidoreductase, 2Fe‑2S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913302,2270AT 6.0% 227.3 / 2.8 65L116H (CTT→CAT) paoAPaoABC aldehyde oxidoreductase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base A (33/28);  new base T (2/2);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCAC  >  NC_000913/302079‑302375
                                                                                                                                                    |                                                                                                                                                    
ttGCCCGGACGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTaa                                                                                                                                                      <  3:523406/149‑1 (MQ=255)
  gCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGc                                                                                                                                                    <  3:369667/149‑1 (MQ=255)
      ggAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGc                                                                                                                                                <  1:614905/149‑1 (MQ=255)
          gTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGCCAGGCTGGCATTAAGCCTGCGAcc                                                                                                                                            >  2:91956/1‑149 (MQ=255)
                tAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAc                                                                                                                                      <  4:97027/149‑1 (MQ=255)
                tAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAc                                                                                                                                      <  4:93850/149‑1 (MQ=255)
                tAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGAc                                                                                                                                      <  4:474694/149‑1 (MQ=255)
                         tGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGt                                                                                                                             >  1:218726/1‑149 (MQ=255)
                          ggAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCt                                                                                                                                                  <  3:291930/127‑1 (MQ=255)
                          ggAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCt                                                                                                                                                  >  4:291930/1‑127 (MQ=255)
                            aaGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACa                                                                                                                          <  4:220075/149‑1 (MQ=255)
                                 atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTc                                                                                                                     <  3:335914/149‑1 (MQ=255)
                                 atcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTc                                                                                                                     <  2:74441/149‑1 (MQ=255)
                                  tcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTcc                                                                                                                    <  4:323689/149‑1 (MQ=255)
                                  tcatGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTcc                                                                                                                    >  4:319855/1‑149 (MQ=255)
                                       cTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACt                                                                                                               >  4:139721/1‑149 (MQ=255)
                                       cTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACt                                                                                                               >  1:560178/1‑149 (MQ=255)
                                          gATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTc                                                                                                            >  2:206789/1‑149 (MQ=255)
                                             aaaGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGt                                                                                                         >  1:500558/1‑149 (MQ=255)
                                              aaGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTg                                                                                                        >  2:565625/1‑149 (MQ=255)
                                                ggCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCTTAATGCCTGcct                                                                                                                                              >  2:229250/1‑107 (MQ=255)
                                                ggCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGAt                                                                                                      >  4:100310/1‑149 (MQ=255)
                                                ggCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGAt                                                                                                      <  2:459954/149‑1 (MQ=255)
                                                ggCCGCCTGCATGGAGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCTTAATGCCTGcct                                                                                                                                              <  1:229250/109‑3 (MQ=255)
                                                 gccgccTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATc                                                                                                     <  3:319855/149‑1 (MQ=255)
                                                      cTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAAc                                                                                                >  1:15059/1‑149 (MQ=255)
                                                        gCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCt                                                                                              <  2:87365/149‑1 (MQ=255)
                                                              ggTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTg                                                                                        <  2:560178/149‑1 (MQ=255)
                                                               gTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTgg                                                                                       <  2:212785/149‑1 (MQ=255)
                                                                tGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                      >  1:384203/1‑149 (MQ=255)
                                                                tGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGt                                                                                      <  1:2634/149‑1 (MQ=255)
                                                                   aGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTAcc                                                                                   <  3:139721/149‑1 (MQ=255)
                                                                    gATTATCTGGCGAACCCAGGCCTTCAATGTTGGTGATCTCGGCCCCCTGATGCATAACTGCAAGCGTCAGGAAGGCATTAAGCATGCGACCATTGACCAGCACGGTACAGGCTACGCACTGACTGTGATCGCAACCTTTCTTGGTACCg                                                                                  >  2:598724/1‑149 (MQ=255)
                                                                      ttATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGAt                                                                                <  1:565625/149‑1 (MQ=255)
                                                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAtt                                                                                                                                         >  4:344071/1‑83 (MQ=255)
                                                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCAtt                                                                                                                                         <  3:344071/83‑1 (MQ=255)
                                                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAg                                                                       >  4:174538/1‑149 (MQ=255)
                                                                               gAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAg                                                                       >  2:100310/1‑149 (MQ=255)
                                                                                  cccAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAtt                                                                    >  2:154646/1‑149 (MQ=255)
                                                                                   ccAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGAttt                                                                   >  2:186608/1‑149 (MQ=255)
                                                                                       gCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCAc                                                               <  1:206789/149‑1 (MQ=255)
                                                                                          ttCAATGGTGGTGAGCTCGGCCCCCTGATGCATGACTGCGAGCGTAAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCa                                                            >  4:75677/1‑149 (MQ=255)
                                                                                             aaTGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCaaa                                                                            <  4:312857/130‑1 (MQ=255)
                                                                                             aaTGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCaaa                                                                            >  3:312857/1‑130 (MQ=255)
                                                                                                 gtggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCaaa                                                                            <  1:584902/126‑1 (MQ=255)
                                                                                                 gtggtgATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCaaa                                                                            >  2:584902/1‑126 (MQ=255)
                                                                                                      gATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTg                                                                                                        >  3:342880/1‑93 (MQ=255)
                                                                                                      gATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTg                                                                                                        <  4:342880/93‑1 (MQ=255)
                                                                                                        tctcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGc                                              >  4:470867/1‑149 (MQ=255)
                                                                                                        tctcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGc                                              <  2:218726/149‑1 (MQ=255)
                                                                                                        tctcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAACTGTCCAGTAGc                                              >  3:290494/1‑149 (MQ=255)
                                                                                                          tcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGt                                            >  3:138165/1‑149 (MQ=255)
                                                                                                          tcGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGt                                            >  2:525371/1‑149 (MQ=255)
                                                                                                                   tGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGt                                   <  2:500558/149‑1 (MQ=255)
                                                                                                                    gATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTa                                  <  1:358036/149‑1 (MQ=255)
                                                                                                                     aTGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTAt                                 >  4:55242/1‑149 (MQ=255)
                                                                                                                        cATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCa                              >  1:214132/1‑149 (MQ=255)
                                                                                                                         aTGACTGCAAGCGTCAAGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCAc                             >  2:156525/1‑149 (MQ=255)
                                                                                                                           gACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCt                           >  3:434338/1‑149 (MQ=255)
                                                                                                                                   gCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAgctgc                   >  1:5976/1‑149 (MQ=255)
                                                                                                                                      tCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCg                <  3:262967/149‑1 (MQ=255)
                                                                                                                                        tcgcaggctTAATGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGt              >  3:453579/3‑149 (MQ=255)
                                                                                                                                         ggcaggcaTTAAGCCTGCGACCATTGACCGGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGtt             >  3:369507/1‑149 (MQ=255)
                                                                                                                                         cgcaggctTAATGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATc                                                                               <  3:1029/82‑1 (MQ=255)
                                                                                                                                                 ttAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCggg                                    <  1:176520/118‑1 (MQ=255)
                                                                                                                                                 ttAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGAGTCCAGTAGCGTGGTTCggg                                    >  2:176520/1‑118 (MQ=255)
                                                                                                                                                    aGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCAc  >  4:385908/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTTTTGCCGTTCAC  >  NC_000913/302079‑302375

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