Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 972,900 T→C R47R (CGT→CGC)  smtA → S‑adenosylmethionine‑dependent methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913972,9000TC84.7% 643.6 / 102.2 281R47R (CGT→CGC) smtAS‑adenosylmethionine‑dependent methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (22/21);  new base C (101/137);  total (123/158)
Fisher's exact test for biased strand distribution p-value = 3.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 287 total aligned reads displayed.

TTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCA  >  NC_000913/972809‑972989
                                                                                           |                                                                                         
 ttACGGCACCACCATAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCgcgt                                                                                         >  2:6447404/1‑93 (MQ=255)
  tACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGtgtg                                                                                        >  2:7213719/1‑93 (MQ=255)
     ggCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTg                                                                                     <  1:7087232/93‑1 (MQ=255)
       caccacCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGa                                                                                   <  1:1311241/93‑1 (MQ=255)
         ccaccaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATg                                                                                 >  2:1782473/1‑93 (MQ=255)
          caccaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGc                                                                                <  2:4372359/93‑1 (MQ=255)
           accaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGGCCGCAAAAACTGCGCGTGCTGGATGCt                                                                               >  1:6287639/1‑93 (MQ=255)
           accaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCt                                                                               >  2:786536/1‑93 (MQ=255)
           accaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCt                                                                               >  1:4516837/1‑93 (MQ=255)
           accaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCt                                                                               <  2:5852494/93‑1 (MQ=255)
            ccaAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTg                                                                              <  2:2732936/93‑1 (MQ=255)
              aaaGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGc                                                                            >  1:3333765/1‑93 (MQ=255)
              aaaGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGc                                                                            <  2:400706/93‑1 (MQ=255)
              aaaGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGc                                                                            <  2:3819696/93‑1 (MQ=255)
              aaaGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGc                                                                            <  2:3365680/93‑1 (MQ=255)
               aaGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCg                                                                           >  1:1879747/1‑93 (MQ=255)
                 gggCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGt                                                                         >  1:7232164/1‑93 (MQ=255)
                 gggCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCggg                                                                         >  1:657035/1‑92 (MQ=255)
                  ggCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTg                                                                        >  2:2027016/1‑93 (MQ=255)
                  ggCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTg                                                                        <  2:522143/93‑1 (MQ=255)
                   gCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTgg                                                                       <  1:3492869/93‑1 (MQ=255)
                    cAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGga                                                                      <  2:6601156/93‑1 (MQ=255)
                    cAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGga                                                                      <  1:3443385/93‑1 (MQ=255)
                    cAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGga                                                                      <  2:229826/93‑1 (MQ=255)
                    cAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGga                                                                      <  2:35688/93‑1 (MQ=255)
                     aGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGgag                                                                     >  1:2657306/1‑93 (MQ=255)
                     aGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGGGgag                                                                     >  1:5876869/1‑93 (MQ=255)
                      gCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGgaga                                                                    <  2:3890832/93‑1 (MQ=255)
                      gCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGgaga                                                                    >  1:2212963/1‑93 (MQ=255)
                      gCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGgaga                                                                    >  2:1550052/1‑93 (MQ=255)
                       cTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGaa                                                                   <  1:1568340/93‑1 (MQ=255)
                         tCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAgg                                                                 <  1:635273/93‑1 (MQ=255)
                         tCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAgg                                                                 <  1:5774712/93‑1 (MQ=255)
                           gACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGc                                                               <  2:286512/93‑1 (MQ=255)
                           gACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGc                                                               >  2:4580557/1‑93 (MQ=255)
                           gACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGc                                                               >  1:974598/1‑93 (MQ=255)
                           gACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGc                                                               >  2:5493287/1‑93 (MQ=255)
                           gACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGc                                                               >  1:6635061/1‑93 (MQ=255)
                            aCAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCa                                                              >  1:2794256/1‑93 (MQ=255)
                                gCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGAcc                                                          <  2:805153/93‑1 (MQ=255)
                                   aTTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCa                                                       >  2:2921679/1‑93 (MQ=255)
                                   aTTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCa                                                       >  1:600389/1‑93 (MQ=255)
                                    ttCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCaa                                                      <  2:6768203/93‑1 (MQ=255)
                                     tCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAAt                                                     <  2:5575068/93‑1 (MQ=255)
                                     tCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAAt                                                     >  1:3645003/1‑93 (MQ=255)
                                        gtgGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaa                                                  <  1:5932148/93‑1 (MQ=255)
                                        gtgGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaa                                                  <  2:1415877/93‑1 (MQ=255)
                                         tgGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaaa                                                 <  1:2828670/93‑1 (MQ=255)
                                          ggCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCaaaa                                                <  2:3989805/93‑1 (MQ=255)
                                           gCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAAt                                               >  2:6108537/1‑93 (MQ=255)
                                             aGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATgg                                             <  2:1025803/93‑1 (MQ=255)
                                             aGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATgg                                             <  2:2455968/93‑1 (MQ=255)
                                              ggATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                                            >  1:489859/1‑93 (MQ=255)
                                              ggATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                                            <  2:4924607/93‑1 (MQ=255)
                                              ggATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGc                                            <  1:1749879/93‑1 (MQ=255)
                                               gATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGcc                                           >  2:2266040/1‑93 (MQ=255)
                                               gATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGcc                                           >  2:3051001/1‑93 (MQ=255)
                                                 tctcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGCCCGCAATAAAAATGGCCGa                                         <  1:7354248/93‑1 (MQ=255)
                                                 tctcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGa                                         >  2:5699049/1‑93 (MQ=255)
                                                   tcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGc                                       <  2:3863602/93‑1 (MQ=255)
                                                   tcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGc                                       <  2:3679259/93‑1 (MQ=255)
                                                   tcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGc                                       <  2:6159606/93‑1 (MQ=255)
                                                   tcGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGc                                       >  1:7466596/1‑93 (MQ=255)
                                                     gATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGt                                     <  2:4824612/93‑1 (MQ=255)
                                                     gATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGt                                     >  1:3204909/1‑93 (MQ=255)
                                                      aTCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTg                                    <  1:4667760/93‑1 (MQ=255)
                                                      aTCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTg                                    >  2:2997410/1‑93 (MQ=255)
                                                      aTCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTg                                    <  1:2311969/93‑1 (MQ=255)
                                                      aTCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTg                                    >  1:6692908/1‑93 (MQ=255)
                                                        tgcgTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg                                  <  1:288752/92‑1 (MQ=255)
                                                        cgcgTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTggg                                  >  1:5191969/1‑93 (MQ=255)
                                                          cgTGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCa                                >  2:7357227/1‑93 (MQ=255)
                                                            tGCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATc                              <  2:5615945/93‑1 (MQ=255)
                                                             gCTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCa                             <  2:5115789/93‑1 (MQ=255)
                                                              cTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAg                            >  2:6313012/1‑93 (MQ=255)
                                                              cTGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAg                            >  2:1235391/1‑93 (MQ=255)
                                                               tGGCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAgg                           <  1:4421643/93‑1 (MQ=255)
                                                                 gCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTc                         >  2:4161996/1‑93 (MQ=255)
                                                                 gCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTc                         <  2:2650/93‑1 (MQ=255)
                                                                 gCGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTc                         <  2:692138/93‑1 (MQ=255)
                                                                  cGGAAATGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCa                        <  2:4110230/93‑1 (MQ=255)
                                                                     aaaTGGGCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCAttt                     <  1:3224817/93‑1 (MQ=255)
                                                                         gggCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATg                 <  1:2577033/93‑1 (MQ=255)
                                                                          ggCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGc                >  2:2582117/1‑93 (MQ=255)
                                                                           gCCCGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCg               >  1:6914976/1‑93 (MQ=255)
                                                                            cccGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGa              <  1:1096838/93‑1 (MQ=255)
                                                                             ccGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGAt             >  2:2538948/1‑93 (MQ=255)
                                                                             ccGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGAt             <  1:1147567/93‑1 (MQ=255)
                                                                              cGCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATc            <  1:3759989/93‑1 (MQ=255)
                                                                               gCAAAAACTGCGTGTGCTGGATGCTGGCGGGGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCt           >  1:6182640/1‑93 (MQ=255)
                                                                               gCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCt           <  1:4214175/93‑1 (MQ=255)
                                                                               gCAAAAACTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCt           <  1:5255912/93‑1 (MQ=255)
                                                                                   aaaCTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCt       <  1:1184062/93‑1 (MQ=255)
                                                                                    aaCTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTg      >  1:5586064/1‑93 (MQ=255)
                                                                                     aCTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTgc     <  2:4223680/93‑1 (MQ=255)
                                                                                     aCTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTgc     <  2:7442452/93‑1 (MQ=255)
                                                                                     aCTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTgc     <  1:3064595/93‑1 (MQ=255)
                                                                                     aCTGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTgc     <  2:5009820/93‑1 (MQ=255)
                                                                                       tGCGCGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTgcgc   <  2:5341358/93‑1 (MQ=255)
                                                                                        gcgcGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCa  >  2:665495/1‑93 (MQ=255)
                                                                                           |                                                                                         
TTTACGGCACCACCAAAGGGCAGCTTCGACAGGCTATTCTGTGGCAGGATCTCGATCGCGTGCTGGCGGAAATGGGCCCGCAAAAACTGCGTGTGCTGGATGCTGGCGGTGGAGAAGGGCAGACCGCAATCAAAATGGCCGAGCGTGGGCATCAGGTCATTTTATGCGATCTTTCTGCGCA  >  NC_000913/972809‑972989

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