Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,104,249 A→T 5.1% N100I (AAC→ATC)  csgB → curlin nucleator protein, minor subunit in curli complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,104,2490AT 5.1% 242.4 / 3.3 78N100I (AAC→ATC) csgBcurlin nucleator protein, minor subunit in curli complex
Reads supporting (aligned to +/- strand):  ref base A (33/41);  new base T (2/2);  total (35/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

CTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGT  >  NC_000913/1104105‑1104393
                                                                                                                                                |                                                                                                                                                
cTGGGACTAATACTAGTGCTCAGTGACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGAt                                                                                                                                              <  6:3841/149‑1 (MQ=255)
  gggACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGc                                                                                                                                            >  4:171465/1‑149 (MQ=255)
    gACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCa                                                                                                                                          <  1:147344/149‑1 (MQ=255)
         ataataGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAtt                                                                                                                                     <  2:396310/149‑1 (MQ=255)
         ataataGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATACATTGATCAGGCGGGCAGTGCCAACGATGCCAGTAtt                                                                                                                                     >  7:80555/1‑149 (MQ=255)
            ataGTGCTCAGTTACGGTAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTc                                                                                                                                   <  4:200153/148‑1 (MQ=255)
            ataGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCg                                                                                                                                  <  8:229979/149‑1 (MQ=255)
             taGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGTCGGGCAGTGCCAACGATGCCAGTAtttag                                                                                                                                  >  3:199716/1‑146 (MQ=255)
              aGTGCTCAGTTACGGCAGGGAGGCTCAAAACTATTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACATTGCATAAATTGATCAGGCGGGCAGTGACAACGATGCCAGAATTTCGc                                                                                                                                 >  1:114475/1‑148 (MQ=255)
               gTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTACAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCaa                                                                                                                               <  6:66216/149‑1 (MQ=255)
                tGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAg                                                                                                                              <  4:41065/149‑1 (MQ=255)
                  cTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCATCGATGCCAGTATTTCGCAAGGt                                                                                                                            >  3:169558/1‑149 (MQ=255)
                  cTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCATCGATGCCAGTATTTCGCAAGGt                                                                                                                            >  3:163981/1‑149 (MQ=255)
                  cTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGt                                                                                                                            >  5:250902/1‑149 (MQ=255)
                  cTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGt                                                                                                                            >  8:262109/1‑149 (MQ=255)
                   tCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGAACAGACAGGACATTATAACCATGCATAAATTGATCAGGCGGGCAGTGCCAACGATGCAAGTATTTAGCAAGTTg                                                                                                                           >  6:112157/1‑149 (MQ=255)
                         aCGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTAt                                                                                                                      <  4:256385/148‑1 (MQ=255)
                          cGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATg                                                                                                                     <  5:298669/148‑1 (MQ=255)
                          aggCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGAGAAAGAAGGTAGTATCATAAGGGCATAGATTGACCAGACAGGTGTTTATTACCTTGGATATATTGATCAGGCGTGCATTGCCAACGATGCCAGTTTTTATCAAGGTGCTTAtt                                                                                                                     <  4:199716/147‑2 (MQ=255)
                             cAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTaa                                                                                                                 <  2:21749/149‑1 (MQ=255)
                             cAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTaa                                                                                                                 <  1:331067/149‑1 (MQ=255)
                                    gCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGc                                                                                                           <  6:250902/148‑1 (MQ=255)
                                      tCAAAACTTTTGGCGTTTGTTGCGCAAGAAGGTAGTAGCTACCGGGAAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATgg                                                                                                                    <  4:31979/137‑1 (MQ=255)
                                      tCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgat                                                                                                        >  6:203323/1‑149 (MQ=255)
                                      tCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgat                                                                                                        >  6:335766/1‑149 (MQ=255)
                                      tCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgat                                                                                                        >  3:234819/1‑149 (MQ=255)
                                      tCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACAAGACAGGAGATTATAACCTTGAATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATgg                                                                                                                    >  3:31979/1‑137 (MQ=255)
                                        aaaaCTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgat                                                                                                        <  2:262830/147‑1 (MQ=255)
                                        aaaaCTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgat                                                                                                        <  1:7112/147‑1 (MQ=255)
                                          aaCTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCGACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCgatgat                                                                                                     >  6:226934/1‑148 (MQ=255)
                                           aCTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTa                                                                                                   <  6:186780/149‑1 (MQ=255)
                                               ttGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGGGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATc                                                                                                 <  1:213792/147‑1 (MQ=255)
                                               ttGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATACATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATc                                                                                                 <  8:80555/147‑1 (MQ=255)
                                                   cGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGa                                                                                             >  7:323743/1‑147 (MQ=255)
                                                    ggTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATACCCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAg                                                                                          <  2:222804/149‑1 (MQ=255)
                                                         ttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTc                                                                                      <  1:198172/148‑1 (MQ=255)
                                                         ttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTc                                                                                      <  7:262109/148‑1 (MQ=255)
                                                         ttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCt                                                                                     >  1:152698/1‑149 (MQ=255)
                                                         ttgCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCt                                                                                     <  4:234819/149‑1 (MQ=255)
                                                            agcAAGAAGGTAGTAGCAACAGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTgg                                                                                   <  6:285420/147‑1 (MQ=255)
                                                               aagaagGTAGTAGCAACCGGGCAAAGATTGACCACACAGGAGATTATAACCTAGCATATATTGATCAGGCGGGCAGTGCCAACGATCCCAGTAATTCGCAAGGTGCTTATGGTAATACTGCGATGATAATCCAGAAAGGTTCTGGtaa                                                                                >  7:56874/1‑148 (MQ=255)
                                                                      tagtagCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGtaataa                                                                             <  1:162973/144‑1 (MQ=255)
                                                                      tagtagCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGtaataa                                                                             >  2:162973/1‑144 (MQ=255)
                                                                           gCAACCTGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGAAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATAtt                                                                   <  5:295179/149‑1 (MQ=255)
                                                                             aaCCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTa                                                                  >  5:103208/1‑148 (MQ=255)
                                                                              aCCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTaca                                                                <  5:83534/149‑1 (MQ=255)
                                                                              aCCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTaca                                                                <  5:221123/149‑1 (MQ=255)
                                                                                cGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTAcaca                                                              <  2:165776/149‑1 (MQ=255)
                                                                                     aaaGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATg                                                         >  7:257717/1‑149 (MQ=255)
                                                                                     aaaGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATgg                                                        <  4:315299/149‑1 (MQ=255)
                                                                                             cccagccagGAGATTATAACCTTGCAGATATTGATCAGGCGGGCAGTGCCAACGATGCCCGTATTACGCAAGGTGCTTATCGAAAAACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATAATACACAGTAAGCTACTCa                                                  >  6:100723/2‑148 (MQ=255)
                                                                                             aCCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCa                                                  >  2:156615/1‑148 (MQ=255)
                                                                                             aCCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCa                                                  >  1:53984/1‑148 (MQ=255)
                                                                                               cAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCATCGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCaaaa                                               <  4:163981/149‑1 (MQ=255)
                                                                                               cAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCATCGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCaaaa                                               <  4:169558/149‑1 (MQ=255)
                                                                                                  aCAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAgg                                                                                         >  3:262052/1‑104 (MQ=255)
                                                                                                  aCAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAgg                                                                                         <  4:262052/104‑1 (MQ=255)
                                                                                                    aGGAGATTATAACCTTGCATATATTGATCAGTCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACgg                                          <  2:296293/149‑1 (MQ=255)
                                                                                                     ggAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACgg                                          >  2:314624/1‑148 (MQ=255)
                                                                                                      gagaTTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGc                                         >  2:270755/1‑148 (MQ=255)
                                                                                                         aTTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAAtt                                     >  3:323912/1‑149 (MQ=255)
                                                                                                            ataACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTg                                    >  2:37584/1‑147 (MQ=255)
                                                                                                            ataACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTg                                    >  4:330804/1‑147 (MQ=255)
                                                                                                             taACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAg                                 >  6:369222/1‑149 (MQ=255)
                                                                                                               aCCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTg                               >  8:18531/1‑149 (MQ=255)
                                                                                                                 cTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTg                               >  5:380330/1‑147 (MQ=255)
                                                                                                                   tGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTa                                                      >  4:218419/1‑122 (MQ=255)
                                                                                                                    gCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCaga                           <  7:18531/148‑1 (MQ=255)
                                                                                                                     cATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAgaga                         <  3:330804/149‑1 (MQ=255)
                                                                                                                       tatatTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACa                       <  8:257717/149‑1 (MQ=255)
                                                                                                                         tatTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAg                      <  5:203323/148‑1 (MQ=255)
                                                                                                                          atTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACTGc                                         >  7:339473/1‑128 (MQ=255)
                                                                                                                          atTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGc                                         <  8:339473/128‑1 (MQ=255)
                                                                                                                           ttGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTc                    <  1:37584/148‑1 (MQ=255)
                                                                                                                           ttGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCg                   <  2:152698/149‑1 (MQ=255)
                                                                                                                                 aGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATg             >  2:96197/1‑149 (MQ=255)
                                                                                                                                  ggCGGGCAGTGCCAACGATTCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATgg            <  5:335766/149‑1 (MQ=255)
                                                                                                                                        cAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATCCTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTc      <  2:311865/149‑1 (MQ=255)
                                                                                                                                             ctaaCGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTTGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGt  >  4:239732/3‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
CTGGGACTAATAATAGTGCTCAGTTACGGCAGGGAGGCTCAAAACTTTTGGCGGTTGTTGCGCAAGAAGGTAGTAGCAACCGGGCAAAGATTGACCAGACAGGAGATTATAACCTTGCATATATTGATCAGGCGGGCAGTGCCAACGATGCCAGTATTTCGCAAGGTGCTTATGGTAATACTGCGATGATTATCCAGAAAGGTTCTGGTAATAAAGCAAATATTACACAGTATGGTACTCAAAAAACGGCAATTGTAGTGCAGAGACAGTCGCAAATGGCTATTCGCGT  >  NC_000913/1104105‑1104393

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