Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,100,793 A→C 8.2% intergenic (+6/+58) 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,7930AC 8.2% 516.6 / 36.1 159intergenic (+6/+58)ycdZ/csgGDUF1097 family inner membrane protein/curli production assembly/transport outer membrane lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (73/73);  new base C (8/5);  total (81/78)
Fisher's exact test for biased strand distribution p-value = 5.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.18e-01

Only 100 of 162 total aligned reads displayed.

GGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCA  >  NC_000913/1100645‑1100938
                                                                                                                                                    |                                                                                                                                                 
 gCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCAc                                                                                                                                                  >  1:164734/1‑149 (MQ=255)
   aGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTc                                                                                                                                                >  4:1646691/1‑149 (MQ=255)
       tGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGc                                                                                                                                            <  3:274107/149‑1 (MQ=255)
       tGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGc                                                                                                                                            <  2:1658158/149‑1 (MQ=255)
           tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                        <  1:366630/149‑1 (MQ=255)
           tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTc                                                                                                                                        >  2:756165/1‑149 (MQ=255)
           tGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGAGAGCAGGAAATCAAAAATAAAGCGTGAGTGCCCCTCACGCttta                                                                                                                                        <  4:811954/149‑2 (MQ=255)
             gAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGc                                                                                                                                      >  1:323579/1‑149 (MQ=255)
                  tGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                 >  4:1385774/1‑149 (MQ=255)
                    gtgtTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACa                                                                                                                               >  2:126606/1‑149 (MQ=255)
                     tgtTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAg                                                                                                                              >  1:1554087/1‑149 (MQ=255)
                       ttACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                 >  1:956215/1‑144 (MQ=255)
                         aCCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGTCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaa                                                                                                                          >  2:522409/1‑149 (MQ=255)
                           cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCAc                                                                                                                                                  >  2:721107/1‑123 (MQ=255)
                           cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCAc                                                                                                                                                  <  1:721107/123‑1 (MQ=255)
                           cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                        >  4:1643544/1‑149 (MQ=255)
                           cTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                        >  1:888387/1‑149 (MQ=255)
                            ttCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATg                                                                                                                       >  4:132104/1‑149 (MQ=255)
                                 tGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                  <  2:1445842/149‑1 (MQ=255)
                                 tGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                  <  1:1071746/149‑1 (MQ=255)
                                 tGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGga                                                                                                                  <  1:768231/149‑1 (MQ=255)
                                  ggCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgat                                                                                                                 <  3:561190/149‑1 (MQ=255)
                                    cgcTTGGGCGGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGGGCAAGGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatga                                                                                                               <  1:480859/149‑1 (MQ=255)
                                     gcTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                        <  3:1356937/139‑1 (MQ=255)
                                     gcTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAAt                                                                                                                        >  4:1356937/1‑139 (MQ=255)
                                     gcTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGgatgat                                                                                                              >  2:646048/1‑149 (MQ=255)
                                           gCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                        <  3:1469111/149‑1 (MQ=255)
                                           gCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTcc                                                                                                        >  1:908108/1‑149 (MQ=255)
                                            cTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCg                                                                                                       <  2:1378478/149‑1 (MQ=255)
                                             tGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCgg                                                                                                      >  4:743881/1‑149 (MQ=255)
                                              gATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGc                                                                                                     <  2:294298/149‑1 (MQ=255)
                                                tatTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                 >  4:92277/1‑119 (MQ=255)
                                                tatTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTaaa                                                                                                                                 <  3:92277/119‑1 (MQ=255)
                                                     ggTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCt                                                                                                                                     <  2:1301329/110‑1 (MQ=255)
                                                          cGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGt                                                                                         <  2:996410/149‑1 (MQ=255)
                                                          cGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGt                                                                                         <  1:391322/149‑1 (MQ=255)
                                                           gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGGAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                        >  2:251521/1‑149 (MQ=255)
                                                           gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCa                                                                                                                                                   >  3:490266/1‑90 (MQ=255)
                                                           gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCa                                                                                                                                                   <  4:490266/90‑1 (MQ=255)
                                                           gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                        >  2:65509/1‑149 (MQ=255)
                                                           gCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTc                                                                                        >  1:813762/1‑149 (MQ=255)
                                                            cAATGACAAACAGTGGCCAGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAGAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACTGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCa                                                                                       <  4:813473/149‑1 (MQ=255)
                                                             aaTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCg                                                                                                       >  2:1167538/1‑132 (MQ=255)
                                                             aaTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAg                                                                                      >  1:879362/1‑149 (MQ=255)
                                                              aTGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAgg                                                                                     >  2:403020/1‑149 (MQ=255)
                                                                      cAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             <  4:276239/149‑1 (MQ=255)
                                                                      cAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             <  2:323579/149‑1 (MQ=255)
                                                                      cAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             <  4:60222/149‑1 (MQ=255)
                                                                      cAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             <  2:1535973/149‑1 (MQ=255)
                                                                           gCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAc                                                                        >  1:1227440/1‑149 (MQ=255)
                                                                            ccTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAcc                                                                       <  1:1289870/149‑1 (MQ=255)
                                                                                  gctggcGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCGTTTTTATCTGTCAGGATTCCGGTGGAACCGACata                                                                 >  4:1098809/1‑149 (MQ=255)
                                                                                       cggcgCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCg                                                            <  2:1006914/149‑1 (MQ=255)
                                                                                           gcgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                        >  3:469745/1‑149 (MQ=255)
                                                                                            cgTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTaaaa                                                                                                                         >  2:613339/1‑83 (MQ=255)
                                                                                                gTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCAc                                                   <  3:1525759/149‑1 (MQ=255)
                                                                                                  gCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCa                                                 <  1:1345824/149‑1 (MQ=255)
                                                                                                  gCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCa                                                 >  2:524176/1‑149 (MQ=255)
                                                                                                   cAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAg                                                >  1:761730/1‑149 (MQ=255)
                                                                                                      aGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTc                                                                                                         <  1:1494590/89‑1 (MQ=255)
                                                                                                       gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                            >  3:489904/1‑149 (MQ=255)
                                                                                                       gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                            >  4:1091260/1‑149 (MQ=255)
                                                                                                       gACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATg                                            >  1:709836/1‑149 (MQ=255)
                                                                                                        aCGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGt                                           >  1:1321747/1‑149 (MQ=255)
                                                                                                          ggCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCa                                         <  1:646048/149‑1 (MQ=255)
                                                                                                            cgcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCAtt                                       >  1:1029879/1‑149 (MQ=255)
                                                                                                            cgcACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCAtt                                       <  2:709836/149‑1 (MQ=255)
                                                                                                                     aGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTc                              <  2:697730/149‑1 (MQ=255)
                                                                                                                        aGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGc                           <  1:499005/149‑1 (MQ=255)
                                                                                                                         ggAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCt                          >  4:1172245/1‑149 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                        >  3:1617672/1‑112 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTAt                                                        <  4:1617672/112‑1 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                   >  4:33255/1‑149 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                   <  2:879362/149‑1 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                   >  1:1673292/1‑149 (MQ=255)
                                                                                                                                aaaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTAtttt                   >  1:1281063/1‑149 (MQ=255)
                                                                                                                                 aaaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTg                  <  2:501810/149‑1 (MQ=255)
                                                                                                                                  aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACatat                                                                <  1:462246/102‑1 (MQ=255)
                                                                                                                                  aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACatat                                                                >  2:462246/1‑102 (MQ=255)
                                                                                                                                  aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                             >  3:1687385/1‑105 (MQ=255)
                                                                                                                                  aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGc                                                             <  4:1687385/105‑1 (MQ=255)
                                                                                                                                  aaaTAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGc                 <  2:1321747/149‑1 (MQ=255)
                                                                                                                                       aaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTg                  <  3:1175832/143‑1 (MQ=255)
                                                                                                                                       aaGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATc            >  2:855108/1‑149 (MQ=255)
                                                                                                                                       aaGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATc            <  3:821160/149‑1 (MQ=255)
                                                                                                                                        aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAc                                                                        <  4:1258344/88‑1 (MQ=39)
                                                                                                                                        aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAAc                                                                        >  3:1258344/1‑88 (MQ=39)
                                                                                                                                        aGCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATcc           >  3:1399653/1‑149 (MQ=255)
                                                                                                                                        aGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             <  1:1402021/83‑1 (MQ=255)
                                                                                                                                        aGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGt                                                                             >  2:1402021/1‑83 (MQ=255)
                                                                                                                                         gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc          >  2:481626/1‑149 (MQ=255)
                                                                                                                                         gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc          >  1:140887/1‑149 (MQ=255)
                                                                                                                                         gCGTGAGTGCCCCTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc          >  3:1018034/1‑149 (MQ=255)
                                                                                                                                         gCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATccc          <  2:900463/149‑1 (MQ=255)
                                                                                                                                          cGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCca         <  2:1632044/149‑1 (MQ=255)
                                                                                                                                            tGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCATAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca       <  3:1098809/149‑1 (MQ=255)
                                                                                                                                                gtgCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGGTTCCGGTGGAACCGACCTATGGCGGTATTTCACCAGGATGTCATTCTGCCGTTCTGCTTTATTTTGCGAATCCCACAGAcc   >  2:612105/3‑149 (MQ=255)
                                                                                                                                                gggCACTCAGGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca       <  2:1539310/145‑1 (MQ=255)
                                                                                                                                                gggCACTCAGGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCcaca       >  1:1539310/1‑145 (MQ=255)
                                                                                                                                                 ggCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCa  <  2:113200/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGCAGGGTGACTGGAAACTGGTGTTACCTTCGCTGGCGCTTGGGCTGATATTTGGTTACGCAATGAAAAACAGTGGCCTGTGGCTGGCGGCGCGTAGTGCAAAGACGGCGCACCGTGAGCAGGAAATCAAAAATAAAGCGTGAGGGGCACTCACGCTTTCGCTTAAACAGTAAAATGCCGGATGATAATTCCGGCTTTTTTATCTGTCAGGATTCCGGTGGAACCGACATATGGCGGTATTTCACCAGAATGTCATTCTGCCGTTCTGCTTTATTTTGCAAATCCCACAGACCA  >  NC_000913/1100645‑1100938

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