Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,292,003 G→A 5.1% S221S (TCG→TCA yraK → putative fimbrial‑like adhesin protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,292,0030GA 5.1% 312.6 / 4.2 97S221S (TCG→TCAyraKputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base G (44/48);  new base A (3/2);  total (47/50)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.90e-01

CCATTCCGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAA  >  NC_000913/3291876‑3292149
                                                                                                                               |                                                                                                                                                  
ccATTCCGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGAt                                                                                                                               >  6:237875/1‑149 (MQ=255)
 cATTCCGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATg                                                                                                                              >  5:192277/1‑149 (MQ=255)
     ccGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGtta                                                                                                                           <  2:404951/148‑1 (MQ=255)
      cGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGttat                                                                                                                          >  1:382049/1‑148 (MQ=255)
      cGAAAACGGATATTCCCGTCATTTAAGGCTCGTGGGTCGATTGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGttatt                                                                                                                         <  4:214681/149‑1 (MQ=255)
        aaaaCGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTaa                                                                                                                       <  7:371199/149‑1 (MQ=255)
         aaaCGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAg                                                                                                                      <  4:1576/149‑1 (MQ=255)
         aaaCGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAg                                                                                                                      >  7:36635/1‑149 (MQ=255)
          aaCGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAgg                                                                                                                     <  6:192277/149‑1 (MQ=255)
             ggATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGtt                                                                                                                   >  8:109772/1‑148 (MQ=255)
              gATATTGCCGTCATTAAAGGCGCGGGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGTCGATGTTATTAAGGTTaa                                                                                                                 <  2:266791/149‑1 (MQ=255)
               atatTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGc                                                                                                                                  >  7:341460/1‑131 (MQ=255)
               atatTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGc                                                                                                                                  <  8:341460/131‑1 (MQ=255)
                 atTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAttt                                                                                                              >  7:4603/1‑149 (MQ=255)
                  ttGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAtttt                                                                                                             <  8:36635/149‑1 (MQ=255)
                    gCCGTCCTTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTg                                                                                                            <  1:96300/148‑1 (MQ=255)
                    gCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTg                                                                                                            <  3:237706/148‑1 (MQ=255)
                     ccGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGa                                                                                                          <  3:319297/149‑1 (MQ=255)
                     ccGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGa                                                                                                          <  1:198733/149‑1 (MQ=255)
                     ccGGCATTAACGGCGCGTGGGGCGATGGAAGGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGGTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGa                                                                                                          <  5:349863/149‑1 (MQ=255)
                          aTTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATtcat                                                                                                     <  3:100578/149‑1 (MQ=255)
                          aTTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATtca                                                                                                      >  3:172433/1‑148 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCa                                                                                                                                      <  5:159705/115‑1 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCa                                                                                                                                      >  6:159705/1‑115 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAtttt                                                                                                             >  1:384543/1‑140 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAtttt                                                                                                             <  2:384543/140‑1 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATtcat                                                                                                     >  4:219109/1‑148 (MQ=255)
                           ttAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATtcat                                                                                                     >  1:228416/1‑148 (MQ=255)
                              aaGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCaa                                                                                                  <  4:45095/148‑1 (MQ=255)
                                ggCGCGTGGGTCGATGGAATGGGAAGCCAGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCAAAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATgg                                                                                               <  6:320787/149‑1 (MQ=255)
                                         gTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGtt                                                                                      <  4:208577/149‑1 (MQ=255)
                                           cGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTa                                                                                    <  2:382049/149‑1 (MQ=255)
                                             aTGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTac                                                                                   >  8:396554/1‑148 (MQ=255)
                                             aTGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTac                                                                                   >  5:146211/1‑148 (MQ=255)
                                               ggAATGGCAAGCCCGTCGACAGGTGACTTCCATGATTTACTCAAGTTATCGATTCAGGGAAATCTCCCCGCCCCACAGTCGAGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGACATCAGTTTacgca                                                                                >  4:4865/1‑146 (MQ=255)
                                                 aaTGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCAccc                                                                              >  4:225934/1‑149 (MQ=255)
                                                      gAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCAACCGCAat                                                                         >  6:356522/1‑149 (MQ=255)
                                                           ccGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAatgccat                                                                    <  8:4603/149‑1 (MQ=255)
                                                           ccGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAatgccat                                                                    >  7:289147/1‑149 (MQ=255)
                                                            cGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAatgccatg                                                                   <  8:289147/149‑1 (MQ=255)
                                                            cGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACAGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAtt                                                                                                               <  5:5755/105‑1 (MQ=255)
                                                            cGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACAGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAAtt                                                                                                               >  6:5755/1‑105 (MQ=255)
                                                             gTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCATGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAatgccatgc                                                                  >  2:136661/1‑149 (MQ=255)
                                                                  cAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCaaaa                                                                                                                                             >  3:348395/1‑69 (MQ=255)
                                                                  cAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCaaaa                                                                                                                                             <  4:348395/69‑1 (MQ=255)
                                                                  cAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAAtgccatgcca                                                                <  7:396554/146‑1 (MQ=255)
                                                                    ggTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACgg                                                           >  3:362978/1‑149 (MQ=255)
                                                                     gTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGt                                                          >  3:204839/1‑149 (MQ=255)
                                                                        aCTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGttt                                                        >  1:89536/1‑148 (MQ=255)
                                                                         cTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTa                                                      <  5:237875/149‑1 (MQ=255)
                                                                         cTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTa                                                      <  2:228416/149‑1 (MQ=255)
                                                                          ttCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCCCCCGCAATGCACTGCCAGACGGTTTTAc                                                     >  5:364481/1‑149 (MQ=255)
                                                                          ttCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTAc                                                     >  7:286336/1‑149 (MQ=255)
                                                                          ttCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGAAAAATTAATCAGGGCGATGTTAGTAAGGTTAATTTTGGATTCATCAATGGACAGAAGTTTTCCACCCGCAATGCCAGGCCAGACGCTTTTAc                                                     >  2:267212/1‑149 (MQ=255)
                                                                          ttCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCATGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTAc                                                     >  2:392592/1‑149 (MQ=255)
                                                                          ttCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCATGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTAc                                                     >  8:4299/1‑149 (MQ=255)
                                                                             cATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTc                                                   <  8:286336/148‑1 (MQ=255)
                                                                             cATGAGTTAGTCAATTTATCGATTCAGGGAAATCTCACCGCCCAACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTc                                                   <  4:489/148‑1 (MQ=255)
                                                                              aTGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTc                                                   >  1:27275/1‑147 (MQ=255)
                                                                                gATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAg                                                >  7:380217/1‑148 (MQ=255)
                                                                                gATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTAATCCAg                                                >  1:246832/1‑148 (MQ=255)
                                                                                 aTTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAg                                                >  2:212818/1‑147 (MQ=255)
                                                                                       tCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACttt                                        <  4:362978/149‑1 (MQ=255)
                                                                                       tCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACttt                                        >  2:79017/1‑149 (MQ=255)
                                                                                            ttATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGatat                                   >  2:344065/1‑149 (MQ=255)
                                                                                               tCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCAc                                >  1:272171/1‑149 (MQ=255)
                                                                                               tCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCAc                                <  7:266046/149‑1 (MQ=255)
                                                                                                   ttCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTAt                            >  4:120567/1‑149 (MQ=255)
                                                                                                      aGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGa                          <  7:218240/148‑1 (MQ=255)
                                                                                                      aGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGAc                         >  3:369900/1‑149 (MQ=255)
                                                                                                       gggAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCAccc                                                                              <  8:232581/95‑1 (MQ=255)
                                                                                                       gggAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCAccc                                                                              >  7:232581/1‑95 (MQ=255)
                                                                                                       gggAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACt                        >  3:15335/1‑149 (MQ=255)
                                                                                                          aaaTCTCACCGCCGCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTTAGGTTAATTTTGGATTCATCAATGTTCAGAAGTTTACCACCCGCAATGACATGCCAGACGGTTTTACTCCAGTTGACTTTGATATCACTTATGACtgtg                     <  3:4865/149‑1 (MQ=255)
                                                                                                            aTCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTgtggt                   <  7:53728/149‑1 (MQ=255)
                                                                                                             tctcACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTgtggt                   >  5:243196/1‑148 (MQ=255)
                                                                                                               tcACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGa                 >  3:278476/1‑148 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCGTGCAAAATTAATCAGGGTGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACtg                       <  7:158904/140‑1 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCGTGCAAAATTAATCAGGGTGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACtg                       >  8:158904/1‑140 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGa                 <  7:109772/146‑1 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATa               <  5:356522/148‑1 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATAc              >  7:255792/1‑149 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCATGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGa                 <  7:4299/146‑1 (MQ=255)
                                                                                                                 aCCGCCCCACAGTCATGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGa                 <  1:392592/146‑1 (MQ=255)
                                                                                                                   cGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACtt            <  4:172433/149‑1 (MQ=255)
                                                                                                                    gCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTc           <  4:15335/149‑1 (MQ=255)
                                                                                                                    gCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTGATGCCTGTGGTGATACtt            >  1:210420/1‑148 (MQ=255)
                                                                                                                    gCCCCACAGCCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTc           <  2:422495/149‑1 (MQ=255)
                                                                                                                      cccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCa          >  8:28324/1‑148 (MQ=255)
                                                                                                                       ccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGAc                                           <  2:255584/114‑1 (MQ=255)
                                                                                                                       ccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGAc                                           >  1:255584/1‑114 (MQ=255)
                                                                                                                       ccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGTTACTTCaaa        <  7:410347/149‑1 (MQ=255)
                                                                                                                       ccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCaaa        <  8:156596/149‑1 (MQ=255)
                                                                                                                       ccACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCaaa        <  1:79017/149‑1 (MQ=255)
                                                                                                                         acaGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGa      <  1:344065/149‑1 (MQ=255)
                                                                                                                             tCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaa  <  8:380217/149‑1 (MQ=255)
                                                                                                                             tCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTaa  <  1:212818/149‑1 (MQ=255)
                                                                                                                               |                                                                                                                                                  
CCATTCCGAAAACGGATATTGCCGTCATTAAAGGCGCGTGGGTCGATGGAATGGGAAGCCCGTCTACAGGTGACTTCCATGATTTAGTCAAGTTATCGATTCAGGGAAATCTCACCGCCCCACAGTCGTGCAAAATTAATCAGGGCGATGTTATTAAGGTTAATTTTGGATTCATCAATGGTCAGAAGTTTACCACCCGCAATGCCATGCCAGACGGTTTTACTCCAGTAGACTTTGATATCACTTATGACTGTGGTGATACTTCAAAGATTAA  >  NC_000913/3291876‑3292149

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