Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,376,262 C→A 8.0% intergenic (+26/+17) dcuD → / ← sspB putative transporter DcuD/ClpXP protease specificity‑enhancing factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,376,2620CA 8.0% 371.9 / 13.2 112intergenic (+26/+17)dcuD/sspBputative transporter DcuD/ClpXP protease specificity‑enhancing factor
Reads supporting (aligned to +/- strand):  ref base C (59/44);  new base A (4/5);  total (63/49)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

Only 100 of 113 total aligned reads displayed.

GGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGC  >  NC_000913/3376116‑3376407
                                                                                                                                                  |                                                                                                                                                 
ggtggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTg                                                                                                                                                >  1:1872981/1‑150 (MQ=255)
  tggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTggg                                                                                                                                              <  2:836057/150‑1 (MQ=255)
  gggtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTggg                                                                                                                                              <  1:234132/149‑1 (MQ=255)
    gtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGcc                                                                                                                                            >  1:1370997/1‑150 (MQ=255)
    gtCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGcc                                                                                                                                            <  1:1982033/150‑1 (MQ=255)
       gTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGt                                                                                                                                         <  2:453963/150‑1 (MQ=255)
         tACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGttt                                                                                                                                       >  1:1736272/1‑150 (MQ=255)
         tACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGttt                                                                                                                                       <  2:835322/150‑1 (MQ=255)
           cGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTg                                                                                                                                     >  1:490376/1‑150 (MQ=255)
           cGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTg                                                                                                                                     <  1:654433/150‑1 (MQ=255)
           cGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTg                                                                                                                                     <  1:842812/150‑1 (MQ=255)
              cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTg                                                                                                                                     >  2:1898155/1‑147 (MQ=255)
              cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTg                                                                                                                                     <  1:1898155/147‑1 (MQ=255)
              cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTGtct                                                                                                                                  <  1:1365569/150‑3 (MQ=255)
              cAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTGtct                                                                                                                                  <  1:162014/150‑3 (MQ=255)
               aGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTAtt                                                                                                                                 <  2:265641/150‑1 (MQ=255)
               aGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCCAGGCGGCCTGTTTTGTCtt                                                                                                                                 <  2:1474355/150‑1 (MQ=255)
                 gAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTAc                                                                                                                               <  1:789306/150‑1 (MQ=255)
                     tGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTca                                                                                                                           <  1:1366331/150‑1 (MQ=255)
                         ggCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAc                                                                                                                       >  2:1113403/1‑150 (MQ=255)
                            gTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAcgcg                                                                                                                    >  2:156066/1‑150 (MQ=255)
                            gTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAcgcg                                                                                                                    <  1:322851/150‑1 (MQ=255)
                            gTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAAcgcg                                                                                                                    <  1:587904/150‑1 (MQ=255)
                             tttCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGt                                                                                                                   <  1:1977258/150‑1 (MQ=255)
                             tttCGCCTTTTGCGGTGGGGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGt                                                                                                                   <  1:1863507/150‑1 (MQ=255)
                                cGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAAt                                                                                                                >  1:1393213/1‑150 (MQ=255)
                                    ttttGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCg                                                                                                            >  1:1505042/1‑150 (MQ=255)
                                       tGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTcaca                                                                                                                         <  1:576566/134‑1 (MQ=255)
                                       tGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTcaca                                                                                                                         >  2:576566/1‑134 (MQ=255)
                                        gCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCg                                                                                                        >  1:1682482/1‑150 (MQ=255)
                                         cGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGa                                                                                                       >  1:1111543/1‑150 (MQ=255)
                                           gtggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCCCACCGCGTAATGCCGGTCGacc                                                                                                     >  2:1056314/1‑150 (MQ=255)
                                             ggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGt                                                                                                                   >  2:1053740/1‑134 (MQ=255)
                                             ggtgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGt                                                                                                                   <  1:1053740/134‑1 (MQ=255)
                                               tgAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGaccaccc                                                                                                 >  1:483242/1‑149 (MQ=255)
                                                gAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACcgc                                                                                                >  2:1143387/1‑150 (MQ=255)
                                                gAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACcgc                                                                                                <  2:517703/150‑1 (MQ=255)
                                                gAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATAGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACcgc                                                                                                >  1:1203427/1‑150 (MQ=255)
                                                   gCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCATGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCgtg                                                                                             >  2:826451/1‑150 (MQ=255)
                                                       aCAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGtggtgg                                                                                         >  2:1558899/1‑150 (MQ=255)
                                                          gCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTg                                                                                      >  2:1709596/1‑150 (MQ=255)
                                                          gCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTg                                                                                      >  1:162307/1‑150 (MQ=255)
                                                          gCGATCCACATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTATTGCCGGTCGACCACCGCGTGGTGGCTg                                                                                      >  2:1369466/1‑150 (MQ=255)
                                                           cGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGc                                                                                     >  2:285502/1‑150 (MQ=255)
                                                           cGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGc                                                                                     >  2:958028/1‑150 (MQ=255)
                                                             aTCCCCATGGCAGGCGGTTTCGTGGCGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCgg                                                                                   <  2:1255100/150‑1 (MQ=255)
                                                              tCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGa                                                                                  <  2:27383/150‑1 (MQ=255)
                                                                cccATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAgg                                                                                >  1:224992/1‑150 (MQ=255)
                                                                  cATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGtt                                                                              >  2:27457/1‑150 (MQ=255)
                                                                   aTGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTc                                                                             <  2:558481/150‑1 (MQ=255)
                                                                    tGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCa                                                                            >  1:463400/1‑150 (MQ=255)
                                                                       cAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCg                                                                         <  2:235097/150‑1 (MQ=255)
                                                                          tCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCa                                                                      <  2:1682482/150‑1 (MQ=255)
                                                                           cGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAg                                                                     >  1:58940/1‑150 (MQ=255)
                                                                            ggTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAgg                                                                    <  2:483242/150‑1 (MQ=255)
                                                                               ttCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATg                                                                 <  2:1807022/150‑1 (MQ=255)
                                                                                 cGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAg                                                               <  2:1737040/150‑1 (MQ=255)
                                                                                 cGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAg                                                               <  1:1143387/150‑1 (MQ=255)
                                                                                  gtggtgAATATTATTGCCAAAAACACGCTATTTTATTAAGTCATTAAAAAGAAAAAACAGTCCGCCTGCGCCTGTTTTGTATTACTTAACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgt                                                              >  2:1000341/1‑150 (MQ=255)
                                                                                  gtggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgt                                                              >  1:287519/1‑150 (MQ=255)
                                                                                   tggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgtg                                                             >  1:1142863/1‑150 (MQ=255)
                                                                                   tggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAgtg                                                             <  1:1398364/150‑1 (MQ=255)
                                                                                    ggtgAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGtt                                                                              <  2:941042/132‑1 (MQ=255)
                                                                                            tGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCg                                                    <  2:1505042/150‑1 (MQ=255)
                                                                                            tGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCg                                                    >  1:1265302/1‑150 (MQ=255)
                                                                                              aTTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGACTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTg                                                  >  1:1149840/1‑150 (MQ=255)
                                                                                                 gCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATc                                               >  2:718070/1‑150 (MQ=255)
                                                                                                   cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTg                                             >  2:1269344/1‑150 (MQ=255)
                                                                                                   cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTg                                             <  2:224992/150‑1 (MQ=255)
                                                                                                   cacaATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTg                                             >  2:1310066/1‑150 (MQ=255)
                                                                                                     caATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGc                                           <  2:1736272/150‑1 (MQ=255)
                                                                                                     caATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGc                                           <  1:998018/150‑1 (MQ=255)
                                                                                                      aaTCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCt                                          >  2:1963514/1‑150 (MQ=255)
                                                                                                      aaTCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCt                                          <  1:1113403/150‑1 (MQ=255)
                                                                                                         cACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACAACCGCGAGGTGGCTCCGGAGGTTCATCGTCAGGATGAGTGTCATAATCGTGATCTGGCTTGt                                       >  2:1987924/1‑150 (MQ=255)
                                                                                                           cGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACAACCGCGTGGTGGCTGAGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCg                                     >  2:227537/1‑150 (MQ=255)
                                                                                                                ttttATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATc                                <  2:1393213/150‑1 (MQ=255)
                                                                                                                   tattaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAAt                             >  2:1372993/1‑150 (MQ=255)
                                                                                                                    attaAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATa                            >  2:1547750/1‑150 (MQ=255)
                                                                                                                      taAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAAc                          >  1:1767299/1‑150 (MQ=255)
                                                                                                                        aGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCg                        >  1:471792/1‑150 (MQ=255)
                                                                                                                          tCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGAc                      <  1:285502/150‑1 (MQ=255)
                                                                                                                           cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACa                     >  1:678914/1‑150 (MQ=255)
                                                                                                                           cATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACa                     >  1:1641630/1‑150 (MQ=255)
                                                                                                                             ttAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATa                   <  2:162307/150‑1 (MQ=255)
                                                                                                                             ttAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATa                   >  2:472187/1‑150 (MQ=255)
                                                                                                                              tAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATaa                  >  1:1776089/1‑150 (MQ=255)
                                                                                                                                gtaataCAAAACAGGCCCAGGCGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACg                >  1:1103074/6‑150 (MQ=255)
                                                                                                                                aaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACg                >  1:177326/1‑150 (MQ=255)
                                                                                                                                aaaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACg                >  1:670956/1‑150 (MQ=255)
                                                                                                                                 aaaGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACgg               <  2:463400/150‑1 (MQ=255)
                                                                                                                                    taCAAAACAGGCCCAGGCGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGttt            >  1:995770/2‑150 (MQ=255)
                                                                                                                                    gACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGttt            >  2:839166/1‑150 (MQ=255)
                                                                                                                                      cAAAACAGGCCGCCTGGGACTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGATGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg          >  2:1891732/1‑150 (MQ=255)
                                                                                                                                      cAAAACAGGCCCAGGCGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCg          >  2:1546965/1‑150 (MQ=255)
                                                                                                                                       aaaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGt         <  1:472187/150‑1 (MQ=255)
                                                                                                                                       aaaaCACGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGt         <  2:667762/150‑1 (MQ=255)
                                                                                                                                        aaaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGtt        >  2:1397545/1‑150 (MQ=255)
                                                                                                                                         aaCAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTg       <  2:1370997/150‑1 (MQ=255)
                                                                                                                                              gccgccTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGc  <  1:156066/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTGGTCGTTACGTCAGGAATTGCAGGCGTTTCGCCTTTTGCGGTGGTGAAGCGGACAGCGATCCCCATGGCAGTCGGTTTCGTGGTGAATATGATTGCCACAATCACGCTATTTTATTAAGTCATTAAAAAGACAAAACAGGCCGCCTGGGCCTGTTTTGTATTACTTCACAACGCGTAATGCCGGTCGACCACCGCGTGGTGGCTGCGGAGGTTCATCGTCAGGATGAGTGTCATCATCGTGATCTGGCTTGTCGCCATCAATAACCGACATAACGGTTTCGTTGTCTGC  >  NC_000913/3376116‑3376407

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