Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,040,218 G→A A127T (GCA→ACA)  appA → phosphoanhydride phosphorylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,040,2180GA100.0% 762.0 / NA 208A127T (GCA→ACA) appAphosphoanhydride phosphorylase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (111/97);  total (111/97)

Only 100 of 220 total aligned reads displayed.

AGGGCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGGCACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAA  >  NC_000913/1040129‑1040307
                                                                                         |                                                                                         
aGGGCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGacac                                                                                        >  2:39638/1‑93 (MQ=255)
aGGCCTCCCCGAAGTCTGGTCAGGTCGCATTTATTGCTGATGTCGACGAGCGTCCCCGAAAAACAGGCGAAGCCTTCGCCGCCGGGCTGacac                                                                                        <  2:3786740/93‑1 (MQ=255)
 gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTCCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       <  2:5475706/93‑1 (MQ=255)
 gggCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACAcc                                                                                       >  2:2258386/1‑93 (MQ=255)
      gCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACt                                                                                  >  2:300504/1‑93 (MQ=255)
       cccGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtg                                                                                 <  1:4660013/93‑1 (MQ=255)
        ccGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtgt                                                                                <  1:2124297/93‑1 (MQ=255)
        ccGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtgt                                                                                >  2:5362370/1‑93 (MQ=255)
        ccGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTAACCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtgt                                                                                >  2:372574/1‑93 (MQ=255)
         cGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACtgtg                                                                               >  1:5984692/1‑93 (MQ=255)
          gCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGc                                                                              >  2:2027984/1‑93 (MQ=255)
          gCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGACTTCGCCGCCGGGCTGACACCTGACTGTGc                                                                              >  1:3786975/1‑93 (MQ=255)
           cAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCa                                                                             >  1:4507430/1‑93 (MQ=255)
           cAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCa                                                                             >  2:2601467/1‑93 (MQ=255)
             gTCTGGTCAGTTCGCGATTATTGCTGATGTCGACGAGCGTCCCCGAAAAATAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           <  1:4434334/93‑1 (MQ=255)
             gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  1:2361877/1‑93 (MQ=255)
             gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  2:1123114/1‑93 (MQ=255)
             gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           <  2:4853955/93‑1 (MQ=255)
             gTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAAt                                                                           >  1:498194/1‑93 (MQ=255)
               cgggTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATaa                                                                         <  2:4403573/91‑1 (MQ=255)
                tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        <  1:1046449/93‑1 (MQ=255)
                tGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAc                                                                        <  1:1541025/93‑1 (MQ=255)
                 ggTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAAcc                                                                       >  2:813194/1‑93 (MQ=255)
                  gTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCg                                                                      >  2:5251494/1‑93 (MQ=255)
                   tCAGGTCTCGATTATTTCTGATGTCGACGAGCGTACCCTTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGt                                                                     <  2:199122/93‑1 (MQ=255)
                     aGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTAc                                                                   >  2:4249428/1‑93 (MQ=255)
                      ggTCGCGATTATTGCTGAGGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACa                                                                  <  1:555185/93‑1 (MQ=255)
                       gacgcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACAt                                                                 >  2:3115371/3‑93 (MQ=255)
                       gTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAACCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACAt                                                                 <  1:2425547/93‑1 (MQ=255)
                       gTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACAt                                                                 >  1:4039534/1‑93 (MQ=255)
                        tCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATa                                                                >  1:3821433/1‑93 (MQ=255)
                         cgcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAc                                                               >  1:5388147/1‑93 (MQ=255)
                          gcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAcc                                                              >  2:2139011/1‑93 (MQ=255)
                          gcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAcc                                                              >  1:1608844/1‑93 (MQ=255)
                          gcgATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAAGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATAcc                                                              >  2:89104/1‑93 (MQ=255)
                            gATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCa                                                            <  2:1533855/93‑1 (MQ=255)
                             attattGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAg                                                           >  1:4872159/1‑93 (MQ=255)
                             attattGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAg                                                           >  2:2699983/1‑93 (MQ=255)
                             attattGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAg                                                           <  1:3424317/93‑1 (MQ=255)
                                  tGCTGATGTCGACGAGCGTCCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  1:5477806/93‑1 (MQ=255)
                                  tGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  2:1503976/93‑1 (MQ=255)
                                  tGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGa                                                      <  2:5618394/93‑1 (MQ=255)
                                    cTGATGTCGACGAGCTTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATa                                                    <  2:3205609/93‑1 (MQ=255)
                                    cTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATa                                                    <  2:5594386/93‑1 (MQ=255)
                                      gATGTCGAGGAGCGTACCCGTAAAGCAGGCGAAGCCTTCGCCGCGGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGACAGATACg                                                  <  2:6118537/93‑1 (MQ=255)
                                         gTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTcc                                               <  2:689393/93‑1 (MQ=255)
                                          tCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCa                                              <  1:5167608/93‑1 (MQ=255)
                                             acgaGCGTCCCCGAAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTc                                           <  2:286488/93‑1 (MQ=255)
                                             acgaGCGTATCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTc                                           <  2:3902104/93‑1 (MQ=255)
                                             acgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTc                                           >  2:2435692/1‑93 (MQ=255)
                                              cgaGCGTACCGGAAAAACGGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTcc                                          <  2:3424082/93‑1 (MQ=255)
                                              cgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTcc                                          >  1:6076488/1‑93 (MQ=255)
                                              cgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTcc                                          <  2:1737037/93‑1 (MQ=255)
                                              cgaGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTACAGTcc                                          >  2:1398371/1‑93 (MQ=255)
                                                aGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCg                                        >  1:4322562/1‑93 (MQ=255)
                                                aGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCg                                        >  1:4780528/1‑93 (MQ=255)
                                                aGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCg                                        <  2:2540663/93‑1 (MQ=255)
                                                 gCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGa                                       <  1:6246641/93‑1 (MQ=255)
                                                 gCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGa                                       >  1:2744483/1‑93 (MQ=255)
                                                   gTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATc                                     >  1:4298428/1‑93 (MQ=255)
                                                    tccccGAAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATcc                                    <  1:64230/91‑1 (MQ=255)
                                                     aCCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCg                                   >  2:843429/1‑93 (MQ=255)
                                                     aCCCGAAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCg                                   <  1:5335258/93‑1 (MQ=255)
                                                      cccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGt                                  >  1:2713237/1‑93 (MQ=255)
                                                      cccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGt                                  <  2:5737154/93‑1 (MQ=255)
                                                       ccGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGtt                                 <  2:4658515/93‑1 (MQ=255)
                                                       cagTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGtt                                 <  1:6118068/91‑1 (MQ=255)
                                                        cGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTa                                <  2:1087907/93‑1 (MQ=255)
                                                         gTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAGTAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAt                               >  2:2083784/1‑93 (MQ=255)
                                                         gTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAt                               >  2:1667824/1‑93 (MQ=255)
                                                          tAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAtt                              >  2:4354699/1‑93 (MQ=255)
                                                          tAAAACAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAtt                              <  1:3613584/93‑1 (MQ=255)
                                                           aaaaCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAttt                             <  2:4896952/93‑1 (MQ=255)
                                                           aaaaCAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTAttt                             <  2:1825503/93‑1 (MQ=255)
                                                               cAGGCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATc                         >  1:5054064/1‑93 (MQ=255)
                                                                 ggCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCct                       >  1:2967705/1‑93 (MQ=255)
                                                                 ggCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCct                       >  1:1452589/1‑93 (MQ=255)
                                                                 ggCGAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCct                       >  1:773515/1‑93 (MQ=255)
                                                                    gAAGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTa                    >  1:950612/1‑93 (MQ=255)
                                                                      aGCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaa                  <  1:1236794/93‑1 (MQ=255)
                                                                       gCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaa                 >  1:2384973/1‑93 (MQ=255)
                                                                       gCCTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaa                 >  1:1951165/1‑93 (MQ=255)
                                                                        ccTTCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTaaaaa                >  1:5364778/1‑93 (MQ=255)
                                                                          ttCGCCGCCGGGCTGACACCTGACTGTGTAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACt              >  1:3898848/1‑93 (MQ=255)
                                                                          ttCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACt              >  1:3014509/1‑93 (MQ=255)
                                                                          ttCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACt              >  1:1366947/1‑93 (MQ=255)
                                                                           tCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTg             <  1:2836970/93‑1 (MQ=255)
                                                                           tCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTg             <  2:2831564/93‑1 (MQ=255)
                                                                           tCGCCGCCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTg             >  1:2308789/1‑93 (MQ=255)
                                                                            cgccgcCGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTgg            <  2:4588343/93‑1 (MQ=255)
                                                                             gccgccGGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGc           >  2:5947280/1‑93 (MQ=255)
                                                                              ccgccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCg          >  1:1097025/1‑93 (MQ=255)
                                                                               cgccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGt         >  2:406605/1‑93 (MQ=255)
                                                                                gccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGtt        <  2:158312/93‑1 (MQ=255)
                                                                                 ccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGttt       <  2:3250882/93‑1 (MQ=255)
                                                                                 ccgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGttt       >  1:5920029/1‑93 (MQ=255)
                                                                                  cgGGCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTg      <  2:1695059/93‑1 (MQ=255)
                                                                                   gggCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGc     >  2:2217049/1‑93 (MQ=255)
                                                                                    ggCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGcc    >  1:2214773/1‑93 (MQ=255)
                                                                                     gCTGACACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCa   >  2:5220644/1‑93 (MQ=255)
                                                                                         |                                                                                         
AGGGCTGCCCGCAGTCTGGTCAGGTCGCGATTATTGCTGATGTCGACGAGCGTACCCGTAAAACAGGCGAAGCCTTCGCCGCCGGGCTGGCACCTGACTGTGCAATAACCGTACATACCCAGGCAGATACGTCCAGTCCCGATCCGTTATTTAATCCTCTAAAAACTGGCGTTTGCCAA  >  NC_000913/1040129‑1040307

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