Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 302,225 T→A 6.2% N117Y (AAT→TAT)  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,2250TA 6.2% 220.1 / 4.7 64N117Y (AAT→TAT) paoAPaoABC aldehyde oxidoreductase, 2Fe‑2S subunit
Reads supporting (aligned to +/- strand):  ref base T (32/28);  new base A (2/2);  total (34/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 = 9.49e-01

TTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTT  >  NC_000913/302079‑302364
                                                                                                                                                  |                                                                                                                                           
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)
                                                                                                                                                  |                                                                                                                                           
TTGCCCGGAGGTGCAGTAGCCGCACTGGAAGCCATCATGCTTGATAAAGGCCGCCTGCATGGGGTGAAGATTATCTGGCGAACCCAGGCCTTCAATGGTGGTGATCTCGGCCCCCTGATGCATGACTGCAAGCGTCAGGCAGGCATTAAGCCTGCGACCATTGACCAGCACGGTACAGGCTCCGCACTGTCCGTGATCGCAACCTTTCTTGGTACCGATCAAATGCAGATTTTCACGCAAAGTGTCCAGTAGCGTGGTTCGGGTATCCACCTCAAGCTGCTCGGTT  >  NC_000913/302079‑302364

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: ‑