Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,101,563 T→C 6.3% H41R (CAT→CGT)  csgG ← curli production assembly/transport outer membrane lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,101,5630TC 6.3% 283.7 / 6.2 94H41R (CAT→CGT) csgGcurli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base C (2/4);  total (42/52)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.14e-01

AATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGC  >  NC_000913/1101417‑1101710
                                                                                                                                                  |                                                                                                                                                   
aaTCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATgg                                                                                                                                                   <  1:1058333/149‑1 (MQ=255)
  tCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTATATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAGATGGGt                                                                                                                                                 <  2:1142691/149‑1 (MQ=255)
     ttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGAT‑GGTCaaa                                                                                                                                             <  2:854786/149‑1 (MQ=255)
          tGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCtt                                                                                                                                         <  2:1318344/149‑1 (MQ=37)
           gCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCttt                                                                                                                                        <  1:530566/149‑1 (MQ=37)
            ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                                       >  1:796738/1‑149 (MQ=255)
            ccGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTGTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTg                                                                                                                                       >  1:791461/1‑149 (MQ=255)
            ccGTGACCAGCATTGCCGTGGCGCTTTGAGTATCAGAATTTGAGAAGTTACTTGCCGTGTAGGGTTTAAATTGCCCGTTTTCGTACTGAATGTTGTATACCTAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTg                                                                                                                                       <  2:1133068/149‑1 (MQ=16)
               tGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAg                                                                                                                                    >  1:888215/1‑149 (MQ=255)
                 aCCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGct                                                                                                                                  >  2:1310940/1‑149 (MQ=255)
                  ccAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGtttcgg                                                                                                                                 <  2:316443/149‑6 (MQ=25)
                   cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                >  2:825497/1‑149 (MQ=255)
                   cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                >  2:238807/1‑149 (MQ=255)
                   cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCACGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                                >  1:1440911/1‑149 (MQ=255)
                     gCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGa                                                                                                                              <  2:701210/149‑1 (MQ=255)
                           ccGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGa                                                                                                                        >  2:573542/1‑149 (MQ=255)
                            cGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                       <  1:1494849/149‑1 (MQ=255)
                              tGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGTGTTAAAATTGCCCGGATTCGTCCTGAATGTTGTATACCGAAACAAAGAATTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCAGAGCACGAGGc                                                                                                                     >  1:1436090/1‑149 (MQ=255)
                                 cgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCAtt                                                                                                                  <  2:536470/149‑1 (MQ=255)
                                 cgcTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGGTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCAtt                                                                                                                  <  2:1289931/149‑1 (MQ=255)
                                  gcTTTGCGGAACAGCAGTGGAGAAGTTACTTGTCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTa                                                                                                                 >  1:575127/1‑149 (MQ=255)
                                    tttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAt                                                                                                               <  1:927550/149‑1 (MQ=255)
                                    tttGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAt                                                                                                               >  1:440670/1‑149 (MQ=255)
                                      tGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGt                                                                                                             <  2:443652/149‑1 (MQ=255)
                                       gCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTc                                                                                                            <  2:1391983/149‑1 (MQ=255)
                                        cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                       >  2:1157460/1‑137 (MQ=255)
                                        cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAg                                                                                                                       <  1:1157460/137‑1 (MQ=255)
                                        cGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTATAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGTCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGAATTAATTTCg                                                                                                           >  1:966985/1‑149 (MQ=255)
                                                       taaGGTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACATAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGATCACGAGGCATTAATGTCGGTCTGGCGGcttctt                                                                                            <  1:967060/148‑1 (MQ=255)
                                                       gaAGTTACATGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGcttctt                                                                                            >  2:1532351/1‑149 (MQ=255)
                                                        aaGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                           <  2:1312584/149‑1 (MQ=255)
                                                        aaGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCttt                                                                                           >  2:738446/1‑149 (MQ=255)
                                                         aGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTa                                                                                          <  1:665018/149‑1 (MQ=255)
                                                          gTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAg                                                                                         >  1:1242716/1‑149 (MQ=255)
                                                           ttACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTTAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAgg                                                                                        >  1:1110498/1‑149 (MQ=255)
                                                              cTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcgg                                                                                     <  1:899948/149‑1 (MQ=255)
                                                               ttGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAggcggg                                                                                    >  2:50406/1‑149 (MQ=255)
                                                                    gtgtattgTTTAAATTTCCCGGTTTCGTCATGAATGTTGTACACCTAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGCGNNATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGt                                                                               <  2:705264/142‑1 (MQ=255)
                                                                      gTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTa                                                                             >  2:1275379/1‑149 (MQ=255)
                                                                       tAGTGTGTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTAACCGTCGGAGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:295789/149‑1 (MQ=255)
                                                                       tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:764512/149‑1 (MQ=255)
                                                                       tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:415099/149‑1 (MQ=255)
                                                                       tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:1519687/149‑1 (MQ=255)
                                                                       tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:672084/149‑1 (MQ=255)
                                                                       tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:573542/149‑1 (MQ=255)
                                                                       tAGGGGTTAAATTGCCCGGTGTCGTCCTGACTGTTGTATACCGAAACAAAGATTTTAACCGTCGGCGCTGGCAGATGTGTCAAATCTTTGTAGCTCTGAGCAAGAGTCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  1:262551/149‑1 (MQ=255)
                                                                           gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                        >  2:293570/1‑149 (MQ=255)
                                                                           gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                        <  2:888215/149‑1 (MQ=255)
                                                                           gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                        >  2:709979/1‑149 (MQ=255)
                                                                            tttAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATc                                                                       >  1:197521/1‑149 (MQ=255)
                                                                             ttAAATTGCCCGGTTTCGTCCAGAATGTTGTATTCCGAAACAAAGATTTTACCAGTCGGCTCTTGCAGATGGGTCAAATCTTAGTAGCTCTGAGCACGAGGCATTAATGTAGGTCTGGCGGCTTCTTTAGGCGGGGCTGTTTAGAATcc                                                                      >  1:1390330/1‑149 (MQ=255)
                                                                               aaaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGc                                                                    <  2:1357662/149‑1 (MQ=255)
                                                                                aaTTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCCTCCGCt                                                                   >  1:1075117/1‑149 (MQ=255)
                                                                                      ccGGTTTCGCCCTGAATGTTGTATACCGAAACAAAGATTATACCCGGCGGCGCTGGCAGATGGGTCTAATCTTTGTAGATCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTaa                                                             <  1:300232/149‑1 (MQ=255)
                                                                                             cGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACg                                                      <  2:1045149/149‑1 (MQ=255)
                                                                                              gTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACTAGGCATTAATGTCGGTCTGGCGCCTTCTTTAGGCGGGCCGGTTAAGCATCCGCTCAGTAACATGACgg                                                     >  2:1230366/1‑149 (MQ=255)
                                                                                              gTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACgg                                                     >  1:1335333/1‑149 (MQ=255)
                                                                                                 cTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCaa                                                  >  1:38479/1‑149 (MQ=255)
                                                                                                  tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                 >  1:1007928/1‑149 (MQ=255)
                                                                                                  tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                 >  2:249954/1‑149 (MQ=255)
                                                                                                  tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                 <  2:796738/149‑1 (MQ=255)
                                                                                                  tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                                 <  2:791461/149‑1 (MQ=255)
                                                                                                     aTGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCCCGAGGCAGTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGCCGGCAACCaa                                              >  2:190841/1‑149 (MQ=255)
                                                                                                           tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            <  2:138511/113‑1 (MQ=255)
                                                                                                           tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                            >  1:138511/1‑113 (MQ=255)
                                                                                                            ataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAAt                                       <  1:322013/149‑1 (MQ=255)
                                                                                                              aCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAATTCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTTAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAATTaa                                     <  2:1451993/149‑1 (MQ=255)
                                                                                                                cGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATTTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAgc                                   <  1:238807/149‑1 (MQ=255)
                                                                                                                   taCAAAGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCt                                >  2:548515/2‑149 (MQ=37)
                                                                                                                    aCAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTg                               <  2:1563349/149‑1 (MQ=255)
                                                                                                                    aCAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTg                               <  2:1569502/149‑1 (MQ=255)
                                                                                                                      aaaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTAGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                             >  1:46539/1‑149 (MQ=255)
                                                                                                                      aaaGATTTGACCCATCTGCCAGCGCCGACGGGTAAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCa                             >  1:525367/1‑149 (MQ=37)
                                                                                                                       aaGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCAt                            >  1:1221277/1‑149 (MQ=255)
                                                                                                                              tACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGattatt                     <  1:577841/149‑1 (MQ=255)
                                                                                                                               aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt                    <  2:147391/149‑1 (MQ=255)
                                                                                                                               aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt                    <  1:825497/149‑1 (MQ=255)
                                                                                                                                accttCGGCGCTTGCAGATGGGTAAAATCTTTGTAGCTCTGAGCACGATGCATTAATGTCGGTATGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCATCCATAAGAAATAAGCGCTGCATGATTAtttt                   <  2:1390330/145‑1 (MQ=255)
                                                                                                                                 ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                  <  2:923131/149‑1 (MQ=255)
                                                                                                                                 ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                  <  2:1440911/149‑1 (MQ=255)
                                                                                                                                 ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                  <  2:38479/149‑1 (MQ=255)
                                                                                                                                 ccGTCGGCGCGGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc                  <  1:288962/149‑1 (MQ=255)
                                                                                                                                  cGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTcc                 >  1:1395016/1‑149 (MQ=255)
                                                                                                                                   gTCTGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGGCACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt                <  1:25947/149‑1 (MQ=255)
                                                                                                                                   gTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCGCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCt                <  2:1007928/149‑1 (MQ=255)
                                                                                                                                    tCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAAAGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCtt               <  1:332031/149‑1 (MQ=255)
                                                                                                                                           tGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGc        >  1:253703/1‑149 (MQ=255)
                                                                                                                                           tGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGc        >  1:1119161/1‑149 (MQ=255)
                                                                                                                                             gCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTg      >  2:418796/1‑149 (MQ=255)
                                                                                                                                              cAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTgg     <  1:1275379/149‑1 (MQ=255)
                                                                                                                                               cgATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTggg    <  2:1168510/148‑1 (MQ=255)
                                                                                                                                                gATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTgggg   >  1:1332062/1‑149 (MQ=255)
                                                                                                                                                 aTGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGc  >  2:762181/1‑149 (MQ=255)
                                                                                                                                                 aTGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGc  >  2:766533/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
AATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGCTGGGGC  >  NC_000913/1101417‑1101710

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