Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,807,091 T→C 18.3% intergenic (+17/‑41) proW → / → proX glycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,807,0910TC18.3% 202.9 / 42.3 104intergenic (+17/‑41)proW/proXglycine betaine/proline ABC transporter permease/glycine betaine/proline ABC transporter periplasmic binding protein; cold shock protein
Reads supporting (aligned to +/- strand):  ref base T (42/43);  new base C (2/17);  total (44/60)
Fisher's exact test for biased strand distribution p-value = 1.79e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

Only 100 of 105 total aligned reads displayed.

GCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTT  >  NC_000913/2806993‑2807190
                                                                                                  |                                                                                                   
gcgcGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTg                                                                                                    <  2:2373368/100‑1 (MQ=255)
  gcgACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGAc                                                                                                  <  2:1412428/100‑1 (MQ=255)
     aCTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGcc                                                                                               >  2:2447253/1‑100 (MQ=255)
       tCACGCAGTCGCGGCAACCGTCGCTGGTAGACTGATGGCCCTGTTGGTCTGCTGACCCGGCCATTCATTAAGTCCCTCTGCACTTGCCCGGTGACGc                                                                                                >  2:828718/1‑97 (MQ=255)
        cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCggg                                                                                            >  2:1844868/1‑100 (MQ=255)
        cACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTAACCCGCCCATCCATTAAGTACCTCTGCACTTGCCCGGTGACGCCggg                                                                                            <  2:1584030/100‑1 (MQ=255)
          cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                          <  1:797899/100‑1 (MQ=255)
          cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                          <  2:1283271/100‑1 (MQ=255)
          cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                          <  2:1777856/100‑1 (MQ=255)
          cGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCa                                                                                          <  1:2211292/100‑1 (MQ=255)
           gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                         >  2:1284637/1‑100 (MQ=255)
           gCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCAt                                                                                         >  1:2106587/1‑100 (MQ=255)
            ctgTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTGGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  2:1007526/98‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  2:806933/100‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  1:459905/100‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  1:410224/100‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  1:1368457/100‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  1:2549721/100‑3 (MQ=255)
            cAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaag                                                                                        <  1:1129956/100‑3 (MQ=255)
             aGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTa                                                                                       <  1:614889/100‑1 (MQ=255)
              gTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtg                                                                                      <  2:2086471/100‑5 (MQ=255)
                cgcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCa                                                                                    <  1:552051/100‑1 (MQ=255)
                cgcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtgca                                                                                    <  1:2456679/100‑7 (MQ=255)
                 gcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAc                                                                                   <  2:603845/100‑1 (MQ=255)
                 gcgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAc                                                                                   <  1:1997023/100‑1 (MQ=255)
                  cgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAcc                                                                                  <  1:2205165/100‑1 (MQ=255)
                  cgGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGACGCCGGGCATTATCAcc                                                                                  >  1:2411450/1‑100 (MQ=255)
                   ggCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGTCCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtgcagag                                                                                 <  1:113216/100‑10 (MQ=37)
                   ggCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtgcagag                                                                                 <  1:790049/100‑10 (MQ=38)
                   ggCAACCGTCGCTGGTACACCACTGGCCCTGTTGGCCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGCGTCACCGGGCaagtgcagag                                                                                 <  1:2276296/100‑10 (MQ=37)
                      aaCCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGc                                                                              <  1:1809105/100‑1 (MQ=255)
                         cGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaa                                                                           >  1:1338213/1‑100 (MQ=255)
                             gCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCaaaaaa                                                                       >  2:2203341/1‑100 (MQ=255)
                               tagTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGTCATTATCACCAAAACAAAAAAAg                                                                     <  2:1993006/98‑1 (MQ=255)
                               tGGTACACCACTGGCCCTGTTTGTCTGCTGACCCGCCCATTCATTACGTAACTCTCCACTTGCCCGTTGACGCCGGGCATTATCACCCTGCCaaa                                                                          <  2:1057083/95‑1 (MQ=255)
                               tGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg                                                                     <  1:616678/100‑1 (MQ=255)
                               tGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAg                                                                     <  2:1078120/100‑1 (MQ=255)
                                ngTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                    >  2:1948361/2‑100 (MQ=255)
                                ggTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAgg                                                                    <  1:1090288/100‑1 (MQ=255)
                                    caccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATa                                                                >  1:2529817/1‑100 (MQ=255)
                                    caccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATa                                                                >  1:2221509/1‑100 (MQ=255)
                                     accacTGGCCCTGTTGGTCTGCTGACCCGACCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATaa                                                               <  1:2163110/100‑1 (MQ=255)
                                      ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAAc                                                              >  2:1805676/1‑100 (MQ=255)
                                      ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAAc                                                              <  1:1380786/100‑1 (MQ=255)
                                      ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAAc                                                              >  1:806919/1‑100 (MQ=255)
                                      ccacTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAAc                                                              >  2:504674/1‑100 (MQ=255)
                                           ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGc                                                         >  2:1716563/1‑100 (MQ=255)
                                           ggCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATAACCCTAAAAAAAAAAGGAATAACAATGc                                                         <  2:806181/100‑1 (MQ=255)
                                              ccTGTTGGTCTGGTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCAAAATAACAAAAAAAGGAATAACAATGCGAc                                                      <  2:2544356/100‑1 (MQ=255)
                                              ccTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGAc                                                      <  1:2057513/100‑1 (MQ=255)
                                                tGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACAt                                                    <  1:907604/100‑1 (MQ=255)
                                                   tGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGc                                                 >  2:859690/1‑100 (MQ=255)
                                                    ggTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTCTCACCCAGCCAAAAAAAGGAATAACAATGCGACATAGCg                                                >  2:1801707/1‑100 (MQ=255)
                                                     gTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGt                                               <  1:737841/100‑1 (MQ=255)
                                                      tCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACGTAGCGTa                                              <  2:1706026/100‑1 (MQ=255)
                                                       ctgctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTAc                                             >  2:1952851/1‑100 (MQ=255)
                                                       ctgctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTAc                                             <  1:614429/100‑1 (MQ=255)
                                                        tgctgACCCGCCCATTCATTAAGTAACTCTGCACTTTCCCGGTGAAGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTAc                                             >  2:513357/1‑99 (MQ=255)
                                                        tgctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACt                                            >  2:1763654/1‑100 (MQ=255)
                                                        tgctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCCGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACt                                            <  1:1200003/100‑1 (MQ=255)
                                                         gctgACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACtt                                           >  1:1889211/1‑100 (MQ=255)
                                                            gACCCGCCCATTCATTAAGTAACTCTGCACTAGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACttttt                                        >  2:2465033/1‑100 (MQ=255)
                                                              cccGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGc                                      >  2:1018302/1‑100 (MQ=255)
                                                                cGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCg                                     <  2:1133230/99‑1 (MQ=255)
                                                                 gCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAATAGGAATAACAATGCGACATAGCGTACTTTTTGCGAc                                   <  1:1801721/100‑1 (MQ=255)
                                                                 gCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGAc                                   >  1:1031723/1‑100 (MQ=255)
                                                                 gCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGAc                                   <  1:2054458/100‑1 (MQ=255)
                                                                 gCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGAc                                   <  2:498036/100‑1 (MQ=255)
                                                                   ccATTCATTAATTAACTCTGCACTTGCTCGGTGACAACAGGCATTATGACCATGCCAAAAATATGAATAACATTACGATAATGCGTAATTTTTGCAACa                                  >  2:483604/1‑99 (MQ=255)
                                                                      ttcattAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGt                              >  1:489557/1‑100 (MQ=255)
                                                                           tAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGcc                         <  1:2328935/100‑1 (MQ=255)
                                                                            aaGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCa                        >  1:1763189/1‑100 (MQ=255)
                                                                               tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGc                     <  1:894724/100‑1 (MQ=255)
                                                                               tAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGACACGc                     >  2:1020749/1‑100 (MQ=255)
                                                                                aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                    >  2:730113/1‑100 (MQ=255)
                                                                                aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                    <  1:2303592/100‑1 (MQ=255)
                                                                                aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                    <  1:1530955/100‑1 (MQ=255)
                                                                                aaCTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                    <  2:483477/100‑1 (MQ=255)
                                                                                aaCTCTGCACTTACCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCt                    >  2:2533901/1‑100 (MQ=255)
                                                                                  ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                  <  1:1940689/100‑1 (MQ=255)
                                                                                  ctctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                  >  1:2092206/1‑100 (MQ=255)
                                                                                   tctGCACTTGCCGGGTGACGCAGGGCATTATTACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                  <  2:2388798/99‑1 (MQ=255)
                                                                                   tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAt                 >  1:515837/1‑100 (MQ=255)
                                                                                   tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAt                 >  1:2153612/1‑100 (MQ=255)
                                                                                   tctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAt                 >  1:112508/1‑100 (MQ=255)
                                                                                   tctGCACTTGCCCGGTGACGCCGGGCAGTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTa                  <  2:2411436/99‑1 (MQ=255)
                                                                                    ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtc                >  1:2167736/1‑100 (MQ=255)
                                                                                    ctGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtc                <  2:650148/100‑1 (MQ=255)
                                                                                     tgataatGCCCGGCGTCACCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtct               >  2:609487/7‑100 (MQ=255)
                                                                                      gCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTAtctc              <  1:2186273/100‑1 (MQ=255)
                                                                                        aCTTGCCCGGTGGCGCCGGGCATTATCACCCTACCAGAAAAAGGAATAACCATGCGACACAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTCTCTCTa            >  2:2491170/1‑100 (MQ=255)
                                                                                          ttGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTaca          <  2:489552/100‑1 (MQ=255)
                                                                                           tGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTacac         >  2:2286257/1‑100 (MQ=255)
                                                                                            gCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGTGTTAGTGTTTGCCACGCTTATCTCTAcaca        >  2:157567/1‑100 (MQ=255)
                                                                                            gCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTAcaca        <  1:805570/100‑1 (MQ=255)
                                                                                             cccGGTGACGCCGAGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaa       >  2:1210500/1‑100 (MQ=255)
                                                                                              ncGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACaaa      >  1:72073/2‑100 (MQ=255)
                                                                                                ggTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACt    >  2:1516439/1‑100 (MQ=255)
                                                                                                ggTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCCACAGCGTTTGCCACGCTTATCTCTACACAAACt    <  2:1692453/100‑1 (MQ=255)
                                                                                                  tGACGCCGGGCATTATCACCCTGCCAAAAAAAGTAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACttt  <  1:297093/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
GCGCGACTCACGCAGTCGCGGCAACCGTCGCTGGTACACCACTGGCCCTGTTGGTCTGCTGACCCGCCCATTCATTAAGTAACTCTGCACTTGCCCGGTGACGCCGGGCATTATCACCCTGCCAAAAAAAGGAATAACAATGCGACATAGCGTACTTTTTGCGACAGCGTTTGCCACGCTTATCTCTACACAAACTTT  >  NC_000913/2806993‑2807190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: