Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,789 G→T 12.8% intergenic (+2/+62) ycdZ → / ← csgG DUF1097 family inner membrane protein/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,100,7890GT12.8% 252.1 / 22.2 78intergenic (+2/+62)ycdZ/csgGDUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base G (33/35);  new base T (6/4);  total (39/39)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGA  >  NC_000913/1100641‑1100935
                                                                                                                                                    |                                                                                                                                                  
gCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAgg                                                                                                                                                    <  2:717434/149‑1 (MQ=255)
 ccGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAggg                                                                                                                                                   <  3:175163/149‑1 (MQ=255)
 ccGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAggg                                                                                                                                                   <  4:115791/149‑1 (MQ=255)
        gggTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCa                                                                                                                                            <  3:451789/149‑1 (MQ=255)
         ggTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCAc                                                                                                                                           <  3:477080/149‑1 (MQ=255)
          gTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACg                                                                                                                                          >  4:5238/1‑149 (MQ=255)
           tGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCAGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGc                                                                                                                                         <  2:83022/149‑1 (MQ=255)
            gACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCt                                                                                                                                        >  2:473693/1‑149 (MQ=255)
             aCTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCtt                                                                                                                                       <  2:144078/149‑1 (MQ=255)
               tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                     >  3:80549/1‑149 (MQ=255)
               tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                     <  1:286779/149‑1 (MQ=255)
                ggAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCg                                                                                                                                    >  2:808123/1‑149 (MQ=255)
                ggAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCg                                                                                                                                    >  4:812124/1‑149 (MQ=255)
                 gAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCt                                                                                                                                        >  4:460501/1‑144 (MQ=255)
                 gAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCt                                                                                                                                        <  3:460501/144‑1 (MQ=255)
                          gtTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGt                                                                                                                          >  3:33117/1‑149 (MQ=255)
                           ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATATAGCGTGAGGGGCACTCACGCTTTCGc                                                                                                                                   <  4:561818/139‑1 (MQ=255)
                           ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGc                                                                                                                                   >  3:561818/1‑139 (MQ=255)
                           ttACCTTCGCTGGCCCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTa                                                                                                                         >  3:63361/1‑149 (MQ=255)
                             aCCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaa                                                                                                                       >  3:545274/1‑149 (MQ=255)
                              ccTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaaa                                                                                                                      <  2:378031/149‑1 (MQ=255)
                               cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                     >  3:328882/1‑149 (MQ=255)
                                ttCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATg                                                                                                                    >  3:361703/1‑149 (MQ=255)
                                  cGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGGCTGTGGCTGGCGGCGAGTAGTGCAAAGGCTGCGCCCAGTGAGCAGGAAAACAAAAATAACGCGTGAGGGGCACTCCCGCTTTCGCTTAAACAGTAAAATGcc                                                                                                                  <  2:291014/149‑1 (MQ=255)
                                  cGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGcc                                                                                                                  <  1:244013/149‑1 (MQ=255)
                                   gCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCg                                                                                                                 <  4:863741/149‑1 (MQ=255)
                                         gcgtGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGGAATCAAAAATAAAGCGTGCGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatgat                                                                                                           >  3:202456/4‑149 (MQ=255)
                                           ttGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGCAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCttt                                                                                                                                      <  3:35694/120‑1 (MQ=255)
                                           ttGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGCAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCttt                                                                                                                                      >  4:35694/1‑120 (MQ=255)
                                              ggCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTc                                                                                                      <  4:604758/149‑1 (MQ=255)
                                               gCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                     >  3:316232/1‑149 (MQ=255)
                                               gCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                     >  3:288909/1‑149 (MQ=255)
                                                     atTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttt                                                                                               <  2:552496/149‑1 (MQ=255)
                                                      tttGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCtttt                                                                                              >  1:380522/1‑149 (MQ=255)
                                                      tttGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCtttt                                                                                              >  4:222343/1‑149 (MQ=255)
                                                       ttGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCttttt                                                                                             >  2:829723/1‑149 (MQ=255)
                                                         ggTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATCAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTa                                                                                           <  4:63361/149‑1 (MQ=255)
                                                          gTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTa                                                                                                                         >  4:695638/1‑118 (MQ=255)
                                                          gTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTa                                                                                                                         <  3:695638/118‑1 (MQ=255)
                                                           ttACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAATACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACTGTAAAATGCCGGATCATACTTCCGGCTTTTTTATc                                                                                         <  2:195947/149‑1 (MQ=255)
                                                              cGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGt                                                                                      <  1:839454/149‑1 (MQ=255)
                                                               gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                     >  4:486642/1‑149 (MQ=255)
                                                                cAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                    <  4:357485/149‑1 (MQ=255)
                                                                cAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                    <  4:545274/149‑1 (MQ=255)
                                                                cAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                    <  3:812124/149‑1 (MQ=255)
                                                                cAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                    <  1:43204/149‑1 (MQ=255)
                                                                cAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTCCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                    <  2:436389/149‑1 (MQ=255)
                                                                 aaTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGt                                                                                                                          <  3:866967/110‑1 (MQ=255)
                                                                 aaTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGt                                                                                                                          >  4:866967/1‑110 (MQ=255)
                                                                 aaTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAg                                                                                   >  3:707160/1‑149 (MQ=255)
                                                                    gAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGAt                                                                                <  1:863510/149‑1 (MQ=255)
                                                                     aaaaaCAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGAtt                                                                               >  2:362084/1‑149 (MQ=255)
                                                                         aCAGTGGCCTGTGGATGGCGGCGCGTAGTGCAAAGACGGCGCCCCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCgg                                                                           <  2:550105/149‑1 (MQ=255)
                                                                          cAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTTAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                          <  1:577481/149‑1 (MQ=255)
                                                                           aGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAGACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTg                                                                         <  1:808123/149‑1 (MQ=255)
                                                                                 cTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCg                                                                   >  2:454102/1‑149 (MQ=255)
                                                                                   gtgGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGAc                                                                 <  4:328882/149‑1 (MQ=255)
                                                                                        tggcGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATg                                                            >  2:867550/1‑149 (MQ=255)
                                                                                         ggcggcGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTACGGTGGAACCGACATATgg                                                           <  1:766633/149‑1 (MQ=255)
                                                                                          gcggcgCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                          >  4:5595/1‑149 (MQ=255)
                                                                                           cggcgCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCg                                                         <  4:769047/149‑1 (MQ=255)
                                                                                           cggcgCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCg                                                         <  2:313028/149‑1 (MQ=255)
                                                                                                   aGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGCTTCCGGTGGAACCGACATATGGCGGTATTTCa                                                 >  1:388550/1‑149 (MQ=255)
                                                                                                     tGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCAcc                                               <  1:627681/149‑1 (MQ=255)
                                                                                                         aaGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGaa                                           >  1:538887/1‑149 (MQ=255)
                                                                                                            aCGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCAGGTGGACCCGACATATGGCGGTATTTCACCAGAATGt                                        >  2:227295/1‑149 (MQ=255)
                                                                                                                       tGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGt                             >  4:809181/1‑149 (MQ=255)
                                                                                                                                  tCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtt                  >  1:856195/1‑149 (MQ=255)
                                                                                                                                   cAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAttt                 <  2:270341/149‑1 (MQ=255)
                                                                                                                                    aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTc                                       <  4:455263/126‑1 (MQ=255)
                                                                                                                                    aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTc                                       >  3:455263/1‑126 (MQ=255)
                                                                                                                                    aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                <  3:5595/149‑1 (MQ=255)
                                                                                                                                     aaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTg               >  1:5231/1‑149 (MQ=255)
                                                                                                                                       gcgaaaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCa             >  2:10911/4‑149 (MQ=255)
                                                                                                                                           aaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATc         >  2:383504/1‑149 (MQ=255)
                                                                                                                                            aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGGAAATcc        >  2:657195/1‑149 (MQ=255)
                                                                                                                                            aGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATcc        <  3:870682/149‑1 (MQ=255)
                                                                                                                                             gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc       >  3:69227/1‑149 (MQ=255)
                                                                                                                                             gCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc       <  1:473693/149‑1 (MQ=255)
                                                                                                                                               gTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCTGTGGAACCGACATATGGCGGTATTTAACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcac     >  2:550115/1‑149 (MQ=255)
                                                                                                                                                  aGGGGCACTCACGCTTTAGCTTAAACAGTCAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGa  <  4:515244/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCCGGGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGA  >  NC_000913/1100641‑1100935

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