Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,909,296 A→T 6.4% intergenic (‑76/+93) spuD ← / ← spuC‑II spermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,909,2960AT 6.4% 251.4 / 3.4 87intergenic (‑76/+93)spuD/spuC‑IIspermidine/putrescine ABC transporter substrate‑binding protein/polyamine:pyruvate transaminase
Reads supporting (aligned to +/- strand):  ref base A (39/42);  new base T (4/2);  total (43/44)
Fisher's exact test for biased strand distribution p-value = 4.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909149‑5909443
                                                                                                                                                       |                                                                                                                                                    
ttCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGaaa                                                                                                                                                     >  6:300130/1‑149 (MQ=255)
ttCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGaa                                                                                                                                                      >  7:235152/1‑148 (MQ=255)
      cTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTc                                                                                                                                                >  1:261048/1‑149 (MQ=255)
       tGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACa                                                                                                                                                   >  2:224452/1‑144 (MQ=255)
       tGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACa                                                                                                                                                   <  1:224452/144‑1 (MQ=255)
        gAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGc                                                                                                                                              >  6:160829/1‑148 (MQ=255)
         aaCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGc                                                                                                                                              >  2:263456/1‑147 (MQ=255)
         aaCAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGc                                                                                                                                              <  7:335315/147‑1 (MQ=255)
               ggAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGC‑‑CTGAAACAGTCGCTAGTCTgg                                                                                                                                      >  1:111396/1‑149 (MQ=255)
                  ggCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                   <  6:194585/149‑1 (MQ=255)
                  ggCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCaa                                                                                                                                   <  1:263456/149‑1 (MQ=255)
                     cATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGt                                                                                                                                 >  5:34874/1‑148 (MQ=255)
                      aTGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTa                                                                                                                               <  5:300130/149‑1 (MQ=255)
                        ggCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTAc                                                                                                                              >  7:271050/1‑148 (MQ=255)
                         gCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACaa                                                                                                                            <  5:52638/149‑1 (MQ=255)
                          cACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAg                                                                                                                           <  8:235152/149‑1 (MQ=255)
                           aCCCATCAGGGCGGCGGCCAGGAACGTGTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAgg                                                                                                                          >  1:224129/1‑149 (MQ=255)
                             ccATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGcc                                                                                                                        >  8:118792/1‑149 (MQ=255)
                             ccATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGcc                                                                                                                        >  3:214690/1‑149 (MQ=255)
                              cATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGccc                                                                                                                       >  3:16335/1‑149 (MQ=255)
                              cATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGccc                                                                                                                       <  7:316170/149‑1 (MQ=255)
                                tCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTATTTCATGCGTTATGATCGAGAATTTTTATTAAGCAACTGTCGAGCAGCCAGATCCCAAGGTGCAGAGCGCGCGGCGA‑‑CTGAATCAGTCTCTAGACGGGCAAGTTACAAggacgtg                                                                                                                     >  8:331068/1‑144 (MQ=255)
                                    ggcggcggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCt                                                                                                                      >  8:332118/1‑144 (MQ=255)
                                    ggcggcggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCt                                                                                                                      <  7:332118/144‑1 (MQ=255)
                                      cggcggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCaa                                                                                                               <  2:224129/149‑1 (MQ=255)
                                          ggcCAGCAACGTCTTGCCCATTTTATTCCTGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCAAGCAGCCAGAACCCACGGTGCAGAGAGCGCGGCGA‑‑CGGTTAAAGTCGCTAGTCTGGCCAGTTACACGGCCGTGGATCAAGGa                                                                                                            >  6:345672/1‑148 (MQ=255)
                                          ggcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGaaa                                                                                                          <  4:163517/149‑1 (MQ=255)
                                           gcCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGaaa                                                                                                          <  3:85968/149‑1 (MQ=255)
                                             cAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGATACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTa                                                                                                                   <  8:151189/138‑1 (MQ=255)
                                             cAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAg                                                                                                                                  <  8:328606/123‑1 (MQ=255)
                                             cAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAg                                                                                                                                  >  7:328606/1‑123 (MQ=255)
                                             cAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTa                                                                                                                   >  7:151189/1‑138 (MQ=255)
                                             cAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGc                                                                                                        <  5:160829/149‑1 (MQ=255)
                                                  aaCGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAgg                                                                                                     <  8:84399/148‑1 (MQ=255)
                                                  aaCGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAgg                                                                                                     <  3:232890/148‑1 (MQ=255)
                                                   aCGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGc                                                                                                   >  1:66844/1‑149 (MQ=255)
                                                    cGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTg                                                                                                                                       <  6:234875/112‑1 (MQ=255)
                                                    cGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGcc                                                                                                  <  8:271050/149‑1 (MQ=255)
                                                      tCTTGACCATTTTATTCATGCGTAATGCTACAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCAt                                                                                                >  6:333435/1‑149 (MQ=255)
                                                           cccATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATATAAc                                                                                           <  2:111396/149‑1 (MQ=255)
                                                             cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGa                                                                                         <  7:118792/149‑1 (MQ=255)
                                                             cATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGa                                                                                         <  2:180755/149‑1 (MQ=255)
                                                                        tGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGGAGGCTGAAACAGTCGCTAGTCTGGCAAGTTAAAAGGCCCTGTATCAAGAAAAGCAGGGCCATGTAACGACTTAATGt                                                                                 >  1:299969/1‑149 (MQ=255)
                                                                             aaTGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGaaa                                                                                                          >  6:8783/1‑116 (MQ=255)
                                                                             aaTGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGaaa                                                                                                          <  5:8783/116‑1 (MQ=255)
                                                                               tGCTCCAGAATTTTTATTTAGCAACTCTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACACTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGcta                                                                       <  5:90834/149‑1 (MQ=255)
                                                                                gCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctag                                                                      >  5:309360/1‑149 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                    >  2:288058/1‑149 (MQ=255)
                                                                                  tCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                    >  8:92222/1‑149 (MQ=255)
                                                                                   ccAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcc                                                                    >  3:364086/1‑148 (MQ=255)
                                                                                    cAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGTCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagc                                                                     <  5:333435/146‑1 (MQ=255)
                                                                                      gAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGctagcctagc                                                                <  1:306147/149‑1 (MQ=255)
                                                                                         tttttATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCaaa                                                             >  1:121558/1‑149 (MQ=255)
                                                                                           tttatttaGCAACTTTCCAGCAGCCAGAACCCACTGTTCAGAGCGCGCGGCGA‑‑CTGAAACATTAGGTAGTCTGGCAATTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAAc                                                            <  5:345672/148‑1 (MQ=255)
                                                                                             tatttaGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCt                                                         >  7:231301/1‑149 (MQ=255)
                                                                                                 taGCAACTGTCCAGCAGACAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAACGTCACCACGCTAGCCTAGCAAACGCTTAtt                                                     >  2:188045/1‑149 (MQ=255)
                                                                                                      aCTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                                >  4:59623/1‑149 (MQ=255)
                                                                                                             agcagcCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCaa                                         <  7:92222/149‑1 (MQ=255)
                                                                                                              gcagcCAGAACACACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAg                                        <  3:59623/149‑1 (MQ=255)
                                                                                                               cagcCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAgg                                       >  2:230686/1‑149 (MQ=255)
                                                                                                                agcCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGCCCCTGTATCAAGGAAAGCAGGGCCTTGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAggg                                      >  5:296696/1‑149 (MQ=255)
                                                                                                                  ccAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTc                                    <  8:231301/149‑1 (MQ=255)
                                                                                                                    aGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAg                                 <  6:309360/149‑1 (MQ=255)
                                                                                                                         ccACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTc                              >  5:168421/1‑147 (MQ=255)
                                                                                                                          cACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAg                           <  4:364086/149‑1 (MQ=255)
                                                                                                                                 cAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAACGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTg                    <  4:214690/149‑1 (MQ=255)
                                                                                                                                  agagCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgc                   >  7:1946/1‑149 (MQ=255)
                                                                                                                                  agagCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGACTCAAGGGTCAGGTCCAGGCACTTgc                   >  2:184844/1‑149 (MQ=255)
                                                                                                                                    agCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgc                 >  7:82856/1‑149 (MQ=255)
                                                                                                                                      agactagcgA‑‑CTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTgcgcgc               >  1:355043/7‑149 (MQ=255)
                                                                                                                                       gactagcgA‑‑CTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGc                         >  5:2483/6‑138 (MQ=255)
                                                                                                                                       gactagcgA‑‑CTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGc                         <  6:2483/133‑1 (MQ=255)
                                                                                                                                         gcggcgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCtt            >  5:43570/1‑149 (MQ=255)
                                                                                                                                          tagcgA‑‑CTGTTTCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCttt           >  7:349602/3‑149 (MQ=255)
                                                                                                                                           ggcgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                                >  7:267040/1‑112 (MQ=255)
                                                                                                                                           ggcgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATc                                                <  8:267040/112‑1 (MQ=255)
                                                                                                                                             cgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCATGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         <  5:375267/148‑1 (MQ=255)
                                                                                                                                             cgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCATGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  6:375267/1‑148 (MQ=255)
                                                                                                                                             cgA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTcc        <  2:299969/149‑1 (MQ=255)
                                                                                                                                                a‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      <  1:288058/149‑1 (MQ=255)
                                                                                                                                                a‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      <  4:16335/149‑1 (MQ=255)
                                                                                                                                                a‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      <  6:43570/149‑1 (MQ=255)
                                                                                                                                                a‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      <  2:121558/149‑1 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTc                                    >  3:316410/1‑118 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTc                                    <  4:316410/118‑1 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      >  1:160192/1‑147 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTccac      <  2:160192/147‑1 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCAGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         >  8:297442/1‑144 (MQ=255)
                                                                                                                                                    tGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCAGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTc         <  7:297442/144‑1 (MQ=255)
                                                                                                                                                      aaaCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGt  >  3:46751/1‑149 (MQ=255)
                                                                                                                                                      aaaCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGt  >  7:207168/1‑149 (MQ=255)
                                                                                                                                                       aaCAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCTAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGt  <  1:230686/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TTCAGCCTGAACAGCGGAGGCCATGGCACCCATCAGGGCGGCGGCCAGCAACGTCTTGCCCATTTTATTCATGCGTAATGCTCCAGAATTTTTATTTAGCAACTGTCCAGCAGCCAGAACCCACGGTGCAGAGCGCGCGGCGA‑‑CTGAAACAGTCGCTAGTCTGGCAAGTTACAAGGCCCTGTATCAAGGAAAGCAGGGCCATGTAACGACTTAATGTCACAACGCTAGCCTAGCAAACGCTTATTGAATCGCCTCAAGGGTCAGGTCCAGGCACTTGCGCGCCTTTTCCACCAGT  >  NC_002947/5909149‑5909443

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