Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,817,948 G→T 5.2% E63D (GAG→GAT dinD → DNA damage‑inducible protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,817,9480GT 5.2% 1297.0 / 38.4 360E63D (GAG→GATdinDDNA damage‑inducible protein
Reads supporting (aligned to +/- strand):  ref base G (185/156);  new base T (5/14);  total (190/170)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 374 total aligned reads displayed.

AATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGA  >  NC_000913/3817802‑3818096
                                                                                                                                                  |                                                                                                                                                    
  tGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATcc                                                                                                                                                  >  3:80538/1‑149 (MQ=255)
     aaaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATcc                                                                                                                                                  <  3:680690/146‑1 (MQ=255)
     aaaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCg                                                                                                                                               >  1:257679/1‑149 (MQ=255)
      aaCGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGa                                                                                                                                              <  1:72395/149‑1 (MQ=255)
        cGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGATCTCCTCGaaa                                                                                                                                            <  4:280115/149‑3 (MQ=255)
        cGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGata                                                                                                                                            <  4:700081/149‑1 (MQ=255)
          gAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATg                                                                                                                                          <  1:377037/149‑1 (MQ=255)
           aGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATgg                                                                                                                                         <  3:370779/149‑1 (MQ=255)
             cAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTc                                                                                                                                       <  3:163608/149‑1 (MQ=255)
             cAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATATCGAGGATCTCCTCGAAATGGTc                                                                                                                                       <  1:517941/149‑1 (MQ=255)
                 aCAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCaaaa                                                                                                                                   <  1:106627/149‑1 (MQ=255)
                      ggTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGc                                                                                                                              <  4:113612/149‑1 (MQ=255)
                        tCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTc                                                                                                                            >  3:453455/1‑149 (MQ=255)
                          aGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCaa                                                                                                                          >  4:236906/1‑149 (MQ=255)
                            cAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAAt                                                                                                                        <  1:66993/149‑1 (MQ=255)
                               gACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGcc                                                                                                                     <  3:622939/149‑1 (MQ=255)
                                   aCTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCaaaa                                                                                                                 <  3:343133/149‑1 (MQ=255)
                                     ttGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                               >  1:562194/1‑149 (MQ=255)
                                     ttGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                               >  1:530522/1‑149 (MQ=255)
                                     ttGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                               >  1:694689/1‑149 (MQ=255)
                                     ttGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAaga                                                                                                               >  1:194076/1‑149 (MQ=255)
                                      tGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAagag                                                                                                              <  2:87864/149‑1 (MQ=255)
                                       gggAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGc                                                                                                             >  3:73937/1‑149 (MQ=255)
                                           aaCTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTa                                                                                                         >  1:297339/1‑149 (MQ=255)
                                              tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                      <  2:294685/149‑1 (MQ=255)
                                              tactGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTaaaa                                                                                                      >  1:15810/1‑149 (MQ=255)
                                               actGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAg                                                                                                     <  4:269197/149‑1 (MQ=255)
                                                ctGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGa                                                                                                    <  1:359431/149‑1 (MQ=255)
                                                ctGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGa                                                                                                    <  2:149584/149‑1 (MQ=255)
                                                   ggTTATTCAGAGTATCGTCACTTTATACCTGTACTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACAt                                                                                                 >  3:213718/1‑149 (MQ=255)
                                                     ttATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCg                                                                                               >  3:690932/1‑149 (MQ=255)
                                                        ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAc                                                                                            >  2:667765/1‑149 (MQ=255)
                                                        ttCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTAc                                                                                            >  1:581206/1‑149 (MQ=255)
                                                            gagTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACtctc                                                                                        >  2:571926/1‑149 (MQ=255)
                                                              gTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTccc                                                                                      <  2:318730/149‑1 (MQ=255)
                                                              gTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTccc                                                                                      <  2:340092/149‑1 (MQ=255)
                                                               tATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCg                                                                                     >  3:266727/1‑149 (MQ=255)
                                                                    tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATg                                                                                >  3:645604/1‑149 (MQ=255)
                                                                    tCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATg                                                                                >  2:226836/1‑149 (MQ=255)
                                                                     cgcTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGc                                                                               <  2:372384/147‑1 (MQ=255)
                                                                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCt                                                                             <  3:194318/149‑1 (MQ=255)
                                                                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCt                                                                             >  4:175113/1‑149 (MQ=255)
                                                                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCt                                                                             <  1:409136/149‑1 (MQ=255)
                                                                       cTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCt                                                                             <  1:750217/149‑1 (MQ=255)
                                                                        tttATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTg                                                                            <  3:406113/149‑1 (MQ=255)
                                                                          tataCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGtt                                                                          <  4:358376/149‑1 (MQ=255)
                                                                             aCCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTAcc                                                                       >  1:386248/1‑149 (MQ=255)
                                                                              ccTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCt                                                                      >  3:425508/1‑149 (MQ=255)
                                                                              ccTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCt                                                                      <  2:131247/149‑1 (MQ=255)
                                                                               cTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTgg                                                                                                                               <  1:649342/91‑1 (MQ=255)
                                                                                tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                    <  2:470765/149‑1 (MQ=255)
                                                                                tGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTgg                                                                    >  4:576055/1‑149 (MQ=255)
                                                                                 gTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgt                                                                   >  1:174339/1‑149 (MQ=255)
                                                                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGcc                                                                                                                     >  3:220055/1‑98 (MQ=255)
                                                                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGcc                                                                                                                     <  4:220055/98‑1 (MQ=255)
                                                                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgta                                                                  >  2:411094/1‑149 (MQ=255)
                                                                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgta                                                                  <  4:433115/149‑1 (MQ=255)
                                                                                  tATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgta                                                                  >  4:45258/1‑149 (MQ=255)
                                                                                    ttAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagt                                                                >  4:601404/1‑149 (MQ=255)
                                                                                     tAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGgtagta                                                               >  4:684893/1‑149 (MQ=255)
                                                                                           gcgcCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAAc                                                         <  4:734129/149‑1 (MQ=255)
                                                                                              ccAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                                      <  4:328964/149‑1 (MQ=255)
                                                                                               cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCg                                                     >  4:52138/1‑149 (MQ=255)
                                                                                               cAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCg                                                     >  1:624256/1‑149 (MQ=255)
                                                                                                  agaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCa                                                                                                            >  2:190067/1‑91 (MQ=255)
                                                                                                  agaagCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGAcc                                                  >  4:395553/1‑149 (MQ=255)
                                                                                                      gCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCGGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGc                                              <  1:225215/149‑1 (MQ=255)
                                                                                                        cTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGa                                            >  3:305065/1‑149 (MQ=255)
                                                                                                        cTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGa                                            >  1:325366/1‑149 (MQ=255)
                                                                                                             agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                                       <  2:64192/149‑1 (MQ=255)
                                                                                                             agaACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCg                                       >  2:545469/1‑149 (MQ=255)
                                                                                                               aaCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGc                                                      >  2:517857/1‑132 (MQ=255)
                                                                                                                aCAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTc                                    <  1:612951/149‑1 (MQ=255)
                                                                                                                 cAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATCGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCa                                   >  1:498252/1‑149 (MQ=255)
                                                                                                                   gTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAAcc                                        >  4:79367/1‑142 (MQ=255)
                                                                                                                   gTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCAtt                                 >  4:474961/1‑149 (MQ=255)
                                                                                                                    tGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTg                                <  2:499585/149‑1 (MQ=255)
                                                                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCgg                             <  1:538943/149‑1 (MQ=255)
                                                                                                                       tCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCgg                             <  4:62945/149‑1 (MQ=255)
                                                                                                                         acacaATTGATGACCATTTCGAGGATATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggca                           <  2:696153/149‑1 (MQ=255)
                                                                                                                           acaATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagg                         >  2:706966/1‑149 (MQ=255)
                                                                                                                             aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggc                       >  4:208800/1‑149 (MQ=255)
                                                                                                                             aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggc                       >  1:456260/1‑149 (MQ=255)
                                                                                                                             aaTTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagggc                       >  2:611616/1‑149 (MQ=255)
                                                                                                                                tgatgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGa                    <  3:247358/149‑1 (MQ=255)
                                                                                                                                tgatgaCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGa                    <  2:171314/149‑1 (MQ=255)
                                                                                                                                     aCCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGa                    <  3:729444/144‑1 (MQ=255)
                                                                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagg                         >  2:765838/1‑138 (MQ=255)
                                                                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCggcagg                         <  1:765838/138‑1 (MQ=255)
                                                                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt              >  3:575961/1‑149 (MQ=255)
                                                                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt              >  3:545903/1‑149 (MQ=255)
                                                                                                                                      ccATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtt              >  2:263351/1‑149 (MQ=255)
                                                                                                                                        aTTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTAtttt            >  1:137533/1‑149 (MQ=255)
                                                                                                                                         tttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg           >  3:128186/1‑149 (MQ=255)
                                                                                                                                         tttCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTg           <  2:233469/149‑1 (MQ=255)
                                                                                                                                         tatCGAGGATCTCCTCGAAATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATc                                                                                                <  3:596810/62‑1 (MQ=25)
                                                                                                                                              aggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATc      <  4:13822/149‑1 (MQ=255)
                                                                                                                                               ggagATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAACAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGTTAGTACAAATCGGCGACCCTGCTAAACCGGTCATTGCGGCAGGGCAGTCTTATTTTGCTATcc     <  1:99435/149‑1 (MQ=255)
                                                                                                                                                gagaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCa    <  2:748489/149‑1 (MQ=255)
                                                                                                                                                  gaTCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGa  <  1:734074/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AATGCAAACGGAGCAGAACAATGGTCAGCAAGACAACTTGGGAAACTACTGGGTTATTCAGAGTATCGTCACTTTATACCTGTATTAACGCGCGCCAAAGAAGCCTGTGAGAACAGTGGTCACACAATTGATGACCATTTCGAGGAGATCCTCGATATGGTCAAAATTGGCTCAAATGCCAAAAGAGCATTAAAAGACATCGTACTCTCCCGCTATGCCTGTTACCTGGTAGTACAAAACGGCGACCCTGCGAAACCGGTCATTGCGGCAGGGCAGACTTATTTTGCTATCCAGA  >  NC_000913/3817802‑3818096

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