Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,734,835 G→C 7.9% T295T (ACC→ACG hutH ← histidine ammonia‑lyase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,734,8350GC 7.9% 199.9 / 7.5 63T295T (ACC→ACGhutHhistidine ammonia‑lyase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (3/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

CAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTT  >  NC_002947/5734688‑5734980
                                                                                                                                                   |                                                                                                                                                 
cagggccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGt                                                                                                                                                 >  2:496116/1‑150 (MQ=255)
    gccaggTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTcagg                                                                                                                                             <  2:234178/150‑1 (MQ=255)
              cggtggCCATGGCCGCCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCa                                                                                                                                   >  2:401803/1‑150 (MQ=255)
               ggcggcCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCAt                                                                                                                                  <  1:2240/150‑1 (MQ=255)
                  ggcCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGGAGGCACCCATGAc                                                                                                                               <  1:1111620/150‑1 (MQ=255)
                            aCCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTgg                                                                                                                     <  1:672366/150‑1 (MQ=255)
                             ccGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCAGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgc                                                                                                                    >  1:1056027/1‑150 (MQ=255)
                                gTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagc                                                                                                                 >  2:729073/1‑150 (MQ=255)
                                     gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgcagc                                                                                                            <  2:953086/150‑1 (MQ=255)
                                     gCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcgcagc                                                                                                            >  2:1150420/1‑150 (MQ=255)
                                        tGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCATGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAg                                                                                                         >  1:666415/1‑150 (MQ=255)
                                          gAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTa                                                                                                       >  2:907214/1‑150 (MQ=255)
                                             gtcgccgccGGAGACCACATCACCTTCGGTAGCGAACACCAGCGCGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCACCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                                    <  2:247302/147‑1 (MQ=255)
                                             gTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                                    <  2:629751/150‑1 (MQ=255)
                                                   gccgGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTg                                                                                              <  1:401803/150‑1 (MQ=255)
                                                         gaTCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCtt                                                                                        <  1:464465/150‑1 (MQ=255)
                                                          atcacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTg                                                                                       >  1:383143/1‑150 (MQ=255)
                                                           tcacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGt                                                                                      <  1:729073/150‑1 (MQ=255)
                                                            cacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTc                                                                                     <  1:359246/150‑1 (MQ=255)
                                                            cacatcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTc                                                                                     <  1:496116/150‑1 (MQ=255)
                                                              catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                                    <  1:558889/133‑1 (MQ=255)
                                                              catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACgg                                                                                                    >  2:558889/1‑133 (MQ=255)
                                                              catcacCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGc                                                                                   >  1:770915/1‑150 (MQ=255)
                                                                    cTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCacaaga                                                                             >  2:350531/1‑145 (MQ=255)
                                                                     ttCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCtt                                                                            <  2:666415/150‑1 (MQ=255)
                                                                     ttCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCACTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCtt                                                                            <  1:1150420/150‑1 (MQ=255)
                                                                        ggCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg                                                                                                                <  1:651462/111‑1 (MQ=255)
                                                                        ggCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGgcagcg                                                                                                                >  2:651462/1‑111 (MQ=255)
                                                                         gCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgc                                                                        <  1:907214/150‑1 (MQ=255)
                                                                           aGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcga                                                                      >  1:250496/1‑150 (MQ=255)
                                                                                   ccAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgcgacagcgaca                                                              <  2:770915/150‑1 (MQ=255)
                                                                                    cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTca                                                                                                                                               <  2:95331/68‑1 (MQ=255)
                                                                                    cAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTca                                                                                                                                               >  1:95331/1‑68 (MQ=255)
                                                                                      gCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACtt                                                           <  1:26711/150‑1 (MQ=255)
                                                                                         tcaagTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGt                                                                                      <  1:39664/116‑1 (MQ=255)
                                                                                         tcaagTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGt                                                                                      >  2:39664/5‑120 (MQ=255)
                                                                                           ttGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                      <  1:307089/150‑1 (MQ=255)
                                                                                             gTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGa                                                    <  2:660574/150‑1 (MQ=255)
                                                                                              tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCa                                                                                  >  2:655183/1‑119 (MQ=255)
                                                                                              tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCa                                                                                  <  1:655183/119‑1 (MQ=255)
                                                                                              tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                       <  2:651890/146‑1 (MQ=255)
                                                                                              tCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCt                                                       >  1:651890/1‑146 (MQ=255)
                                                                                                    acaGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgc                                                                        <  2:69614/123‑1 (MQ=255)
                                                                                                    acaGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTgc                                                                        >  1:69614/1‑123 (MQ=255)
                                                                                                    aaaGCGTTGGCTTCGCTGCACAGTAACTCGGCAGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCAGCTGGCAGCGCAGCGAGTACAGGTCATGGACCTTGTAGCCGTTCTTGTGCGACCGAGACCCTTCGCTGGAGTCGccc                                             >  2:850247/3‑150 (MQ=255)
                                                                                                          ttGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGATCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCCgg                                       >  2:97692/1‑150 (MQ=255)
                                                                                                             gCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcg                                    >  1:1049905/1‑150 (MQ=255)
                                                                                                              cTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                      <  2:81628/131‑1 (MQ=255)
                                                                                                              cTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTg                                                      >  1:81628/1‑131 (MQ=255)
                                                                                                               ttCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTcgcg                                  >  1:701267/1‑150 (MQ=255)
                                                                                                                      cccAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCATGGGTCGCCCAGCAGGTCGCGGAAGCAg                           <  1:823576/150‑1 (MQ=255)
                                                                                                                       ccAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCa                                                                                  <  2:361090/94‑1 (MQ=255)
                                                                                                                       ccAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCa                                                                                  >  1:361090/1‑94 (MQ=255)
                                                                                                                       ccAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAgg                          >  1:392381/1‑150 (MQ=255)
                                                                                                                              cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGc                                               <  2:240009/122‑1 (MQ=255)
                                                                                                                              cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGc                                               >  1:240009/1‑122 (MQ=255)
                                                                                                                              cTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGgt                   <  1:273044/150‑1 (MQ=255)
                                                                                                                                  cctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCg               >  2:348565/4‑150 (MQ=255)
                                                                                                                                   ctgcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGa              >  2:768820/3‑150 (MQ=255)
                                                                                                                                     gcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGaa                               >  2:258359/1‑131 (MQ=255)
                                                                                                                                     gcCTGACCCAGCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGaa                               <  1:258359/131‑1 (MQ=255)
                                                                                                                                            gCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCcgc     >  1:868629/1‑150 (MQ=255)
                                                                                                                                               gCTGCGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGtt  <  1:768820/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CAGGGCCAGGTTGTCGGCGGCCATGGCCACCGGTTCGGCGTGGAAGTTGCCGCCGGAGATCACATCACCTTCGGCAGCGAACACCAGCGGGTTGTCCGACACAGCGTTGGCTTCGATGCCCAGTACCTCGGCGGCCTGACGCAGCTGGGTCAGGCAGGCACCCATGACCTGCGGCTGGCAGCGCAGCGAGTACGGGTCCTGGACCTTGTCGCAGTTCTTGTGCGACAGCGACACTTCGCTGGAGTCGCCCAGCAGGTCGCGGAAGCAGGCCGCGGTGTCGATCTGGCCGCGTT  >  NC_002947/5734688‑5734980

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