Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,125,244 T→A 7.3% T435S (ACC→TCC) ‡ tdcG‑II ← L‑serine dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,125,2440TA 7.3% 185.2 / 3.6 67T435S (ACC→TCC) ‡tdcG‑IIL‑serine dehydratase
Reads supporting (aligned to +/- strand):  ref base T (31/31);  new base A (2/3);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01

GGGCATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATCGGG  >  NC_002947/1125106‑1125389
                                                                                                                                          |                                                                                                                                                   
gggCATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCAGACCGgcgccga                                                                                                                                           >  8:32122/1‑146 (MQ=255)
  gCATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGt                                                                                                                                         >  6:177065/1‑149 (MQ=255)
  gCATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCAGACCGgcgccgaca                                                                                                                                         <  7:32122/149‑6 (MQ=255)
   cATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCAGACCGgcgccgacat                                                                                                                                        <  1:301670/149‑7 (MQ=255)
              gTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGgcgc                                                                                                                                              <  4:210194/132‑1 (MQ=255)
              gTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGgcgc                                                                                                                                              >  3:210194/1‑132 (MQ=255)
               tCTTGAGCAGTGTTTCGGACATCGATGACTCCTAGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTg                                                                                                                            >  1:203791/1‑149 (MQ=255)
                     gCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATgg                                                                                                                                          <  7:139002/129‑1 (MQ=255)
                     gCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATgg                                                                                                                                          >  8:139002/1‑129 (MQ=255)
                     gCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATgg                                                                                                                                          <  3:252686/129‑1 (MQ=255)
                     gCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGGTGACCTTGTCCAG‑c                                                                                                                     <  5:232829/149‑1 (MQ=255)
                     gCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGACGGTCTGGCGCATgg                                                                                                                                          >  4:252686/1‑129 (MQ=255)
                        gtgtTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGac                                                                                                                  <  6:8860/149‑1 (MQ=255)
                         tgtTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGaca                                                                                                                 >  5:192238/1‑149 (MQ=255)
                             tCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTa                                                                                                             >  1:114012/1‑149 (MQ=255)
                               ggACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAg                                                                                                            <  6:62437/148‑1 (MQ=255)
                                    tCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGtgctgc                                                                                                      <  3:279894/149‑1 (MQ=255)
                                      gATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTAGGCGCCGGTCTGCCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGccc                                                                                                    >  5:69714/1‑149 (MQ=255)
                                       aTGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCg                                                                                                   >  1:63811/1‑149 (MQ=255)
                                           cTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTcgc                                                                                               >  6:168348/1‑149 (MQ=255)
                                                 gCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCcgcg                                                                                           <  3:235403/147‑1 (MQ=255)
                                                       aTTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCa                                                                                   <  5:122017/149‑1 (MQ=255)
                                                              ggTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTAc                                                                            <  5:177065/149‑1 (MQ=255)
                                                                tCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGc                                                                                <  1:6008/143‑1 (MQ=255)
                                                                tCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGc                                                                                >  2:6008/1‑143 (MQ=255)
                                                                   aCATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGtgctgc                                                                                                      <  3:186186/118‑1 (MQ=255)
                                                                   aCATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGtgctgc                                                                                                      >  4:186186/1‑118 (MQ=255)
                                                                    cATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGtt                                                                      <  6:69714/149‑1 (MQ=255)
                                                                     aTTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAGCCGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCAt                                                                                  <  1:49030/137‑1 (MQ=255)
                                                                     aTTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCAt                                                                                  >  2:49030/1‑136 (MQ=255)
                                                                        cGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGAt                                                                   <  3:144667/148‑1 (MQ=255)
                                                                         gatgatGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATgg                                                                 >  2:108530/1‑149 (MQ=255)
                                                                         gatgatGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATgg                                                                 >  4:90301/1‑149 (MQ=255)
                                                                           tgatgTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGcc                                                               <  2:234510/149‑1 (MQ=255)
                                                                              tgTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTc                                                            >  2:76177/1‑149 (MQ=255)
                                                                                ttGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑c                                                         >  1:206190/1‑149 (MQ=255)
                                                                                  gACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCg                                                       >  1:322146/1‑149 (MQ=255)
                                                                                  gACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCg                                                       >  5:44235/1‑149 (MQ=255)
                                                                                     ggCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACGGCGTTGATGGCCTTCA‑CCGAGc                                                    >  5:182473/1‑149 (MQ=255)
                                                                                     ggCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGc                                                    >  7:122993/1‑149 (MQ=255)
                                                                                     ggCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGc                                                    <  7:27047/149‑1 (MQ=255)
                                                                                         aGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCAt                                                <  2:322146/149‑1 (MQ=255)
                                                                                              gccACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCACCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCAt                                           >  2:295003/1‑149 (MQ=255)
                                                                                                       ggTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGa                                  <  6:192238/149‑1 (MQ=255)
                                                                                                          ctctTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATa                                >  2:224283/1‑148 (MQ=255)
                                                                                                               tGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAggg                           >  3:110218/1‑148 (MQ=255)
                                                                                                                   tttGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCt                      <  6:182473/149‑1 (MQ=255)
                                                                                                                   tttGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCt                      >  8:113390/1‑149 (MQ=255)
                                                                                                                    ttGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCt                      <  8:122993/148‑1 (MQ=255)
                                                                                                                    ttGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTg                     >  2:169222/1‑149 (MQ=255)
                                                                                                                    ttGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGACCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCG‑TACAggtgccctg                     >  4:51550/1‑141 (MQ=37)
                                                                                                                           tCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGTGCCCTGTatcga               <  3:90301/148‑5 (MQ=255)
                                                                                                                           tCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAg               <  1:169222/148‑1 (MQ=255)
                                                                                                                            cATGTCGGCACCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAgg              <  1:295003/148‑1 (MQ=255)
                                                                                                                             aTGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCGCCTGTACGAGgcc            >  7:267960/1‑149 (MQ=255)
                                                                                                                             aTGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGgcc            >  2:186837/1‑149 (MQ=255)
                                                                                                                             aTGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGgcc            >  2:191729/1‑149 (MQ=255)
                                                                                                                                tCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCACCCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGgccgc          <  1:76177/149‑1 (MQ=255)
                                                                                                                                tCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑GCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGgccgcc         <  2:63811/149‑1 (MQ=255)
                                                                                                                                tCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGgccgc          <  2:203791/148‑1 (MQ=255)
                                                                                                                                catgcgcCAGACCGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAg                                                     <  4:61933/102‑1 (MQ=25)
                                                                                                                                catgcgcCAGACCGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAg                                                     >  3:61933/4‑105 (MQ=25)
                                                                                                                                    cgcCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATc     <  1:191729/149‑1 (MQ=255)
                                                                                                                                    cgcCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATc     <  1:186837/149‑1 (MQ=255)
                                                                                                                                      ccGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATCgg   >  6:299710/1‑149 (MQ=255)
                                                                                                                                      ccGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATCgg   <  1:224283/149‑1 (MQ=255)
                                                                                                                                       cGGTCGGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATCggg  <  7:113390/149‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                                   
GGGCATGCAGTGGGGTCTTGAGCAGTGTTTCGGACATCGATGACTCCTTGCTGTCATTGTCGGGTCAACATTCGATGATGTTGACGGCCAGGCCGCCACGGGCGGTCTCTTTGTATTTGCTTTTCATGTCGGCGCCGGTCTGGCGCATGGTGCGGATGACCTTGTCCAG‑CGACACGTAGTGCTGCCCGTCGCCGCGCAGGGCCATGCGTACCGCGTTGATGGCCTTCA‑CCGAGCCCATCGCATTGCGCTCGATACAGGGCACCTGTACGAGGCCGCCGATCGGG  >  NC_002947/1125106‑1125389

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