Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,101,556 T→G 7.1% P43P (CCA→CCC) ‡ 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,5560TG 7.1% 251.2 / 5.8 89P43P (CCA→CCC) ‡csgGcurli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/47);  new base G (2/5);  total (37/52)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGC  >  NC_000913/1101408‑1101704
                                                                                                                                                    |                                                                                                                                                    
aCCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATTCCGAAACAAAGATTTTACCCGTCGGCGCt                                                                                                                                                      >  2:8437/1‑149 (MQ=255)
aCCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCt                                                                                                                                                      >  2:322013/1‑149 (MQ=255)
aCCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCt                                                                                                                                                      >  2:332031/1‑149 (MQ=255)
 ccAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTg                                                                                                                                                     <  2:657385/149‑1 (MQ=255)
  cAGCGAGAATCTTTCAGTGACGTGACCAGCAATCCCGTGGAGCTTTGCGGAACAGCAGTGGAGAAGTTACTAGCCGGGTAGGGTTAAAATAGCCCGCTTTCGTCCTGAATGTTGTATACCGAAACAAATATATTACCCGTCGGCGCTgg                                                                                                                                                    >  1:1170598/1‑149 (MQ=255)
      gagaATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGa                                                                                                                                                <  1:1001644/149‑1 (MQ=255)
         aaTCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATgg                                                                                                                                             <  1:1058333/149‑1 (MQ=255)
           tCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTATATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACACGTCGGCGCTGGCAGATGGGt                                                                                                                                           <  2:1142691/149‑1 (MQ=255)
              ttCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGTCaaa                                                                                                                                       <  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:238807/1‑149 (MQ=255)
                            cAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGctct                                                                                                                          >  2:825497/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:672084/149‑1 (MQ=255)
                                                                                tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      <  1:573542/149‑1 (MQ=255)
                                                                                tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      <  1:415099/149‑1 (MQ=255)
                                                                                tAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      <  1:1519687/149‑1 (MQ=255)
                                                                                tAGGGGTTAAATTGCCCGGTGTCGTCCTGACTGTTGTATACCGAAACAAAGATTTTAACCGTCGGCGCTGGCAGATGTGTCAAATCTTTGTAGCTCTGAGCAAGAGTCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      <  1:262551/149‑1 (MQ=255)
                                                                                    gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                  <  2:888215/149‑1 (MQ=255)
                                                                                    gTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCAt                                                                  >  2:293570/1‑149 (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:791461/149‑1 (MQ=255)
                                                                                                           tGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAAc                                           <  2:796738/149‑1 (MQ=255)
                                                                                                              aTGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCCCGAGGCAGTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGCCGGCAACCaa                                        >  2:190841/1‑149 (MQ=255)
                                                                                                                    tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      >  1:138511/1‑113 (MQ=255)
                                                                                                                    tataCCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTaa                                                                      <  2:138511/113‑1 (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              <  1:825497/149‑1 (MQ=255)
                                                                                                                                        aCCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTAttt              <  2:147391/149‑1 (MQ=255)
                                                                                                                                         accttCGGCGCTTGCAGATGGGTAAAATCTTTGTAGCTCTGAGCACGATGCATTAATGTCGGTATGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCATCCATAAGAAATAAGCGCTGCATGATTAtttt             <  2:1390330/145‑1 (MQ=255)
                                                                                                                                          ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc            <  2:38479/149‑1 (MQ=255)
                                                                                                                                          ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc            <  2:1440911/149‑1 (MQ=255)
                                                                                                                                          ccGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTc            <  2:923131/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)
                                                                                                                                                    |                                                                                                                                                    
ACCAGCGAGAATCTTTCAGTGCCGTGACCAGCATTGCCGTGGCGCTTTGCGGAACAGCAGTGGAGAAGTTACTTGCCGGGTAGGGTTTAAATTGCCCGGTTTCGTCCTGAATGTTGTATACCGAAACAAAGATTTTACCCGTCGGCGCTGGCAGATGGGTCAAATCTTTGTAGCTCTGAGCACGAGGCATTAATGTCGGTCTGGCGGCTTCTTTAGGCGGGGCGGTTAAGCATCCGCTCAGTAACATGACGGCAACCAAAAGAAATAAGCGCTGCATGATTATTTTCCTTATGAAGC  >  NC_000913/1101408‑1101704

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: